; ^ 
! c 



! H 

I o 



LAW OFFICES 

SUGHRUE, MiON, ZiNN, MACPEAK & SeAS, PLLC 

2100 PENNSYLVANIA AVENUE, N,W, 
WASHINGTON, DX. 20037-3202 
TELEPHONE (202) 293-7060 
FACSIMILE (202) 293-7860 



CALIFORNIA OFFICE 

lOlOELCAMINO REAL 
MENLO PARK, CA 94025 
TELEPHONE (650) 325-5800 
FACSIMILE (650) 325-6606 

BOX: PATENT APPLICATION 

Assistant Commissioner for Patents 
Washington, D.C. 20231 



January 14, 2000 



xt-^ 



JAPAN OFFICE 

TOEI NISHI SHIMBASHI BLDG. 4F 
13-5 NISHI SHIMBASHI 1-CHOME 
MINATO-KU, TOKYO 105, JAPAN 
TELEPHONE (03) 3503-3760 
FACSIMILE (03) 3503-3756 



Re: Application of Masahiko YAMADA 
METHOD OF AND SYSTEM FOR 
CARRYING OUT THE METHOD 

Our Reference: Q56529 



IMAGE PROCESSING AND RECORDING MEDIUM FOR 



Dear Sir: 

Attached hereto is the application identified above including the specification, claims, executed 
declaration and Power of Attorney, sixty-one (61) sheets of drawings, four (4) priority documents. Information 
l^isclosure Statement and PTO Form 1449 with references, executed Assignment and PTO Form 1595, 

The Government filing fee is calculated as follows: 



Total Claims 
Independent Claims 
Base Filing Fee 
Multiple Dep. Claim Fee 
TOTAL FILING FEE 
Recordation of Assignment Fee 
TOTAL U.S. GOVERNMENT FEE 



99 - 20 = 
21 - 3 = 
($690.00) 
($260.00) 



79x $18 
18x $78 



$ 1422.00 
$ 1404.00 
$ 690.00 
$ 000.00 
$ 3516.00 
$ 40.00 
$ 3556.00 



Checks for the statutory filing fee of $ 3516,00 and Assignment recordation fee of S 40.00 are attached, 
.^ou are also directed and authorized to charge or credit any difference or overpayment to Deposit Account 
jsfo. 1 9-4880. The Commissioner is hereby authorized to charge any fees under 37 C.F.R. 1.16 and 1.17 and 
ahy petitions for extension of time under 37 C.F.R. 1 .1 36 which may be required during the entire pendency of 
the application to Deposit Account No. 19-4880. A duplicate copy of this transmittal letter is attached. 

Priority is claimed from: 



Japanese Patent Application 



Filing Date 



7293/1999 
15904/1999 
17250/1999 
26323/1999 



January 14, 1999 
January 25, 1999 
January 26, 1999 
February 3, 1999 



Respectfully submitted, 
SUGHRUE, MION, ZINN, MACPEAK & SEAS 
Attorneys for Applicant(s) 




DM:tnj 



Registration No. 23,063 



METHOD OF AND SYSTEM FOR IMAGE PROCESSING AND RECORDING MEDIUM 

FOR CARRYING OUT THE METHOD 

BACKGROUND OF THE INVENTION 
Field of the Invention 

This invention relates to a method of and system for 
carrying out an image processing such as a processing for 
enhancing a predetermined frequency component of an image signal . 
This invention further relates to a computer-readable recording 
medium loaded with program for causing a computer to perform 
the image processing in accordance with the method. 
Description of the Related Art 

We have proposed various image processing methods and 
systems for improving diagnostic performance of a radiation 
image by carrying out on a radiation image signal representing 
the radiation image, for instance, a frequency enhancement 
processing or a dynamic range compression processing by use of 
an unsharp mask image signal (will be referred to as ''unsharp 
image signal'', hereinbelow) . See, for instance, Japanese 
Unexamined Patent Publication Nos. 55(1980)-163472, 
55 (1980) -87953, 3 ( 1991 ) -222577 , 10 ( 1998 ) -753 95 , and 
10(1998)-171983. For example, in the frequency enhancement 
processing, a predetermined spatial frequency component of an 
original image signal is enhanced by subtracting an unsharp image 
signal S^^ from the original image signal S^^igf and adding the 
remainder multiplied by a coefficient of enhancement p to the 



original image signal S^^^g. This is represented by the following 
formula ( 1 ) . 

Sproc=Sorg+Px(S,,g-S,J (1) 

wherein Sp^^^ is a frequency-enhanced image signal, S^^g is an 
5 original image signal, S^^ is an unsharp image signal and p is 

a coefficient of enhancement. 

Further, in Japanese Unexamined Patent Publication No. 
10 ( 1998 )-75395 , there is disclosed a method of preventing 
generation of an artifact in the frequency-enhanced image signal 
by adjusting the frequency response characteristic of the add 
% signal to be added to the original image signal. In this method, 

fit a plurality of unsharp image signals, which are different from 

\j% each other in frequency response characteristic, that is, in 

1;^' sharpness, are prepared, differences between two of the original 

|5 image signal and the unsharp image signals are taken, thereby 

'J- making a plurality of band-limited signals respectively 

S representing frequency components in limited frequency bands 

of the original image signal, the band-limited signals thus 
obtained are transformed into signals of desired values by use 
2 0 of different transformation functions, and the add signal is 

made by adding up the suppressed band-limited signals. This is 
represented, for instance, by the following formulae (2). 

25 =fi(S,,g-S,3l)+f2(S,3l-S,,2)+ - 

+ fk(S,3k-l-S,3k)+ - +f^(S,,N-l-S,,N) (2) 
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wherein S^^^^ is a processed image signal, S^^g is an original image 
signal, S^3k (k:=l to N) is an unsharp image signal, fk(k=l to N) 
is a transformation function, and p(S^^g) is a coefficient of 
enhancement determined on the basis of the original image signal . 

Further, in Japanese Unexamined Patent Publication No. 
10 ( 1998) -17 1983, there is disclosed a method of preventing 
generation of an artifact in the processed signal when both the 
frequency enhancement processing and the dynamic range 
compression processing are to be carried out. In this method, 
a plurality of band-limited signals are made in the manner 
described above, a high frequency component signal representing 
high frequency components of the original image signal and a 
low frequency component signal representing low frequency 
components of the original image signal are obtained on the basis 
of the band-limited signals, and the frequency enhancement 
processing and the dynamic range compression processing are 
carried out by adding the high frequency component signal and 
the low frequency component signal to the original image signal. 
This is represented, for instance, by the following formulae 
(3). 

Sproc=Sorg+P(So,g) X F,,,(S,,g, S,^!, S,32 , S,,N) 

= { f nl ( Sorg-Susl ) +f n2 ( ^u.-i^'^.s^ ) "* 

+ f,,(S,,k-l-S,3k)-h - +fuN(Su3N-l-S,3N) 



= { f ( S,,,-S,3l ) +f d2 ( Susl"S,,2 ) + - 

+ fak(S.,k-l-S,3k)+ - +f,^(S,,N-l-S,3N) 

(3) 

wherein S^^^^ is a processed image signal, S^^^ is an original image 
signal, S^3k (k=l to N) are unsharp image signals, f^k(k=l to N) 
are transformation functions for obtaining the high frequency 
component signal, fak(k=l to N) is a transformation function for 
obtaining the low frequency component signal, p(S^^g) is a 
coefficient of enhancement determined on the basis of the 

original image signal, and D{So^g-Fdj.c (S^^g/ S^gl, S^^2 r S^^N)} is 

a coefficient of dynamic range compression determined on the 
basis of the low frequency component signal, D being a function 

for transforming D{S^^g-Fdrc(Sorg/ S^^l, S^^2 ^ S^^N)}. 

In the frequency enhancement processing and the dynamic 
range compression processing (will be representatively referred 
to "as the transformation processing", hereinbelow) , the 
frequency response characteristic of the add signal to be added 
to the original image signal can be adjusted by changing the 
definition of the transformation functions and the like for 
transforming the band-limited signals. Accordingly, a 
processed image signal having a desired frequency response 
characteristic, e.g., suitable for preventing generation of an 
artifact, can be obtained by properly defining the 
transformation functions. However, it is not easy to know how 
to define the transformation functions in order to obtain a 
desired result. Therefore, there has been proposed, in Japanese 



unexamined Patent Publication No. 10 ( 1998 )-63838 , a method in 
which a processed image signal having a desired frequency 
response characteristic is easily obtained by designating a 
desired frequency response characteristic for a processed image 
signal and determining parameters for defining the 
transformation functions (this parameter will be referred to 
as ''transformation function defining parameter", hereinbelow) 
on the basis of the designated frequency response 
characteristic . 

The unsharp image signals used in the aforesaid 
transformation processing are made by thinning picture elements 
by filtering picture elements of the original image signal at 
predetermined intervals and interpolating like number of 
picture elements • As the filtering processing, a processing for 
removing high frequency components from the original image 
signal by use of a low-pass filter, more specifically a 
processing for calculating an average or a weighted average of 
the values of picture elements in the filter, has been carried 
out. In the filtering processing carried out in order to obtain 
a plurality of unsharp image signals in Japanese Unexamined 
Patent Publication No. 10 ( 1998 ) -75395 or the like, the unsharp 
image signals are obtained by filtering the original image signal 
to obtain an image signal with less picture elements, further 
filtering the image signal with less picture elements , and 
interpolating picture elements into the image signal with less 
picture elements obtained at each filtering stage so that the 



number of the picture elements in the image signal becomes equal 
to that in the original image signal. 

Each unsharp image signal is thus made on the basis of 
the original image signal, which is obtained by reading an 
original image at a predetermined read density by use of an image 
reader and digitizing the image signal thus obtained into a 
digital image signal which can reproduce an image at a 
predetermined picture element density- It has been known that 
frequency components lower than a certain frequency determined 
according to the picture element density (a Nyquist rate) can 
be correctly reproduced when a digitized image signal is to be 
reproduced as a printed output. That is, since being determined 
taking into account the level of image quality required upon 
reproduction, the read density, i.e., the picture element 
density is not constant. 

For example, in a radiation image read-out and reproducing 
system, where a radiation image of a human body recorded on a 
stimulable phosphor sheet is read out as a digital image by 
scanning the stimulable phosphor sheet with a laser beam, the 
read density or the picture element density varies depending 
on the size of the stimulable phosphor sheet and can be freely 
set by an user. 

When image signals different in picture element density 
or Nyquist frequency are subjected to the same filtering 
processing using the same low-pass filters and then to the same 
interpolation, the frequency characteristics of the obtained 



band-limited signals (more specifically the frequency bands of 
the obtained band-limited signals) differ depending on the 
picture element density. This means that, for instance, when 
a pair of image signals having different picture element 
densities are obtained by reading an original image at different 
read densities, and a frequency enhancement processing or a 
dynamic range compression processing is carried out on the image 
signals by use of band-limited signals obtained on the basis 
of the same unsharp image signals, the enhanced frequency band 
or the compressed frequency band differ between the two original 
image signals. 

In order to overcome this problem, there has been proposed, 
in Japanese Unexamined Patent Publication No. 10 ( 1998 ) -63 83 7 , 
a method in which unsharp image signal are obtained by selecting 
coefficients of filter from a list of coefficients of filter 
according to information on the picture element density of the 
original image signal and filtering the original image signal 
by use of filters of the selected coefficients of filter. That 
is, when original image signals, for instance, respectively read 
at read densities of 5 lines/mm and 6.7 limes /mm are filtered 
by use of the same low-pass filter, the obtained two band-limited 
signals differ from each other in frequency band. However in 
the method proposed by the above identified Japanese Unexamined 
Patent Publication, the two original image signals are filtered 
by different low-pass filters and accordingly the obtained two 
band-limited signals can be substantially the same in frequency 



band. Accordingly, unsharp image signals of the same frequency 
bands can be obtained irrespective of the picture element density 
of the original image signals, whereby band-limited signals of 
the same frequency characteristics can be made and a desired 
transformation processing, e.g., the aforesaid frequency 
enhancement processing, can be constantly carried out in the 
same manner. 

However, since energy of a band-limited image represented 
by a band-limited signal, that is, the peak value of response 
of a band-limited signal, varies depending on the picture element 
density, the method disclosed in Japanese Unexamined Patent 
Publication No. 10 (1998) -63 83 7 cannot make constant the 
response of the band-limited signals in the same frequency band 
irrespective of the picture element density though can make the 
frequency bands of the band-limited signals equal to each other. 
Accordingly, even if a processing is carried out to enhance a 
band-limited signal in a particular frequency band, the response 
characteristic of the enhanced band-limited signal delicately 
varies depending on the picture element density. 

In the image processing system described above, there is 
sometimes input an object original image signal (an original 
image signal to be processed) representing an image which differs 
in resolution from images which are normally processed by the 
image processing system. (The resolution of the images which 
are normally processed by the image processing system will be 
referred to as "the reference resolution" , hereinbelow) In such 
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a case, when the transformation processing such as the frequency 
enhancement processing is carried out on the object original 
image signal by use of the transformation functions which have 
been determined for original image signals representing images 
at the reference resolution, there is a fear that the frequency 
response characteristic of the image reproduced on the basis 
of the processed image signal obtained from the object original 
image signal becomes different from that of the image reproduced 
on the basis of the original image signal representing an image 
at the reference resolution. This problem may be overcome, for 
instance, by preparing a plurality of groups of transformation 
functions and using one group of the transformation functions 
according to the resolution of the image represented by the 
object original image signal. However, this approach is 
disadvantageous in that the number of the transformation 
functions becomes too large and management of the transformation 
functions becomes too troublesome. 

Further, though there have been known, as formats for 
compressing an original image signal, various formats such as 
JPEG, GIF, TIFF and the like, there is recently proposed a format 
in which an original image signal is hierarchically decomposed 
by resolution into hierarchical data, and the hierarchical data 
in each hierarchy is encoded and compressed. in this 
compression format, specifically, an original image signal is 
decompos^ed by wavelet transformation or the like into a plurality 
of hierarchical image signals, each having a resolution of 1/2^ 



times that of the original image signal, and the hierarchical 
image signals are encoded in the hierarchical sequence and 
compressed into a single file. 

The compression format has the following features. 

( 1 ) Since the image signal is not processed block by block unlike 
the DOT format employed in the conventional JPEG, artifact like 
block distortion is not generated. 

(2) Since the image signals are hierarchically encoded, 
information on necessary resolutions has only to be transferred 
upon transfer of the image signals, which results in efficient 
image transfer. 

(3) Since the image signal is decomposed into multiple 
resolutions, various image processing such as frequency 
enhancement processing can be relatively easily carried out. 

(4) Since space decomposition and frequency decomposition can 
be simultaneously carried out by multiple resolution analysis, 
an orthogonal transformation can be carried out over a wide range 
on a low frequency region, which largely affects the encoding 
efficiency, whereas over a narrow region on a high frequency 
region. Accordingly, even if quantization noise is generated 
around an edge of the image, spatial spread of the noise can 
be suppressed so that the noise becomes less apt to be recognized. 

Further, there have been proposed various file formats 
such as a FlashPix file proposed by Eastman Kodak in which data 
of different kinds can be stored in a single file. The aforesaid 
hierarchical image signals can be also stored in such a FlashPix 
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standard file- 

By decomposing an original image signal into multiple 
resolutions, it is possible to construct the original image 
signal by a plurality of hierarchical image signals, each having 
5 a resolution of 112^ times that of the original image signal. 
This makes it feasible to reconstruct an image on the basis of 
a part of the hierarchical image signals which is selected 
according to the image quality required by the output system. 
That is, when a high quality image is to be reproduced as in 

10 a printer, an image signal which can reproduce a high quality 

image equivalent to the original image in resolution can be 
obtained by reconstructing the image signal on the basis of 
hierarchical image signals up to that of the highest resolution. 
To the contrast, in the case of, for instance, a CRT which cannot 

15 reproduce an image in a quality so high as a printer, an image 

signal which can reproduce an image suitable for the CRT though 
not so high as the original image in resolution can be obtained 
by reconstructing the image signal on the basis of hierarchical 
image signals not including the highest resolution hierarchical 

2 0 image signal and enlarging or contracting the image signal, if 

necessary. 

However since the hierarchical image signals each 
representing an image lower in resolution (such hierarchical 
image signals will be referred to as ''lower hierarchical image 
25 signals" hereinbelow) than that represented by the original 

image signal differ from the original image signal in frequency 
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response characteristic, if the transformation functions for 
the original image signal are employed as they are in frequency 
enhancement processing of the lower hierarchical image signals, 
there is a fear that an image signal which is different in 
5 frequency response characteristic from an image signal obtained 

by carrying out the frequency enhancement processing on the 
original image signal can be obtained. This may be overcome by 
preparing a number of transformation functions conforming to 
various resolutions and employing transformation functions 
AO according to the resolution of the image signal to be processed. 

J: However this approach is disadvantageous in that the number of 

nj transformation functions to be managed becomes too large and 

iji management of transformation functions becomes too troublesome. 

^ This problem occurs not only when carrying out frequency 

enhancement processing on image signals obtained by decomposing 
an original image signal into multiple resolutions but also when 
p carrying out frequency enhancement processing on an image signal 

in order to reproduce an image lower than an original image 
signal . 

20 Further, for example, a radiation image reproduced at 

pitches of 10 lines/mm by doubling and interpolating an original 
image signal read out from a stimulable phosphor sheet at a read 
density of 5 lines/mm is inferior in sharpness to a radiation 
image reproduced at pitches of 10 lines/mm on the basis of an 

25 original image signal read out from a stimulable phosphor sheet 

at a read density of 10 lines/mm though the sizes of images are 
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the same. This is because high frequency components of the 
original image signal are weakened depending on the 
characteristic of the filter employed for changing the picture 
element density of the original image signal^ that is, for 
5 contracting and interpolating the original image signal and the 

characteristic of the filter for doubling the contracted 
original image signal and because the former image is different 
from the latter image in frequency response characteristic. 
Further 

i*0 when a low resolution image is reproduced on the basis of a 

4^ hierarchical image signal obtained by decomposing an original 

nj image signal representing an original image into multiple 

kp resolutions , the obtained image becomes inferior to the original 

rH image in sharpness depending on the wavelet transformation 

functions employed in wavelet transformation. Further when a 
4: low resolution image is enlarged to the same size as the original 

P image, high frequency components of the original image signal 

are weakened depending on the characteristic of the filter for 
enlargement and interpolation and the obtained image becomes 
20 inferior to the original image in sharpness. The same problem 

occurs also when the original image is to be enlarged or 
contracted to a desired size. Accordingly when the aforesaid 
frequency enhancement processing is carried out on an image 
signal representing an image less sharp than the original image, 
25 the obtained image becomes different in impression from an image 
obtained by carrying out the same frequency enhancement 



13 



processing on the original image signal. 

SUMMARY OF THE INVENTION 
In view of the foregoing observations and description, 
the primary object of the present invention is to provide a method 
of and a system for image processing which can constantly carry 
out a desired transformation processing, e.g., a processing for 
enhancing a particular frequency component, so as to result in 
the same result irrespective of the picture element density of 
the original image signal. 

Another object of the present invention is to provide a 
method of and a system for image processing which can constantly 
carry out frequency enhancement processing so as to result in 
the same result irrespective of the resolution of the image to 
be processed and without necessity of storing a large number 
of transformation functions. 

Another object of the present invention is to provide a 
method of and a system for image processing which can constantly 
carry out frequency enhancement processing so as to result in 
the same result irrespective of the resolution of the image to 
be reproduced and without necessity of storing a large number 
of transformation functions. 

Another object of the present invention is to provide a 
method of and system for image processing which can carry out 
frequency enhancement processing on an image which is different 
from an original image in picture element density and/or scale 
so that the frequency response characteristic of the obtained 



image becomes equivalent to that of the original image. 

Still another object of the present invention is to provide 
a computer-readable recording medium loaded with program for 
causing a computer to perform the image processing in accordance 
with the method of the present invention. 

In accordance with a first aspect of the present invention, 
there is provided an image processing method for obtaining a 
processed image signal from an original image signal 
representing an original image having a certain picture element 
density in which a plurality of intermediate image signals which 
are different in frequency band are made on the basis of the 
original image signal, a plurality of transformed image signals 
are obtained by carrying out a transformation processing on the 
respective intermediate image signals on the basis of respective 
transformation functions, and a processed image signal is 
obtained from the transformed image signals, wherein the 
improvement comprises the step of 

defining said transformation functions by determining 
transformation function defining parameters for the 
transformation functions on the basis of the picture element 
density of the original image. 

The "transformation processing" includes, for instance, 
a frequency enhancement processing represented by formulae (2) 
for enhancing a particular frequency component and a dynamic 
range compression processing represented by formulae (3) for 
reducing the contrast of the high density range and/or the low 



density range, thereby narrowing the dynamic range which is the 
difference between the maximum density and the minimum density 
of the original image. 

The ^'picture element density of the original image'' may 
be either input by the operator or may be automatically obtained 
when the original image signal is processed. In the latter case, 
information on the picture element density is attached to the 
original image signal in advance. The ''picture element density" 
may be, for instance, a value of resolution (e.g., in dpi) which 
represents the relation between the size of the original image 
and the sampling intervals for obtaining the original image 
signal as well as the read density at which a radiation image 
recorded on a stimulable phosphor sheet is read. 

Further, the expression ''transformation function 
defining parameters for the transformation functions are 
determined on the basis of the picture element density of the 
original image signal" means, for instance, to obtain the 
parameters by resolving relational expressions of the aforesaid 
plurality of intermediate image signals, a desired frequency 
response characteristic and the transformation function 
defining parameters as simultaneous equations taking the 
plurality of intermediate image signals and the desired 
frequency response characteristic as known values and the 
transformation function defining parameters as variables, or 
to determine the transformation function defining parameters 
by gradually changing the parameters while watching the image 



reproduced on the basis of the processed image signal. As the 
transformation functions, various functions such as linear 
functions, non-linear functions, constants and the like can be 
used. 

It is preferred that said plurality of intermediate image 
signals which are different in frequency band be made by carrying 
out on the original image signal a filtering processing by use 
of filters whose coefficients of filter are determined on the 
basis of the picture element density of the original image, 
thereby making a plurality of unsharp image signals which are 
different in frequency response characteristic, making a 
plurality of band-limited signals representing the signals in 
the respective frequency bands of the original image signal on 
the basis of the unsharp image signals and the original image 
signal, and taking the band-limited signals as the intermediates 
image signals. 

The "unsharp image signals" are image signals which are 
equivalent to the original image signal in the number of picture 
elements but represent images which are lower in sharpness than 
that represented by the original image signal. The unsharp 
image signals are made by thinning picture elements by filtering 
picture elements of the original image signal at predetermined 
intervals, further thinning picture elements by filtering 
picture elements of the thinned image signal, repeating these 
procedures, and interpolating picture elements into the image 
signal with less picture elements obtained at each filtering 



stage so that the number of the picture elements in the image 
signal becomes equal to that in the original image signal. The 
filtering processing is carried out by use of filters having 
coefficients of filter determined on the basis of the picture 
element densities, and is described in detail in the above 
identified Japanese patent publication, Japanese Unexamined 
Patent Publication No* 10 ( 1998 ) -63 837 . 

The "plurality of band-limited signals representing the 
signals in the respective frequency bands of the original image 
signal" may be made, for instance, by taking differences between 
unsharp image signals in adjacent frequency bands, differences 
between the original image signal and the respective unsharp 
image signals, or differences between any other combinations 
of the unsharp image signals. 

In accordance with a second aspect of the present invention, 
there is provided an image processing system for obtaining a 
processed image signal from an original image signal 
representing an original image having a certain picture element 
density comprising an intermediate image signal making means 
which makes a plurality of intermediate image signals which are 
different in frequency band on the basis of the original image 
signal, and a transformation means which obtains a plurality 
of transformed image signals by carrying out a transformation 
processing on the respective intermediate image signals on the 
basis of respective transformation functions, and obtains a 
processed image signal from the transformed image signals, 
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wherein the improvement comprises that 

a transformation function defining means which determines 
transformation function defining parameters for the 
transformation functions on the basis of the picture element 
density of the original image and defines the transformation 
functions . 

It is preferred that the intermediate image signal making 
means comprises an unsharp image signal making means which makes 
a plurality of unsharp image signals which are different in 
frequency response characteristic by carrying out on the 
original image signal a filtering processing by use of filters 
whose coefficients of filter are determined on the basis of the 
picture element density of the original image, and a band-limited 
signal making means which makes, as the intermediates image 
signals, a plurality of band-limited signals representing the 
signals in the respective frequency bands of the original image 
signal on the basis of the unsharp image signals and the original 
image signal. 

In accordance with a third aspect of the present invention, 
there is provided a computer-readable recording medium loaded 
with program for causing a computer to perform the image 
processing method in accordance with the first aspect of the 
present invention. 

In accordance with the method and the system of the first 
and second aspects of the present invention, since the 
transformation function defining parameters are obtained on the 



basis of the picture element density of the original image and 
the transformation processing is carried out according to the 
transformation functions defined by the transformation function 
defining parameters, the parameters can be set so that the 
intermediate image signals different in frequency band are 
transformed into transformed image signals which are the same 
in frequency response characteristic irrespective of the 
picture element density of the original image, whereby the 
transformation processing can be carried out so that a processed 
image signal can reproduce an image having a constant frequency 
response characteristic without affected by the picture element 
density of the original image. 

In accordance with a fourth aspect of the present invention, 
there is provided an image processing method for obtaining a 
processed image signal from an original image signal 
representing an original image at a certain resolution, in which 
a plurality of intermediate image signals which are different 
in frequency band are made on the basis of the original image 
signal, a plurality of transformed image signals are obtained 
by carrying out a transformation processing on the respective 
intermediate image signals on the basis of respective 
transformation functions, and a processed image signal is 
obtained from the transformed image signals, wherein the 
improvement comprises that 

the transformation functions for processing the 
intermediate image signals obtained from an object original 
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image signal to be processed are determined by correcting, 
according to the resolution of the object original image signal, 
reference transformation functions for transforming the 
intermediate image signals obtained from a reference original 
image signal representing an image having a reference resolution, 
the reference transformation functions having been determined 
and stored in advance. 

The "transformation processing" includes, for instance, 
a frequency enhancement processing represented by formulae (2) 
for enhancing a particular frequency component and a dynamic 
range compression processing represented by formula (3) for 
reducing the contrast of the high density range and/or the low 
density range, thereby narrowing the dynamic range which is the 
difference between the maximum density and the minimum density 
of the original image. 

The "resolution of the object original image signal" may 
be either input by the operator or automatically obtained when 
the original image signal is processed. In the latter case, 
information on the resolution is attached to the original image 
signal in advance. The "resolution" may be, for instance, a 
value of resolution (e.g. , in dpi) which represents the relation 
between the size of the original image and the sampling intervals 
for obtaining the original image signal as well as the read 
density at which a radiation image recorded on a stimulable 
phosphor sheet is read. 

It is preferred that said plurality of intermediate image 



signals which are different in frequency band be made by carrying 
out on the original image signal a filtering processing by use 
of filters having predetermined coefficients of filter, thereby 
making a plurality of unsharp image signals which are different 
in frequency response characteristic, making a plurality of 
band-limited signals representing the signals in the respective 
frequency bands of the original image signal on the basis of 
the unsharp image signals and the original image signal, and 
taking the band-limited signals as the intermediates image 
signals . 

The "unsharp image signals" are image signals which are 
equivalent to the original image signal in the number of picture 
elements but represent images which are lower in sharpness than 
that represented by the original image signal. The unsharp 
image signals are made by thinning picture elements by filtering 
picture elements of the original image signal at predetermined 
intervals, further thinning picture elements by filtering 
picture elements of the thinned image signal, repeating these 
procedures, and interpolating picture elements into the image 
signal with less picture elements obtained at each filtering 
stage so that the number of the picture elements in the image 
signal becomes equal to that in the original image signal. The 
filtering processing is carried out by use of filters having 
coefficients of filter determined on the basis of the picture 
element densities, and is described in detail in the above 
identified Japanese patent publication, Japanese Unexamined 



Patent Publication No. 10 ( 1998 ) -63837 . 

The ''plurality of band-limited signals representing the 
signals in the respective frequency bands of the original image 
signal" may be made, for instance, by taking differences between 
unsharp image signals in adjacent frequency bands, differences 
between the original image signal and the respective unsharp 
image signals, or differences between any other combinations 
of the unsharp image signals. 

In accordance with a fifth aspect of the present invention, 
there is provided an image processing system comprising an 
intermediate image signal making means which makes a plurality 
of intermediate image signals which are different in frequency 
band on the basis of an original image signal representing an 
original image having a predetermined resolution, and a 
transformation means which obtains a plurality of transformed 
image signals by carrying out a transformation processing on 
the respective intermediate image signals on the basis of 
respective transformation functions, and obtains a processed 
image signal from the transformed image signals, wherein the 
improvement comprises that 

a transformation function calculating means which 
calculates the transformation functions for processing the 
intermediate image signals obtained from an object original 
image signal to be processed by correcting, according to the 
resolution of the object original image signal, reference 
transformation functions for transforming the intermediate 



image signals obtained from a reference original image signal 
having a reference resolution, the reference transformation 
functions having been determined and stored in advance* 

It is preferred that the intermediate image signal making 
5 means comprises an unsharp image signal making means which makes 
a plurality of unsharp image signals which are different in 
frequency response characteristic by carrying out on the 
original image signal a filtering processing by use of filters 
whose coefficients of filter are determined on the basis of the 
10 U resolution of the original image signal, and a band-limited 
^ signal making means which makes, as the intermediates image 
signals, a plurality of band-limited signals representing the 
signals in the respective frequency bands of the original image 
^ signal on the basis of the unsharp image signals and the original 
15 image signal. 

4;^ In accordance with a sixth aspect of the present invention, 

r.J there is provided a computer-readable recording medium loaded 
with program for causing a computer to perform the image 
processing method in accordance with the fourth aspect of the 
2 0 present invention. 

In accordance with the method and the system of the fourth 
and fifth aspects of the present invention, since the 
transformation functions for transforming the intermediate 
image signals obtained from an object original image signal 
25 representing an image whose resolution is different from the 

reference resolution are obtained by correcting the reference 
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transformation functions according to the resolution of the 
object original image signal, the transformation functions can 
be set so that the intermediate image signals are transformed 
into transformed image signals which are the same in frequency 
response characteristics irrespective of the picture element 
density of the original image signal, whereby the transformation 
processing can be carried out so that a processed image signal 
can reproduce an image having a constant frequency response 
characteristic without affected by the resolution of the 
original image signal • 

Further, since the transformation functions for original 
image signals having resolutions different from the reference 
resolution are obtained by correcting the reference 
transformation functions, it is not necessary to prepare a 
plurality of transformation functions for various resolutions, 
whereby trouble to manage a plurality of transformation 
functions can be saved. 

In accordance with a seventh aspect of the present 
invention, there is provided an image processing method for 
obtaining a processed image signal from an original image signal 
representing an original image having a certain resolution, in 
which a plurality of intermediate image signals which are 
different in frequency band are made on the basis of the original 
image signal, a plurality of transformed image signals are 
obtained by carrying out a transformation processing on the 
respective intermediate image signals on the basis of 



transformation functions which are set according to the 
frequency bands of the respective intermediate image signals 
and a processed image signal is obtained from the transformed 
image signals, wherein the improvement comprises the steps of 

preparing reference transformation functions for 
transforming the intermediate image signals obtained from a 
reference original image signal representing a reference 
original image having a reference resolution which are set 
according to the frequency bands of the respective intermediate 
image signals, and, when an object original image signal to be 
processed represents an original image having a resolution lower 
than the reference resolution, setting the transformation 
functions for transforming the intermediate image signals 
obtained from the object original image signal to be equal to 
the respective reference transformation functions for the 
intermediate image signals in the frequency bands not higher 
than the frequency band corresponding to the resolution of the 
original image represented by the object original image signal. 

The "transformation functions which are set according to 
the frequency bands of the respective intermediate image 
signals" are six in number when the intermediate image signals 
are in six frequency bands. The transformation functions may 
be constants. The expression "the respective reference 
transformation functions for the intermediate image signals in 
the frequency bands not higher than the frequency band 
corresponding to the resolution of the original image 



represented by the object original image signal" means, for 
instance, when the resolution of the original image represented 
by the object original image signal is 1/2^ times that of the 
reference resolution, the reference transformation functions 
for the k-th intermediate image signal as numbered from the high 
frequency side and those in the frequency bands lower than the 
k-th intermediate image signals. For example, when six 
intermediate image signals are obtained from the reference 
original image signal and six transformation functions (f^ to 
fg, fi being in the highest frequency band) are set for the 
respective image intermediate signals with the resolution of 
the original image represented by the object original image 
signal being 1/2 of that of the reference original image, "the 
respective reference transformation functions for the 
intermediate image signals in the frequency bands not higher 
than the frequency band corresponding to the resolution of the 
original image represented by the object original image signal" 
are the transformation functions to fg. 

In the image processing method in accordance with the 
seventh aspect of the present invention, it is preferred that 
said plurality of intermediate image signals be band-limited 
signals which are made by making, on the basis of the original 
image signal (the reference original image signal or the object 
original image signal), a plurality of unsharp image signals 
which are different in frequency response characteristic, and 
making a plurality of band-limited signals representing the 



signals in the respective frequency bands of the original image 
signal on the basis of the unsharp image signals and the original 
image signal - 

The ''unsharp image signals" are image signals which are 
equivalent to the original image signal in the number of picture 
elements but represent images which are lower in sharpness than 
that represented by the original image signal. The unsharp 
image signals are made by thinning picture elements by filtering 
picture elements of the original image signal at predetermined 
intervals, further thinning picture elements by filtering 
picture elements of the thinned image signal, repeating these 
procedures to obtain a plurality of low resolution image signals 
respectively representing images whose resolutions are 1/2" of 
that of the image represented by the reference original image 
signal, and interpolating picture elements into the low 
resolution image signal with less picture elements obtained at 
each filtering stage so that the number of the picture elements 
in the image signal becomes equal to that in the original image 
signal . 

The "plurality of band-limited signals representing the 
signals in the respective frequency bands of the original image 
signal" may be made, for instance, by taking differences between 
unsharp image signals in adjacent frequency bands, differences 
between the original image signal and the respective unsharp 
image signals, or differences between any other combinations 
of the unsharp image signals. 
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Further^ it is preferred that information on the 
resolution of the original image represented by the object 
original image signal be obtained and the step of setting the 
transformation functions for transforming the intermediate 
image signals obtained from the object original image signal 
be executed on the basis of the information. 

The resolution information may be input as a value by the 
operator through a keyboard or by the operator selecting a value 
from a plurality of resolutions displayed on a control screen. 
Further, the resolution information may be attached to the object 
original image signal and may be detected by the image processing 
system. The resolution information may be input in any way so 
long as the system can recognize the resolution of the original 
image represented by the object original image signal. 

In the image processing method in accordance with the 
seventh aspect of the present invention, when the original image 
represented by the object original image signal to be processed 
has a resolution lower than the reference resolution, the 
transformation functions for transforming the intermediate 
image signals obtained from the object original image signal 
are set to be equal to the respective reference transformation 
functions for the intermediate image signals in the frequency 
bands not higher than the frequency band corresponding to the 
resolution of the original image represented by the object 
original image signal, when the resolution of the original 
image represented by the object original image signal is 



of the resolution of the reference original image represented 
by the reference original image signal, the peak frequencies 
of the intermediate image signals obtained from the object 
original image signal are the same as those of a part of the 
reference intermediate image signals (obtained from the 
reference original image signal) which are in frequency bands 
not higher than the frequency band corresponding to the 
resolution of the original image represented by the object 
original image s ignal . Accordingly , when the intermediate 
image signals obtained from the object original image signal 
are transformed on the basis of the reference transformation 
functions for the part of the reference intermediate image 
signals which are in frequency bands not higher than the 
frequency band corresponding to the resolution of the original 
image represented by the object original image signal, the 
obtained transformed image signals are substantially the same 
as the reference transformed image signals in frequency response 
characteristics • The processed image signal can reproduce an 
image which is substantially the same in frequency response 
characteristic as an image reproduced on the basis of the 
processed image signal obtained from the reference original 
image signal- Thus in accordance with the image processing 
method of the seventh aspect of the present invention, a 
processed image signal which can reproduce an image which is 
substantially constant in frequency response characteristic can 
be obtained irrespective of the resolution of the object original 
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image signal* 

Further, since the transformation functions for obtaining 
the processed image signal from the reference original image 
signal are employed to obtain the processed image signal from 
the object original image signal, it is not necessary to prepare 
a plurality of transformation functions for various resolutions , 
whereby the system can be simplified in structure and trouble 
to manage a plurality of transformation functions can be saved. 

In accordance with an eighth aspect of the present 
invention, there is provided an image processing system for 
carrying out the image processing method in accordance with the 
seventh aspect of the present invention. 

In accordance with a ninth aspect of the present invention, 
there is provided a computer-readable recording medium loaded 
with program for causing a computer to perform the image 
processing method in accordance with the seventh aspect of the 
present invention. 

In accordance with a tenth aspect of the present invention, 
there is provided an image processing method for making a 
contraction-processed image signal representing an image having 
a desired resolution lower than that of an image represented 
by a reference processed image signal which is obtained by making 
first to n-th reference low resolution image signals from a 
reference original image signal representing a reference 
original image, the reference low resolution image signals being 
different from each other in frequency band and respectively 



representing images whose resolutions are 1/2^ (k=l to n) of the 
resolution of the reference original image signal; making 
reference low-resolution band-limited signals on the basis of 
the reference low resolution image signals; obtaining a 
5 plurality of reference transformed image signals by carrying 
out a predetermined transformation processing on the reference 
low resolution band-limited signals on the basis of a plurality 
of reference transformation functions which are set according 
to the respective frequency bands of the reference low resolution 
10 O band-limited signals; and carrying out a predetermined 
4- processing on the reference transformed image signals, the image 
OJ processing method characterized by the steps of 

taking as an original image signal one of said first to 
n-th reference low resolution image signals representing an 
15 image whose resolution is the closest to said desired resolution 

4^ in the images represented by the reference low resolution image 
O signals; making first to m-th low resolution image signals from 
the original image signal, the low resolution image signals being 
different from each other in frequency band and respectively 
20 representing images whose resolutions are 1/2^^' (k'=l to m) of 

the resolution of the original image signal; making a low- 
resolution band-limited signals on the basis of the low 
resolution image signals; obtaining a plurality of transformed 
image signals by carrying out said predetermined transformation 
25 processing on the low resolution band-limited signals on the 

basis of a part of said reference transformation functions which 
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are set for the respective frequency bands of the reference low 
resolution band-limited signals not higher than the frequency- 
band corresponding to the resolution closest to said desired 
resolution; carrying out said predetermined processing on the 
transformed image signals, thereby obtaining an intermediate 
processed image signal; enlarging or contracting the 
intermediate processed image so that the resolution of the image 
represented by the intermediate processed image signal becomes 
equal to the desired resolution; and taking the enlarged or 
contracted intermediate processed image signal as the 
contraction-processed image signal. 

In the image processing method in accordance with the tenth 
aspect of the present invention, out of the reference low 
resolution image signals obtained from the reference original 
image signal, a reference low resolution image signal 
representing an image whose resolution is closest to the desired 
resolution is taken as an object original image signal for 
obtaining a contraction-processed image signal, and a plurality 
of low resolution image signals are obtained from the object 
original image signal. Accordingly, the peak frequencies of the 
low resolution band-limited signals obtained from the object 
original image signal are the same as those of a part of the 
reference low resolution band-limited signals which are in 
frequency bands not higher than the frequency band corresponding 
to the resolution closest to the desired resolution. 
Accordingly, when the low resolution band-limited signals 



obtained from the object original image signal are transformed 
on the basis of the reference transformation functions for the 
part of the reference low resolution band-limited signals which 
are in frequency bands not higher than the frequency band 
corresponding to the resolution closest to the desired 
resolution, the obtained transformed image signals are 
substantially the same as the reference transformed image 
signals in frequency response characteristics. The 
contraction-processed image signal obtained by enlarging or 
contracting the intermediate processed image signal obtained 
from the transformed image signals so that the resolution of 
the image represented by the intermediate processed image signal 
becomes equal to the desired resolution can reproduce an image 
which is substantially the same in frequency response 
characteristic as an image reproduced on the basis of the 
reference processed image signal. Thus in accordance with the 
image processing method of the tenth aspect of the present 
invention, a contraction-processed image signal which can 
reproduce an image which is substantially constant in frequency 
response characteristic can be obtained irrespective of the 
resolution of the image to be reproduced. 

Further, since the transformation functions for obtaining 
the reference processed image signal are employed to obtain the 
contraction-processed image signal, it is not necessary to 
prepare a plurality of transformation functions for various 
resolutions, whereby the system can be simplified in structure 
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and trouble to manage a plurality of transformation functions 
can be saved. 

In accordance with an eleventh aspect of the present 
invention, there is provided an image processing system for 
carrying out the image processing method in accordance with the 
tenth aspect of the present invention. 

In accordance with a twelfth aspect of the present 
invention , there is provided a computer-readable recording 
medium loaded with program for causing a computer to perform 
the image processing method in accordance with the tenth aspect 
of the present invention. 

In accordance with a thirteenth aspect of the present 
invention, there is provided an image processing method for 
obtaining a processed image signal from an original image signal 
representing an original image having a certain picture element 
density, in which a plurality of intermediate image signals which 
are different in frequency band are made on the basis of the 
original image signal, a plurality of transformed image signals 
are obtained by carrying out a transformation processing on the 
respective intermediate image signals on the basis of respective 
transformation functions, and a processed image signal in which 
a predetermined frequency component is enhanced is obtained from 
the transformed image signals, wherein the improvement 
comprises that 

when an object original image signal to be processed is 
a contracted image signal obtained by carrying out a picture 



element density transformation processing on a reference 
original image signal representing an image having a reference 
picture element density, the transformation functions for 
processing the intermediate image signals obtained from the 
object original image signal are determined by correcting, 
according to the characteristic of the picture element density 
transformation processing, reference transformation functions 
for transforming the intermediate image signals obtained from 
the reference original image so that the frequency response 
characteristic of a desired frequency component of the processed 
image signal obtained from the object original image signal 
becomes equivalent to that of the processed image signal obtained 
from the reference original image signal, the reference 
transformation functions having been determined and stored in 
advance. 

It is preferred that when the processed image signal 
obtained from the contracted image signal is enlarged at a 
magnification, the transformation functions for processing the 
intermediate image signals obtained from the object original 
image signal be determined by correcting the reference 
transformation functions according to the characteristic of the 
picture element density transformation processing and the 
characteristic of the enlargement. 

In accordance with a fourteenth aspect of the present 
invention, there is provided an image processing method for 
obtaining a processed image signal from an original image signal 



representing an original image, in which a plurality of 
intermediate image signals which are different in frequency band 
are made on the basis of the original image signal, a plurality 
of transformed image signals are obtained by carrying out a 
transformation processing on the respective intermediate image 
signals on the basis of respective transformation functions, 
and a processed image signal in which a predetermined frequency 
component is enhanced is obtained from the transformed image 
signals, wherein the improvement comprises that 

when an enlarged processed image signal is to be obtained 
by enlarging the processed image signal at a desired 
magnification by carrying out an enlargement transformation on 
the processed image signal, the transformation functions for 
transforming the respective intermediate image signals are 
corrected according to the characteristic of the enlargement 
transformation so that the frequency response characteristic 
of a desired frequency component of the enlarged processed image 
signal becomes a predetermined frequency response 
characteristic . 

In the image processing methods in accordance with the 
thirteenth and fourteenth aspects of the present invention, it 
is preferred that said plurality of intermediate image signals 
be band-limited signals which are made by making, on the basis 
of the original image signal, a plurality of unsharp image 
signals which are different in frequency response 
characteristic, and making a plurality of band-limited signals 



representing the signals in the respective frequency bands of 
the original image signal on the basis of the unsharp image 
signals and the original image signal. 

In accordance with fifteenth and sixteenth aspects of the 
present invention, there are provided image processing systems 
for carrying out the image processing methods in accordance with 
the thirteenth and fourteenth aspects of the present invention. 

In accordance with seventeenth and eighteenth aspects of 
the present invention, there is provided a computer-readable 
recording medium loaded with program for causing a computer to 
perform the image processing methods in accordance with the 
thirteenth and fourteenth aspects of the present invention. 
BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a schematic block diagram showing the 
arrangement of an image processing system in accordance with 
a first embodiment of the present invention. 

Figure 2 is a schematic block diagram showing the 
arrangement of the unsharp image signal making means employed 
in the image processing system shown in Figure 1, 

Figures 3A and 3B are views showing examples of the filter 
to be used in the filtering processing. 

Figure 4 shows in detail the low-resolution image signal 
making processing. 

Figure 5 is a view showing an example of the filter employed 
in the interpolation/enlargement processing. 

Figure 6 is a view showing system for executing a frequency 
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enhancement processing as a non-linear processing^ 

Figure 7 shows the frequency response characteristics of 
band-limited signals obtained by carrying out, on an original 
image signal at a picture element density of 10 lines/mm, the 
filtering processing using the one-dimensional filter shown in 
Figure 3A, 

Figure 8 shows the frequency response characteristics of 
band-limited signals obtained by carrying out, on an original 
image signal at a picture element density of 6.7 lines /mm, the 
filtering processing using the one-dimensional filter shown in 
Figure 3B, 

Figure 9 shows the frequency response characteristics of 
band-limited signals obtained by carrying out, on an original 
image signal at a picture element density of 5 lines/mm, the 
filtering processing using the one-dimensional filter shown in 
Figure 3A, 

Figure 10 is a schematic block diagram showing the 
structure of the transformation function defining means. 

Figure 11 is a view showing the frequency response 
characteristic of the processed image signal obtained by 
processing the band-limited signals obtained from an original 
image signal at 10 lines/mm by the image processing system of 
the first embodiment. 

Figure 12 is a view showing the frequency response 
characteristic of the processed image signal obtained by 
processing the band-limited signals obtained from an original 



image signal at 6.7 lines/mm by the image processing system of 
the first embodiment. 

Figure 13 is a view showing the frequency response 
characteristic of the processed image signal obtained by 
processing the band-limited signals obtained from an original 
image signal at 5 lines /mm by the image processing system of 
the first embodiment. 

Figure 14 is a view showing an example of the user interface 
for designating the frequency response characteristics. 

Figure 15 is a flow chart for illustrating operation of 
the image processing system of the first embodiment. 

Figure 16 is a view showing system for executing a 
frequency enhancement processing and a dynamic range 
compression processing as a non-linear processing in an image 
processing system in accordance with a second embodiment of the 
present invention. 

Figure 17 is a view showing an example of the 
transformation functions used in the first transformation 
section in the second embodiment. 

Figure 18 is a view showing another example of the 
transformation functions used in the first transformation 
section in the second embodiment. 

Figure 19 is a view showing an example of the 
transformation functions used in the second transformation 
section in the second embodiment. 

Figure 20 is a view showing the frequency response 
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characteristic of the processed image signal Sp^^^ obtained by 
processing an original image signal at 10 lines/mm in the image 
processing system of the second embodiment. 

Figure 21 is a view showing the frequency response 
characteristic of the processed image signal S^^^^ obtained by 
processing an original image signal at 6.7 lines /mm in the image 
processing system of the second embodiment. 

Figure 22 is a view showing the frequency response 
characteristic of the processed image signal Sp^^^ obtained by 
processing an original image signal at 5 lines/mm in the image 
processing system of the second embodiment. 

Figure 23 is a schematic block diagram showing the 
arrangement of an image processing system in accordance with 
a third embodiment of the present invention. 

Figure 24 is a view showing the frequency response 
characteristics of the band-limited signals obtained from an 
original image signal S^^g having a picture element density of 
10 lines/mm. 

Figure 25 is a view showing the frequency response 
characteristic of the processed image signal S^^^^ obtained by 
carrying out the frequency enhancement processing by use of the 
transformation functions shown in table 3 on the original image 
signal S^^g at 10 lines/mm (solid line) and that of the processed 
image signal Sp^^^ obtained by carrying out the frequency 
enhancement processing by use of the transformation functions 
shown in table 3 on the original image signal S^^gat 300dpi (broken 



line) t 

Figure 2 6 is a view showing the frequency response 
characteristics of the band-limited signals obtained from an 
original image signal at 3 00dpi, 

Figure 27 is a view showing the frequency response 
characteristic of the processed image signal Sp^oc obtained by 
processing an original image signal at 300dpi by the image 
processing of the third embodiment of the present invention. 

Figure 28 is a flow chart for illustrating operation of 
the image processing system of the third embodiment. 

Figures 29A and 2 9B are views for illustrating different 
manners for selecting a reference picture element density. 

Figure 3 0 is a view showing the relation between the number 
of the band-limited signals and the resolution. 

Figure 31 is a view showing the frequency response 
characteristics of the band-limited signals obtained from an 
original image signal at 20 0dpi, 

Figure 32 shows the frequency response characteristics 
of the processed image signal S^^^^ obtained by processing the 
original image signal S^^g at 200dpi by use of the reference 
transformation functions (broken line) and the processed image 
signal Sp^^^^O obtained by processing the reference original image 
signal S^^^gO at the reference picture element density by use of 
the reference transformation functions (solid line). 

Figure 33 is a view showing the relation between the 
response XI [Fi] and the function XI [i]. 



Figure 3 4 is a view showing the relation between the 
functions XI [i]^ and Zl[L], 

Figure 35 is a view showing a system for executing a 
frequency enhancement processing and a dynamic range 
compression processing as a transformation processing in an 
image processing system in accordance with a fifth embodiment 
of the present invention. 

Figure 36 is a view showing an example of the 
transformation functions used in the first transformation 
section in the fifth embodiment. 

Figure 3 7 is a view showing another example of the 
transformation functions used in the first transformation 
section in the fifth embodiment. 

Figure 3 8 is a view showing an example of the 
transformation functions used in the second transformation 
section in the fifth embodiment. 

Figure 39A is a view showing the frequency response 
characteristic of the reference processed image signal obtained 
by transforming the band-limited signals from the reference 
image signal (having a picture element density of 5 lines /mm) 
by use of reference transformation functions which are set so 
that relatively low frequency components are enhanced and the 
frequency response characteristic of the processed image signal 
obtained by transforming the band-limited signals from the 
original image signal to be processed by use of the 
transformation functions which are obtained by correcting the 



reference transformation functions according to the picture 
element density (4 lines/ram) of the original image signal to 
be processed so that the frequency response characteristic of 
the processed image signal conforms to that of the reference 
processed image signal in frequency bands not lower than 1/2 
of the Nyquist frequency of the original image signal to be 
processed. 

Figure 3 9B is a view showing the frequency response 
characteristic of the reference processed image signal obtained 
by transforming the band-limited signals from the reference 
image signal (having a picture element density of 5 lines /mm) 
by use of reference transformation functions which are set so 
that relatively high frequency components are enhanced and the 
frequency response characteristic of the processed image signal 
obtained by transforming the band-limited signals from the 
original image signal to be processed by use of the 
transformation functions which are obtained by correcting the 
reference transformation functions according to the picture 
element density (4 lines /mm) of the original image signal to 
be processed so that the frequency response characteristic of 
the processed image signal conforms to that of the reference 
processed image signal in frequency bands not lower than 1/2 
of the Nyquist frequency of the original image signal to be 
processed. 

Figure 40 is a schematic block diagram showing the 
arrangement of an image processing system in accordance with 



a sixth embodiment of the present invention. 

Figures 41A to 41C are views for illustrating encoding 
the hierarchical image signals, 

Figure 42 is a schematic block diagram showing the 
arrangement of the unsharp image signal making means employed 
in the image processing system shown in Figure 40, 

Figure 43 is a view showing an example of the filter to 
be used in the filtering processing. 

Figure 44 shows in detail the low-resolution image signal 
making processing. 

Figure 45 is a view showing an example of the filter 
employed in the interpolation/enlargement processing. 

Figure 46 is a view showing system for executing a 
frequency enhancement processing. 

Figure 47 shows the frequency response characteristics 
of band-limited signals obtained from the reference original 
image signal. 

Figure 4 8 shows the frequency response characteristics 
of band-limited signals obtained from low frequency original 
image signal having a resolution of 1/2 of that of the reference 
original image signal. 

Figure 4 9 shows the frequency response characteristics 
of band-limited signals obtained from low frequency original 
image signal having a resolution of 1/4 of that of the reference 
original image signal. 

Figure 50 shows the frequency response characteristics 



of band-limited signals obtained from low frequency original 
image signal having a resolution of 1/8 of that of the reference 
original image signal, 

Figure 51 shows the frequency response characteristics 
of band-limited signals obtained from low frequency original 
image signal having a resolution of 1/16 of that of the reference 
original image signal. 

Figure 52 is a flow chart showing operation of the image 
processing system in accordance a sixth embodiment of the present 
invention. 

Figure 53 is a view showing correspondence between the 
unsharp image signals obtained from the reference original image 
signal S^^g and those obtained from an object low resolution image 
signal. 

Figure 54 is a view showing correspondence between the 
band-limited signals obtained from the reference original image 
signal S^^g and those obtained from an object low resolution image 
signal. 

Figure 55 is a view showing correspondence between the 
transformation functions for the band-limited signals obtained 
from the reference original image signal S^^g and those for the 
band-limited signals obtained from an object low resolution 
image signal. 

Figure 56 is a block diagram showing a unsharp image signal 
making means in the image processing system in accordance with 
a seventh embodiment of the present invention. 



Figure 57 is a view showing correspondence between the 
band-limited signals obtained from the reference original image 
signal S^^g and those obtained from object original image signals, 

Figure 5 8 is a view showing correspondence between the 
transformation functions for the band-limited signals obtained 
from the reference original image signal S^^g and those for the 
band-limited signals obtained from object original image 
signals , 

Figure 59 is a schematic block diagram showing the 
arrangement of an image processing system in accordance with 
a seventh embodiment of the present invention. 

Figure 60 is a view showing the frequency response 
characteristic of the picture element density transformation 
filter. 

Figure 61 is a view showing the frequency response 
characteristic of the interpolation filter. 

Figure 62 is a view showing the frequency response 
characteristic of the processed image signal Sp^,^^ obtained by 
carrying out the frequency enhancement processing by use of the 
transformation functions fj, set for the frequency enhancement 
processing of the base original image signal S^^g' on the 5-line 
transformed original image signal S,,g5 and by doubling the 
obtained intermediate processed image signal Sp^^^' by use of the 
interpolation filter whose filtering characteristic is F2 , 

Figure 63 is a view showing the frequency response 
characteristics of the band-limited signals obtained from the 



base original image signal S^^^' ^ 

Figure 64 is a flow chart showing operation of the image 
processing system in accordance the seventh embodiment of the 
present invention. 

Figure 65 is a schematic block diagram showing the 
arrangement of an image processing system in accordance with 
an eighth embodiment of the present invention. 

Figures 66A to 66C are views for illustrating encoding 
the hierarchical image signals. 

Figure 67 is a view showing frequency response 
characteristics of the low resolution image signals which are 
restored up to a resolution equal to 1/2 of that of the base 
original image signal S^^g' and up to a resolution equal to 1/4 
of that of the base original image signal S^^^' when the base 
original image signal S^^g' has been decomposed by use of wavelet 
transformation functions having first coefficients of filter. 

Figure 68 shows the same when the base original image 
signal S^^g' has been decomposed by use of wavelet transformation 
functions having second coefficients of filter. 

Figure 69 is a view showing an example of table data storing 
therein a space frequency characteristic representing the 
characteristic of enlargement or the characteristic of 
contraction, and 

Figures 7 OA and 7 OB are views showing the parameters 
employed in a Gaussian function representing the space frequency 
characteristic which in turn represents the characteristic of 
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enlargement or the characteristic of contraction. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

An image processing system in accordance with a first 
embodiment of the present invention will be described, 
hereinbelow. The image processing system is for carrying out 
a frequency enhancement processing by use of unsharp image 
signals on an original image signal obtained by reading out a 
radiation image of a human body recorded on a stimulable phosphor 
sheet so that an image suitable for diagnosis can be obtained. 
An image reproduced on the basis of the processed image signal 
is mainly recorded on photographic film and used in diagnosis. 

In Figure 1, an image processing system 1 in accordance 
with a first embodiment of the present invention comprises an 
unsharp image signal making means 2 which makes a plurality of 
unsharp image signals from an original image signal S^^.^ which 
is input from an image read-out apparatus or the like and has 
a predetermined picture element dens ity , and a non-linear 
processing means 3 which carries out a frequency enhancement 
processing for enhancing a particular frequency and obtains a 
processed image signal Sp^^^. The image processing system 1 
further comprises a transformation function defining means 4 
and a picture element density information input means 5. The 
transformation function defining means 4 is a means which defines 
transformation functions which the non-linear processing means 
3 uses in the transformation processing, for instance, by 
determining parameters such as inclinations of the functions 



and degrees of non-linearity on the basis of picture element 
density information M input from the picture element density 
information input means 5 and defining the functions on the basis 
of the parameters. The picture element density information 
input means 5 is a means for obtaining picture element density 
information M on the original image signal S^^g. The picture 
element density information M may be input as a value by the 
operator by a keyboard or by the operator selecting a value from 
a plurality of picture element densities displayed on a control 
screen. Further, the picture element density information M may 
be attached to the original image signal S^^.^ in the read-out 
apparatus and may be detected by the image processing system 
1. The picture element density information M may be input in 
any way so long as the transformation function defining means 
4 can recognize the picture element density of the original image 
signal S^^^. 

Making the unsharp image signals will be first described, 
hereinbelow. As shown in Figure 2, the unsharp image signal 
making means 2 comprises first to n-th filtering means 10. The 
first filtering means 10 carries out a filtering processing on 
the original image signal S^^g in x- and y-directions and makes 
a low resolution image signal Bl which is lower than the original 
image signal S^^g in resolution. The second filtering means 10 
carries out a similar filtering processing on the low resolution 
image signal Bl thus obtained and makes a low resolution image 
signal B2 which is lower than the low resolution image signal 
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Bl in resolution, and the third filtering means 10 carries out 
a similar filtering processing on the low resolution image signal 
B2 thus obtained and makes a low resolution image signal B3 which 
is lower than the low resolution image signal B3 in resolution. 
In this manner, the n-th filtering means 10 carries out a similar 
filtering processing on the low resolution image signal Bn- 
1 and makes a low resolution image signal Bn which is lower than 
the low resolution image signal Bn-1 in resolution. The unsharp 
image signal making means 2 further comprises first to n-th 
interpolation means 11. The interpolation means 11 carry out 
an interpolation/enlargement processing on the low resolution 
image signals Bl to Bn obtained at the respective stages of 
filtering, and makes a plurality of unsharp image signals S^^l 
to S^gN which are different in sharpness. 

In this particular embodiment, filters which 
substantially conform to one-dimensional Gaussian distribution 
are used in the filtering processing. That is, coefficients of 
the filters are determined according the following formula (4) 
which relates to a Gaussian signal. 

/(/) = /^-^ (4) 

This is because the Gaussian signal is good in localization in 
both a real space and a frequency space. For example, the 7x1 
one-dimensional filters are as shown in Figure 3A when a=l in 
formula ( 4 ) . 

The filtering is carried out on the original image signal 



S^rg or the low resolution image signals every second picture 
element as shown in Figure 4. When such filtering is carried 
out in both x- and y- directions, the number of picture elements 
in obtained low resolution image signal is reduced to 1/4 of 
the preceding image signal (e.g., the original image signal S^^^ 
in the case of the low resolution image signal Bl, and the low 
resolution image signal Bl in the case of the low resolution 
image signal B2 ) . That is, the number of picture elements in 
each of the low resolution image signals Bk(k stands for 1 to 
n) is 1/2^^ of the original image signal S^^.^, 

The interpolation/enlargement processing to be carried 
out on the low resolution image signals Bk thus obtained will 
be described, hereinbelow. Though various methods of 
interpolation such as B-spline can be employed, a Gaussian signal 
is also employed in the interpolation in this particular 
embodiment since low-pass filters based on Gaussian signals are 
employed in the filtering processing. Specifically, an 
approximation, 0=2*^"^ is employed in the following formula (5). 

/(r) = 2*a-/2^' (5) 

When interpolating the low resolution image signal Bl, 
a=l since k=l. In this case, a one-dimensional filter of 5x1 
such as shown in Figure 5 is employed in the interpolation. In 
this interpolation, a picture element of 0 in value is 
interpolated every second picture element in the low resolution 
image signal Bl, whereby the low resolution image signal Bl is 
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enlarged to a size equal to the size of the original image, and 
the interpolated low resolution image signal Bl is subjected 
to a filtering processing using the one-dimensional filter shown 
in Figure 5 . 

This interpolation/enlargement processing is carried out 
on all the low resolution image signals Bl to Bn. When 
interpolating a low resolution image signal Bk, a filter which 

is 3x2^-1 in length is prepared according to formula (5) and 2^-1 
picture elements of 0 in value are interpolated between each 
pair of adjacent picture elements, whereby the low resolution 
image signal Bk is enlarged to a size equal to the size of the 
original image signal. Then the interpolated low resolution 
image signal Bk is subjected to a filtering processing using 
the filter which is 3x2^-1 in length. 

The frequency enhancement processing has been generally 
based on the assumption that the picture element density at which 
the original image signal S^^g is obtained is constant, and the 
processed image signal S^^^^ has been made without taking into 
account the picture element density. However, in a practical 
image processing system, various image signals are input and 
all the image signals are not the same in picture element density. 
For example, in the radiation image processing system of this 
embodiment, the picture element density differs according to 
the size of the stimulable phosphor sheet, and is 5 lines /mm 
in a half size, 6.7 lines /mm in a quarter size and 10 lines/mm 
in a sixmo size. Further, the user can freely set the picture 
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element density. 

When image signals which are different in picture element 
density and accordingly in Nyquist frequency are subjected to 
a filtering processing using the same low-pass filters and to 
the same interpolation processing, band-limited signals 
obtained from unsharp signals S,,k in a manner to be described 
later become different from each other in frequency response 
characteristic according to the picture element density. For 
example, when band-limited signals are made from unsharp image 
signals S^^k obtained by carrying out a filtering processing 
using a one-dimensional filter shown in Figure 3A on original 
image signals which are respectively read at 10 lines /mm, 5 
lines/mm and 6.1 lines/mm, the band-limited signals from the 
original image signal at 6.7 lines /mm quite differs in frequency 
bands from those from the original image signals at 10 lines /mm 
and 5 lines /mm though the latter band-limited signals are the 
same in frequency bands since the original image signal 5 
lines/mm is a half of that at 10 lines/mm in Nyquist frequency. 
In this embodiment, the filtering processing is carried out using 
the filter shown in Figure 3A on original image signals at 10 
lines /mm and 5 lines /mm and using the filter shown in Figure 
3B on original image signals at 6 . 7 lines /mm so that the frequency 
bands of the obtained band-limited signals become substantially 
the same irrespective of the picture element density of the 
original image signal. The unsharp image signal making means 
2 switches the filters on the basis of the picture element density 
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information M input from the picture element density information 
input means 5 . 

The non-linear processing to be carried out on the unsharp 
image signals S^.k (k stands for 1 to n) thus obtained will be 
5 described, hereinbelow. Figure 6 shows a system for carrying 
out a frequency enhancement processing (an example of the 
non-linear processing) together with the unsharp image signal 
making means 2. As shown in Figure 6, unsharp image signals S^^k 
are made from the original image signal S^^g ^i^d differences 
10^^^ between the original image signal S^^g and the unsharp image 
Jl signals S^^k are taken by respective subtractors 21, whereby 
E band-limited signals (S,,g-S,,l, S,,l-S,,2, and the like) which are 
components in limited frequency bands of the original image 
L:. signal S^^g are made, 
isf ' The band-limited signals are transformed to predetermined 

i amplitudes with different transformation functions f^ to f^ in 
respective transformation circuits 22 and the transformed 
band-limited signals are added up and further added to the 
original image signal S^^g by an operator 23 according to the 
20 following formulae (2), whereby a processed image signal Sp^^^ 

in which a desired frequency components is enhanced to a required 
degree is generated. 

Sp.oc=S,rg+P(Sorg) ^ F,3„.(S.,g, S.^l, S,,2 , S^^N ) 

F.sm(Sorg/ S,,l, S,32, S,,N) 

25 =fi(So,g-S,,l)+f2(S,3l-S,32)+ 

+ f,(S,,k-l-S,3k)+ - +fN(S.3N-l-S,3N) (2) 
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wherein Sp^^^ is a processed image signal in which the high- 
frequency components are enhanced, S^^^ is an original image 
signal, S^gk (k=l to N) are unsharp image signals, f,^(k=l to N) 
is a transformation function, and p(Sorg) is a coefficient of 
enhancement determined on the basis of the original image signal . 

As the transformation functions fj,, non-linear functions 
represented by the following formula (6) are employed, 
f ( S i J =S,,,=S Yx { exp { X / S i J - 1 } / { exp ( X / Si J + 1 } ( 6 ) 

wherein S^^ represents an input signal, S^^^ represents an output 
signal, X represents a parameter which governs the degree of 
non-linearity, and Y represents a parameter which controls the 
inclination of the overall function. By adjusting the 
parameters X and Y, the frequency response characteristic of 
the band-limited signal can be changed. 

The processed image signal Sp,,^,^ is thus made. The problem 
to be solved by this embodiment and the manner of solving the 
problem in this embodiment will be described through an example, 
hereinbelow. When the filtering processing is carried out on 
the original image signal S^^^g using a filter which is determined 
according to the picture element density of the original image 
signal S^^gf thereby obtaining unsharp image signals S^^k, and 
band-limited signals are obtained from the unsharp image signals 
S^gk as disclosed in Japanese Unexamined Patent Publication No. 
10 (1998) -63837, the band-limited signals can be in the same 
frequency bands irrespective of the picture element density of 
the original image signal S^^g* However, in this case, the 



frequency response characteristics of the band-limited signals 
differ according to the picture element density of the original 
image signal S^^g- 

Figure 7 shows the frequency response characteristics of 
5 band-limited signals {S^^^-S^,!, S^,1^S^^2, S^^2~S^^3 , S^,3-S,,4, 

Sus^-S^gS, S^g5-S^g6) obtained by carrying out ^ on an original image 
signal at a picture element density of 10 lines/mm, the filtering 
processing using the one-dimensional filter shown in Figure 3A. 
Figure 8 shows the frequency response characteristics of 
10 J^^ band-limited signals (S^^^-S^^l, 8,31-8^,2, 8^,2-8^,3, S^^S-S^^Ar 
2: 8^34-8^^5, S^g5-S^,s6) obtained by carrying out, on an original image 
signal at a picture element density of 6.7 lines/mm, the 
if filtering processing using the one-dimensional filter shown in 
Figure 3B. Figure 9 shows the frequency response 

15 [7 characteristics of band-limited signals (S^^g-S^^l, 8^^1-8^32, 
J; S,32-S,,3, 8,,3-8,,4, 8,,4-S,,5, S,,5-8,,6 ) obtained by carrying out , 
on an original image signal at a picture element density of 5 
lines/mm/ the filtering processing using the one-dimensional 
filter shown in Figure 3A. As can be seen from Figures 7 to 9, 
20 though the maximum frequency band of the band-limited signals 

can be determined according to the picture element density, the 
peak frequencies of the band-limited signals obtained from the 
respective original image signals other than the maximum 
frequency band are reduced 1/2 by 1/2, e.g. , 1 (only in the case 
25 of picture element density of 10 lines/mm) , 0.5, 0.25, 0.125 

That is, the peak frequency of k-th band-limited signal is 



57 



1/2 of (k-l)-t:h band-limited signal for each picture element 
density. Whereas, the response at each peak frequency, that is, 
energy of each band-limited signal, differs according to the 
picture element density. Accordingly, if the same parameters 
X and y (formula (6)) are used for different picture element 
densities, the frequency response characteristic of the 
obtained processed image signal Sp^^^ differs according to the 
picture element density of the original image signal S^^g. 

In accordance with the present invention, in order to 
overcome this problem, the transformation function defining 
means 4 sets the parameters X and Y of the transformation function 
on the basis of the picture element density information M input 
from the picture element density information input means 5 and 
the band-limited signals are transformed according to the 
transformation function which is defined according to the 
picture element density of the original image signal S^^g. Figure 
10 shows in detail the structure of the transformation function 
defining means 4. As shown in Figure 10, the transformation 
function defining means 4 comprises a parameter setting means 
7 and a memory 8 in which combinations of parameters X and Y 
for each picture element density are stored. 

An example of combination of parameters X and Y is shown 
in the following table 1. In the combinations of parameters X 
and Y shown in table 1, the parameters X and Y for the picture 
element densities of 6.7 lines/mm and 5 lines/mm are set so that 
a processed image signal S^^^^ having a frequency response 
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substantially equal to that obtained by the parameters X and 
Y for the picture element density of 10 lines/mm can be obtained. 
In table 1, the parameters X and Y for each picture element 
density on higher lines are for band-limited signals in higher 
frequency bands, and since no band-limited signal having a 
frequency peak at 0.03 cycles /mm exists for the picture element 
densities of 6.7 lines /mm and 5 lines/mm, no parameter is given. 



Tablg 1 



10 lines /mm 


6.7 lines /mm 


5 lines /mm 


X 


Y 


X 


Y 


X 


Y 


15 


1.00 


20 


1.00 


20 


1.00 


20 


0.99 


30 


0.83 


30 


0.81 


30 


0.80 


40 


0.67 


40 


0.66 


40 


0.66 


50 


0.43 


50 


0.42 


50 


0.41 


60 


0.32 


60 


0.29 


60 


0.29 











The parameter setting means 7 selects the combinations 
of parameters X and Y stored in the memory 8 according to the 
picture element density represented by the picture element 
density information M input from the picture element density 
information input means 5, and the selected combinations of 
parameters X and Y are input into the transformation function 
defining means 4. The transformation function defining means 
4 defines the transformation functions on the basis of the 
selected combinations of parameters X and Y, and the non-linear 
processing means 3 carries out on the band-limited signals 
according to formulae (2) on the basis of the transformation 
functions thus defined. 



Figure 11 shows the frequency response characteristic of 
the processed image signal Sp^„^ obtained by processing the 
band-limited signals obtained from an original image signal at 
10 lines/mm using the transformation functions defined in the 
manner described above. Figure 12 shows the frequency response 
characteristic of the processed image signal Sp^^^ obtained by 
processing the band-limited signals obtained from an original 
image signal at 6.7 lines/mm using the transformation functions 
defined in the manner described above. Figure 13 shows the 
frequency response characteristic of the processed image signal 
Spr^^ obtained by processing the band-limited signals obtained 
from an original image signal at 5 lines/mm using the 
transformation functions defined in the manner described above. 
As can be seen from Figures 11 to 13, when the frequency 
enhancement processing is carried out by use of transformation 
functions defined by the parameters X and Y shown in table 1, 
the frequency response characteristic of the obtained processed 
image signal Sp^^^ is substantially constant irrespective of the 
picture element density of the original image signal. 

The parameters X and Y may be set, for instance, by the 
method disclosed in Japanese Unexamined Patent Publication No. 
10(1998)-63838. This method will be described, hereinbelow. 
AS a system for setting the parameters X and Y, a system which 
prompts the user to enter a desired frequency response 
characteristic by way of a display and recognizes through the 
input the frequency response characteristic which the user 
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desires. Specifically, the system comprises software, an input 
system, a display and the like. 

When setting the parameters X and Y, a screen 7a (Figure 
14) which shows the frequency response characteristic curve of 
5 the processed image signal Sp^^^ obtained when the band-limited 

signals are processed by the current transformation functions 
(initial transformation functions). On the characteristic 
curve, six movable points are displayed. That is, in the case 
where the image processing is a non-linear processing using n 
10 y band-limited signals, the number of the movable points is n- 
5" The user moves necessary point (s) to a position 

corresponding to a desired frequency response characteristic 
S by operating a pointing device such as a mouse. Screen 7b in 
Figure 14 shows this procedure. This is done so that the 
15 f"' responses in the respective frequency bands become equal to each 
% other. After the necessary points are moved, the new points are 
i^J recognized as to indicate desired frequency responses at the 
frequencies designated by the points , and the frequency response 
characteristic curve is recalculated and a new characteristic 
20 curve which passes through the new points is displayed. Screen 

7c in Figure 14 shows this state. When the new characteristic 
curve is a desired one, the user selects the "OK" button shown 
in screen 7d and otherwise selects the "Cancel" button to redo 
setting parameters. The parameters may be set in other various 
25 ways. For example, the parameter setting means may be arranged 
so that a plurality of values of frequency are shown on the screen 
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in sequence and the user inputs a desired frequency response 
characteristic for the shown frequency. 

When instruction that the parameters are to be set is made, 
the input desired frequency response characteristics al to a6 
(in the case where the image processing to be carried out using 
six band-limited signals) are substituted, for instance, in the 
following simultaneous equations (7). 

al=YlxSll+Y2xS21+Y3xS31-hY4xS41+Y5xS51-fY6xS61 
a2=YlxSl2+Y2xS22H-Y3xS32+Y4xS42+Y5xS52+Y6xS62 
a3=YlxSl3+Y2xS23+Y3xS33+Y4xS43-fY5xS53+Y6xS63 
a4=YlxSl4+Y2xS24+Y3xS34+Y4xS44+Y5xS54+Y6xS64 
a5=YlxS15+Y2xS25+Y3xS35+Y4xS45+Y5xS55+Y6xS65 
a6=YlxSl6+Y2xS26+Y3xS36+Y4xS4 6+Y5xS56+Y6xS66 
Yl'-Ye^O (7) 
wherein al to a6 represent the desired frequency response 
characteristics, Sll to S66 represent the frequency response 
characteristics of the band-limited signals made on the basis 
of the original image signal, and Yl to Y6 represent parameters 
representing the inclinations of the transformation functions 
near 0. The degree of the equations is determined according to 
the number of the band-limited signals and the above equations 
are given only as an example. 

By solving the simultaneous equations (7) with the desired 
frequency response characteristics al to a5 and the frequency 
response characteristics of the band-limited signals Sll to S66 
taken as known values and the parameters Yl to Y6 taken as 
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variables, the parameters Yl to Y6 are obtained, and the 
transformation functions are defined on the basis of the 
parameters. The parameters X representing the degree of 
non-linearity may be set to any values corresponding to the 
parameters Yl to Y6 and, for instance, may be set to fixed values 
which are determined according to the picture element densities 
as shown in the aforesaid table 1 , The parameters X and Y 
determined picture element density by picture element density 
are stored in the memory 8 . 

The values of the frequency response characteristics of 
the band-limited signals 811 to 866 are for six band-limited 
signals and 811 to 816 (821 to 826, — 861 to 866) represent the 
frequency response characteristics at six frequencies of one 
band-limited signal . The values of the frequency response 
characteristics of the band-limited signals 811 to 866 can be 
obtained by a method of transformation normally used in frequency 
analysis on the basis of the coefficients of filter and the 
coefficients of interpolation used in making the unsharp image 
signals . 

The parameters X and Y may be set by trial and error. 

Operation of the image processing system 1 of this 
embodiment will be described with reference to the flow chart 
shown in Figure 15. An original image signal S^^^ is first input 
into the image processing system 1, for instance, from an image 
read-out apparatus, (step 81) The original image signal S^^g 
is input into the unsharp image signal making means 2 and unsharp 



image signals S^^k are made according to the picture element 
density of the original image signal S^^g. (step S2) Picture 
element density information M on the original image signal S^^g 
is input into the picture element density information input means 
5 . ( step S3 ) Then the combinations of the parameters are 
selected in the transformation function defining means 4 on the 
basis of the picture element density information M in the manner 
described above, and the transformation functions are defined 
on the basis of the selected combinations of the parameters • 
(step S4) Steps S3 and S4 may be executed before steps SI and 
S2. When steps S3 and S4 and steps SI and 32 are executed in 
parallel, the operation time can be shortened. Frequency- 
limited signals are made on the basis of the unsharp image signals 
S^gk and the non-linear processing (frequency enhancement 
processing) represented by the aforesaid formulae (2) is 
executed by the non-linear processing means 3 on the basis of 
the transformation functions defined by the transformation 
function defining means 4, whereby a processed image signal Sp^^^ 
is obtained, (step S5) 

The frequency enhancement processing is for making add 
signals, which are to be added to the original image signal in 
order to enhance desired frequency components, by use of the 
unsharp image signals. The frequency bands are differently 
processed so that the signals in the respective frequency bands 
forming the respective add signals become desired signals in 
order to prevent generation of artifact due to the frequency 
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enhancement processing. In order to make the desired signals, 
it is preferred that signals have desired frequency response 
characteristics irrespective of the picture element density of 
the original image signal S^^g. In the image processing system 
5 1 of this embodiment, the parameters of the transformation 

functions are defined so that the frequency response 
characteristic of the processed image signal Sp^^^ become 
substantially constant irrespective of the picture element 
density of the original image signal S^^g, and accordingly, a 
10 O processed image signal Sp^^^ having substantially constant 
4^ frequency response characteristic can be obtained irrespective 
flJ of the picture element density of the original image signal S^^g. 

Though, in the embodiment described above, the frequency 
- enhancement processing represented by the aforesaid formulae 
15 (2) is carried out as a non-linear processing, a dynamic range 

4- compression processing and a frequency enhancement processing 
O may be executed simultaneously as a non-linear processing. 
Figure 16 is a view showing a system for executing a frequency 
enhancement processing and a dynamic range compression 
20 processing as a non-linear processing in an image processing 

system in accordance with a second embodiment of the present 
invention. As shown in Figure 16, differences between two of 
the original image signal S^^g and the unsharp image signals S^^k 
made by the filtering means 10 and the interpolation means 11 
25 are taken, whereby a plurality of band-limited signals (e.g., 

S-S^^lr S^3l-S^32 and the like) respectively representing 
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frequency components in limited frequency bands of the original 
image signal are made. The filters used in the filtering means 
19 are determined according to the picture element density of 
the original image signal in the manner similar to that described 
above. The band-limited signals thus obtained are input into 
first and second transformation sections 3a and 3b and processed 
by respective transformation means 22a and 22b in the first and 
second transformation sections 3a and 3b. 

Transformation by the transformation means 22a in the 
first transformation section 3a is executed by use of the 
transformation functions defined on the basis of the picture 
element density information M on the original image signal S^^^ 
as described above. For example^ the transformation functions 
shown in Figure 17 or 18 or combinations of the transformation 
functions shown in Figures 17 and 18 are used. These 
transformation functions are obtained by setting the parameters 
so that the frequency response characteristic of the processed 
image signal Sp^^c becomes substantially constant irrespective 
of the picture element density of the original image signal S^^g. 

In accordance with the transformation functions shown in 
Figure 17 , the band-limited signals are transformed so that those 
which are large in amplitude are suppressed and the degree of 
suppression is increased as the frequency band of the band- 
limited signal becomes higher. This is for taking into account 
the fact that higher frequency components contained in edges 
of an actual radiation image is smaller in amplitude than low 



66 



frequency components. In actual radiation images, even a 
substantially sharp edge is not in the form of a correct step 
and the amplitude often becomes smaller as the frequency becomes 
higher. Accordingly, it is preferred that suppression be made 
from a smaller amplitude as the frequency of the band-limited 
signals becomes higher. The functions shown in Figure 17 serve 
for this purpose. 

The transformation functions shown in Figure 18 are for 
transforming the band-limited signals to those having values 
not larger than absolute values of the band-limited signals, 
which values are determined on the basis of the absolute values 
of the band-limited signals* As the frequency of the frequency 
band to be processed by the function becomes lower, the absolute 
value of a transf oimied image signal obtained by transforming 
a band-limited signal whose absolute value is in a predetermined 
range near 0 becomes smaller. That is, the functions all pass 
through the origin and all have inclinations smaller than 1. 
Further, the inclination near 0 is smaller as the frequency of 
the frequency band to be processed by the function becomes lower. 
When an add signal obtained by adding up the transformed image 
signals is added to the original image signal S^^g, these 
functions contribute to smoothen the joint between the original 
image signal S^^.^ and the add signal, that is, rise of the signal. 

Similarly transformation by each transformation means 22b 
in the second transformation section 3b is executed by use of 
the transformation functions shown in Figure 19 or 18 or 



combinations of the transformation functions shown in Figures 
19 and 18. The transformed image signals output from the first 
and second transformation sections 3a and 3b are respectively 
input into operators 23a and 23b. The operator 23a executes 
operation for making signals necessary for the frequency 
enhancement processing, and the operator 23b executes operation 
for making signals necessary for the dynamic range compression 
processing. 

The operator 23a executes a frequency enhancement 
processing similar to that described above with reference to 
Figure 6. That is, the band-limited signals transformed by the 
first transformation section 3a are added up, and an add signal 
thus obtained is multiplied by a coefficient of enhancement p 
which is determined according to the value of the original image 
signal S^^^. 

The operator 23b executes a dynamic range compression 
processing in the following manner. The band-limited signals 
transformed by the second transformation section 3b are added 
up, and an add signal thus obtained is subtracted from the 
original image signal S^^g. Then the differential signal thus 
obtained is transformed on the basis of a transformation function, 
thereby obtaining a coefficient of dynamic range compression. 
The signals respectively obtained by the operators 23a and 23b 
are added to the original image signal S^^^by an adder 28, whereby 
a processed image signal Sp^^^ is obtained. 

The processing described above is represented by the 
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following formula (3), 
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10 C3 wherein Sp^^^ is a processed image signal, S^^g is an original image 
signal, S^^k (k=l to N) is an unsharp image signal, fyj,(k=l to 
nJ N) is a transformation function used in the first transformation 
=-0 section, fdj,(k=l to N) is a transformation function used in the 
^ second transformation section, P(S^3.g) is a coefficient of 
15 enhancement determined on the basis of the original image signal, 

£ and D{S,^g-Fd^, (3,^^, S,,l , S^^2 , S^,N) } is a coefficient of dynamic 

Q range compression determined on the basis of the low frequency 
component signal, D being a function for transforming D{S^^g- 

Fdrc(Sorg/ S^^l , S^^2 , S^^N ) } . 

20 In the image processing system shown in Figure 16, the 

add signal used in the frequency enhancement processing and that 
used in the dynamic range compression processing are separately 
obtained on the basis of the original image signal S^^g and are 
finally added. However, it is possible to first carry out on 

25 one of the frequency enhancement processing and the dynamic range 



compression processing on the original image signal S^^g and to 
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carry out the other processing on the resultant signal. 
Generally the lower density portion of a radiation image includes 
a relatively large amount of noise since the lower density 
portion is exposed to less radiation upon taking the radiation 
5 image. Since the frequency enhancement processing is a 

density-dependent processing where the degree of enhancement 
is increased with increase in density, the low density portion 
is not enhanced and accordingly noise is not enhanced when the 
frequency enhancement processing is directly carried out on the 
10 O original image signal. To the contrast, when the frequency 
4^ enhancement processing is carried out after the dynamic range 
nJ compression processing is carried out on the original image 
signal, the low density portion is increased in its density by 
1 the dynamic range compression processing, and accordingly, the 
15 low density portion is enhanced by the following frequency 

4- enhancement processing, whereby noise is also enhanced. 
Q Accordingly, it is preferred that each signal be made on the 
basis of the original image signal S^^g. Further also from the 
viewpoint of saving the processing time, it is preferred that 
20 the two processings be executed in parallel. 

By defining the parameters for the transformation 
functions so that the frequency response characteristic of the 
processed image signal Sp^,^^ becomes constant irrespective of the 
picture element density of the original image signal S^^g, the 
25 processed image signal Sp^^^ can have a substantially constant 

frequency response characteristic irrespective of the picture 
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element density of the original image signal S^^^. 

In the second embodiment described above, when the 
parameters of the transformation functions are set as shown in 
the following table 2 in the case where the picture element 
density of the original image signal is 10 lines /mm, the 
frequency response characteristic of the processed image signal 
Sproc becomes as shown in Figure 20. As shown in Figures 7 to 
9, in the case of an original image signal at 10 lines /mm, there 
exists a band-limited signal corresponding to 1 cycle/mm, but 
in the case of an original image signal at 6.7 lines /mm or at 
5 lines /mm, there exists no band-limited signal corresponding 
to 1 cycle/mm. Accordingly, in the case of original image 
signals at 6.7 lines /mm and at 5 lines/mm, the frequency response 
characteristics of the processed image signals S^^^^ cannot 
conform to the processed image signal Sp^^^ for the original image 
signal at 10 lines /mm over the entire frequency bands. In an 
image obtained by reproducing a processed image signal, the low 
frequency components are visually more prominent than the high 
frequency components. Accordingly, in such a case, it is 
preferred that the transformation functions for band-limited 
signals at 6.7 lines/mm and at 5 lines /mm be defined by the 
parameters X and Y shown in table 2 so that the high frequency 
components not lower than 1 cycle/mm conform to those of the 
band-limited signals at 10 lines /mm to a certain extent and the 
low frequency components lower than 1 cycle/mm conform to those 
of the band-limited signals at 10 lines /mm. When the parameters 



of the transformation functions are set as shown in the following 
table 2 in the case where the picture element densities of the 
original image signals are 6.7 lines/mm and 5 lines/mm, the 
frequency response characteristics of the processed image 
signals Sp^^^ become as shown in Figures 21 and 22. As shown in 
Figures 21 and 22, when original image signals at 6.7 lines /mm 
and 5 lines /mm are processed by use of transformation functions 
defined by the parameters shown in table 2, the frequency 
response characteristics of the processed image signals Sp^^^ well 
conform to the processed image signal Sp^^^ for the original image 
signal at 10 lines/mm as for the low frequency components lower 
than 1 cycle/mm though less conform to the processed image signal 
Sproc for the original image signal at 10 lines /mm as for the high 
frequency components not lower than 1 cycle/mm. Accordingly, 
a processed image signal Sp^^c the frequency response 
characteristic of which is substantially constant irrespective 
of the picture element density of the original image signal at 
least for the low frequency components, which are visually more 
prominent than the high frequency components, can be obtained. 



Ta b le 2 



10 lines /nun 


6.7 lines /nun 


5 lines/nim 


X 


y 


X 


Y 


X 


Y 


15 


1.00 


20 


0.93 


20 


0.92 


20 


0.83 


30 


0.44 


30 


0.42 


30 


0.42 


40 


0.37 


40 


0.36 


40 


0.36 


50 


0.18 


50 


0.17 


50 


0.17 


60 


0.13 


60 


0.12 


60 


0.12 











Though, in the embodiments described above, the original 



image signal S^^g is an image signal read out from a stimulable 
phosphor sheet, the original image signal S^^gmay be any image 
signal provided that the resolution of the image represented 
by the image signal, e.g., the relation between the size of the 
5 object and the sampling intervals (dpi), is known. 

Further, though in the embodiments described above, 
non-linear functions are used and a non-linear processing is 
carried out on the band-limited signals, linear functions or 
constants may be employed as the transformation functions. In 
10 Q such a case, the parameters of the linear functions or the values 
4^ of the constants are set according to the picture element density 
fij of the original image signal. 

4^ Further, though, in the embodiments described above, the 

. unsharp image signals are obtained from the original image signal 
15 Uh by filtering and interpolation/enlargement and the band-limited 
4: signals are obtained from the original image signal and the 
Ej unsharp image signals, the band-limited signals may be made, 
for instance, by transforming the original image signal to 
multiple resolution image signals by a wavelet transformation 
20 or a Laplacian pyramid, making the unsharp image signals from 
the image signals at the respective resolutions and making the 
band-limited signals from the unsharp image signals. 

Further, though, in the embodiments described above, the 
parameter setting means 7 shown in Figure 10 is arranged to set 
25 the parameters X and Y on the basis of the original image signal 
and the picture element density thereof each time an original 
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image signal S^^g is input, the parameter setting means 7 may be 
arranged to determine the parameters X and Y for the original 
image signal S^^g by comparing the picture element density of the 
original image signal S^^^ with reference picture element 
5 densities which have been stored in the memory 8 together with 
the parameters X and Y for the respective picture element 
densities and by setting, as the parameters X and Y for the 
original image signal S^^g, the parameters X and y for the 
reference picture element density closest to the picture element 
10 fi density of the original image signal S^^^. The reference picture 
element density closest to the picture element density of the 
nj original image signal S^^g niay be determined, for instance, in 
,|| the following manner. That is, assuming that three reference 
^ picture element densities stored in the memory 8 together with 
15 the parameters X and Y for the respective reference picture 
'jz element densities are ^, h and q, (a<b<c) and the picture element 
^ density of the original image signal S^^g to be processed is x, 
the reference picture element density is q when x>b+(c-b)/2, 

is h when b+(c-b) /2^x>a+(b-a) /2 , and is a when a-l-(b-a) /2s:x. 

20 Specifically the reference picture element densities may be at 

least two of 10 lines /mm, 6.7 lines/mm and 5 lines/imti. In this 
case, when the picture element density of the original image 
signal S^^g to be processed is higher than 6 . 7+( 10-6 . 7 ) /2=8 . 35 , 
the combinations of the parameters X and Y for the picture element 

25 density of 10 lines/mm shown in table 1 or table 2 are set as 

the parameters X and Y for the original image signal S . When 
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the picture element density of the original image signal S^^gto 
be processed is not higher than 6 . 7+{ 10-6 • 7 ) /2=8 . 35 and higher 
than 5+(6,7-5)/2=5*85, the combinations of the parameters X and 
Y for the picture element density of 6.7 lines/imn shown in table 
1 or table 2 are set as the parameters X and Y for the original 
image signal S^^g. When the picture element density of the 
original image signal S^^g to be processed is not higher than 
5+ (6. 7-5)72=5.85, the combinations of the parameters X and Y 
for the picture element density of 5 lines/mm shown in table 
1 or table 2 are set as the parameters X and Y for the original 
image signal S^^g. 

It is preferred that the parameters selected be stored 
related to the processed image signal Sp^^^ and/or the picture 
element density of the original image signal S^^-g. With this 
arrangement, when an image is output (output as a hard copy on 
photographic film or displayed on a CRT or the like) another 
time, an image at a quality equivalent to that of a previously 
output image can be obtained. The parameters selected and the 
picture element density of the original image signal S^^g may be 
stored either in the same storage medium or different storage 
media . 

An image processing system in accordance with a third 
embodiment of the present invention will be described, 
hereinbelow. An image processing system 101 in accordance with 
the third embodiment of the present invention is also for 
carrying out a frequency enhancement processing by use of unsharp 



image signals on an original image signal obtained by reading 
out a radiation image of a human body recorded on a stimulable 
phosphor sheet so that an image suitable for diagnosis can be 
obtained* An image reproduced on the basis of the processed 
image signal is mainly recorded on photographic film and used 
in diagnosis. 

As shown Figure 23, the image processing system 101 in 
accordance with the third embodiment of the present invention 
is substantially the same as the image processing system 1 in 
accordance with the first embodiment of the present invention* 
Accordingly, in Figure 23, the elements analogous to those shown 
in Figure 1 are given the same reference numerals and will not 
be described in detail here. The image processing system of this 
embodiment comprises an unsharp image signal making means 2, 
a processing means 3 which carries out a frequency enhancement 
processing for enhancing a particular frequency and obtains a 
processed image signal Sp^^^, a transformation function 
calculating means 104 and a resolution information input means 
105 . The image processing system of the third embodiment mainly 
differs from that of the first embodiment in arrangement of the 
transformation function calculating means 104. That is, in this 
embodiment, transformation function calculating means 104 
calculates the transformation functions, which the processing 
means 3 uses in the frequency enhancement processing, by 
correcting reference transformation functions which are 
employed to carry out the frequency enhancement processing on 
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an original image signal having a known reference picture element 
density, according to the resolution of the original image signal 
S^rg represented by the resolution information M input from the 
resolution information input means 105, as will become apparent 
later* The transformation functions for the original image 
signal have been determined and stored in advance. 

Unsharp image signals S^^k are made in the same manner as 
described above in conjunction with Figure 2 and the frequency 
enhancement processing is carried out by use of the unsharp image 
signals S^^k in the same manner as described above in conjunction 
with Figure 6 . 

In this embodiment, the transformation functions f,, (k 
stands for 1 to N) for transforming the band-limited signals 
are determined in the following manner. In this particular 
embodiment, the transformation functions f^ are constants and 
the reference resolution is 10 lines/mm. Figure 24 shows the 
frequency response characteristics of the band-limited signals 
obtained from an original image signal S^^g having a resolution 
of 10 lines /mm. As shown in Figure 24, six band-limited signals 
are obtained from the original image signal S^^g- The original 
image signal S^^g having a resolution of 10 lines /mm is 5 cycles /mm 
in Nyquist frequency. The peak frequency of the band-limited 
signal in the highest frequency band falls on the Nyquist 
frequency. The peak frequency of the band-limited signal in the 
second highest frequency band is 1 . 0 cycle/mm, 1/5 of the Nyquist 
frequency. The peak frequency of the band-limited signal is 



reduced 1/2 by 1/2, e.g., 0.5 cycles/mm, 0.25 cycles/mm/ 0.12 
cycles /mm and 0.06 cycles /mm, as the frequency band of the 
band-limited signal becomes lower. In this particular 
embodiment, the transformation functions fj, (k stands for 1 to 
N) for transforming the band-limited signals obtained from the 
original image signal S^^g having a resolution of 10 lines /imn are 
as shown in the following table 3, 



Table 3 



peak frequency (cycles /mm) 


5 


1.0 


0.5 


0.25 


0.12 


0.06 


transformation function 


1.00 


0.90 


0.80 


0.60 


0.40 


0.20 



In order to enhance a desired frequency component of the 
original image signal S^^g, it is necessary to carry out the 
frequency enhancement processing by use of transformation 
functions fj, determined frequency by frequency of the band- 
limited signals. However, when the transformation functions 
for an original image signal S^^^ having a resolution of 10 
lines /mm are used, as they are, for frequency enhancement 
processing of an original image signal S^^g having a resolution 
of 300dpi, the frequency response characteristic of the 
processed image signal Sp^^^ differs from that of the processed 
image signal Sp^^^ obtained from the original image signal S^^g 
having a resolution of 10 lines /mm. In Figure 25, the solid line 
shows the frequency response characteristic of the processed 
image signal Sp^^^ obtained by carrying out the frequency 
enhancement processing by use of the transformation functions 
shown in table 3 on the original image signal S at 10 lines/mm 
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and the broken line shows that of the processed image signal 
Sp^^^ obtained by carrying out the frequency enhancement 
processing by use of the transformation functions shown in table 
3 on the original image signal S^^g at 300dpi. 

In order to overcome this problem, when an original image 
signal S^^g having a resolution different from 10 lines/mm 
(reference resolution) is inputs the transformation functions 
shown in table 3 (reference transformation functions) are 
corrected according to the resolution (300dpi in this particular 
embodiment) of the input original image signal S^^g. 

Figure 26 shows the frequency response characteristic of 
band-limited signals obtained from an original image signal S^^^ 
having a resolution of 300dpi. As can be seen from Figure 2S , 
the original image signal S^^g at 300dpi is 5.9 cycles /mm in 
Nyquist frequency, which falls on the peak frequency of the 
band-limited signal in the highest frequency band. The peak 
frequency of the band-limited signal in the second highest 
frequency band is 1 . 018cycle/mm, 1/5 of the Nyquist frequency. 
The peak frequency of the band-limited signal is reduced 1/2 
by 1/2, e.g., 0.59 cycles/mm, 0.30 cycles/imu/ 0.15 cycles/mm 
and 0.07 cycles /mm, as the frequency band of the band-limited 
signal becomes lower. This relation substantially conforms to 
the relation for the original image signal S^^g at 10 lines/mm 
shown in Figure 24, and holds irrespective of the resolution 
of the original image signal S^^g so long as the filters used in 
making the low-resolution image signals are not changed. 
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The transformation functions for transforming the 
band-limited signals obtained from an original image signal S^^.^ 
at 3 00dpi are obtained by correcting the reference 
transformation functions shown in table 3 on the basis of the 
peak frequencies of the band-limited signals obtained from the 
original image signal S^rg at 10 lines /mm (reference original 
image signal S^rgO) and those of the band-limited signals obtained 
from the original image signal S^^g at 300dpi. This correction 
can be made^ for instance, by interpolation or extrapolation. 
For example, the transformation function for the band-limited 
signals whose peak frequency is 1.18 cycles/mm is obtained by 
determining the relation between the peak frequencies and the 
values of the transformation functions on the basis of the values 
of the transformation functions at 5 cycles/mm and 1 cycle/imn 
in table 3, and substituting the value of the peak frequency 
(1.18) in the following formula (6) representing the determined 
relation between the peak frequencies and the values of the 
transformation functions. 
fx=l/40xC-h0.875 (6) , 

wherein fx=the value of the transformation function and C 
represents the peak frequency. Such linear interpolation is 
carried out on all the peak frequencies of the band-limited 
signals obtained from the original image signal S^^g at 300dpi, 
thereby obtaining the transformation functions for all the 
band-limited signals as shown in the following table 4. 
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Table 4 



peak frequency (cycles/mm) 


5.9 


1.18 


0.59 


0.30 


0.15 


0.07 


transformation function 


1.02 


0.91 


0.82 


0.64 


0.45 


0.24 



Though, in the description above, linear interpolation is used, 
other interpolations may be employed. 

In Figure 27, the solid line shows the frequency response 
characteristic of the processed image signal Sp^^^ obtained by 
carrying out the frequency enhancement processing by use of the 
transformation functions shown in table 3 on the original image 
signal S^^g at 10 lines/mm and the broken line shows that of the 
processed image signal S^^^^ obtained by carrying out the 
frequency enhancement processing by use of the transformation 
functions shown in table 4 on the original image signal S^^^ at 
300dpi. As can be seen from Figure 27, the processed image 
signals S^^^^ are substantially equal to each other in frequency 
response characteristic. 

Operation of the image processing system 101 of this 
embodiment will be described with reference to the flow chart 
shown in Figure 28. An original image signal S^^g is first input 
into the image processing system 101, for instance, from an image 
read-out apparatus, (step Sll) The original image signal S^^g 
is input into the unsharp image signal making means 2 and unsharp 
image signals S^^k are made by the unsharp image signal making 
means 2 • (step S12) Resolution information M on the original 
image signal S^^g is input into the resolution information input 
means 105. (step S13) Then the transformation function 
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calculating means 4 calculates the transformation functions on 
the basis of the resolution information M in the manner described 
above, (step S14) Steps S3 and S4 may be executed before steps 
SI and 32. When steps S3 and S4 and steps SI and 32 are executed 
in parallel,. the operation time can be shortened. 
Frequency-limited signals are made on the basis of the unsharp 
image signals S^^k and the frequency enhancement processing 
represented by the aforesaid formulae (2) is executed by the 
processing means 3 on the basis of the transformation functions 
calculated by the transformation function calculating means 4, 
whereby a processed image signal Sp^^^ is obtained. (step S15) 
The frequency enhancement processing is for making add 
signals, which are to be added to the original image signal in 
order to enhance desired frequency components,, by use of the 
unsharp image signals. The frequency bands are differently 
processed so that the signals in the respective frequency bands 
forming the respective add signals become desired signals in 
order to prevent generation of artifact due to the frequency 
enhancement processing. In order to make the desired signals, 
it is preferred that signals have desired frequency response 
characteristics irrespective of the resolution of the original 
image signal S^^^. In the image processing system 101 of this 
embodiment, the transformation functions are calculated so that 
the frequency response characteristics of the processed image 
signal Sp^^^ become substantially constant irrespective of the 
resolution of the original image signal S^^^, and accordingly. 
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a processed image signal S^^^^ having substantially constant 
frequency response characteristics can be obtained irrespective 
of the resolution of the original image signal S^^g. 

Further, since the transformation functions for original 
image signals S^^^ having resolutions different from the 
reference resolution are obtained by correcting the reference 
transformation functions, it is not necessary to prepare a 
plurality of transformation functions for various resolutions, 
whereby trouble to manage a plurality of transformation 
functions can be saved. 

Though, in the third embodiment described above, only one 
series of reference transformation functions for one reference 
resolution, 10 lines /mm, are prepared and transformation 
functions for other resolution are calculated on the basis of 
the one series of reference transformation functions, it is 
possible to prepare a plurality of series of reference 
transformation functions for a plurality of resolutions which 
are relatively frequently used and to calculate transformation 
functions for an object picture element density on the basis 
of one of the reference resolutions. For example, when first 
and second reference transformation functions for first and 
second resolutions are prepared in advance, transformation 
functions for an object original image signal to be processed 
may be calculated on the basis the first reference transformation 
functions when the resolution of the object original image signal 
is lower than the middle between the first and second resolutions 
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as shown by arrow A in Figure 29A and on the basis the second 
reference transformation functions when the resolution of the 
object original image signal is not lower than the middle between 
the first and second resolutions as shown by arrow B. Otherwise, 
the transformation functions for the object original image 
signal may be calculated on the basis the first reference 
transformation functions when the resolution of the object 
original image signal is not higher than the first resolution 
as shown by arrow C in Figure 29B and on the basis the second 
reference transformation functions when the resolution of the 
object original image signal is higher than the first resolution 
as shown by arrow D. 

Further, though in the third embodiment described above, 
the transformation functions are constants, the transformation 
functions may be non-linear functions having predetermined 
inclinations . 

An image processing system in accordance with a fourth 
embodiment of the present invention, where non-linear functions 
are employed as the transformation functions, will be described, 
hereinbelow. 

In the fourth embodiment, as the transformation functions 
f},, non-linear functions represented by the following formula 
( 8 ) are employed . 

f ( Sin ) =S,,,=Si,xYx { exp ( X/Si J -1 } / { exp ( X/Si J +1 } ( 8 ) 

wherein S^^ represents an input signal, S^^^ represents an output 

signal, X represents a parameter which governs the degree of 
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non-linearity, i.e., the condition of suppression, and Y 
represents a parameter which controls the inclination of the 
overall function, i.e., the frequency response characteristic. 
By adjusting the parameters X and Y, the frequency response 
characteristic of the band-limited signal can be changed. In 
this embodiment, parameters X and Y shown in the following table 
5 are employed. Table 5 shows combinations of parameters X and 
Y for the resolutions of, 10 lines /mm, 6 . 7 lines /mm and 5 lines /mm, 
and in table 5, the parameters X and Y for each resolution on 
higher lines are for band-limited signals in higher frequency 
bands. As shown in Figure 30, as the Nyquist frequency becomes 
lower, the number of the band-limited signals used in the 
non-linear processing is reduced. For example, in the case of 
an original image signal at 10 lines/mm, the number of the 
band-limited signals is 6, whereas in case of original image 
signals at 6.7 lines/mm and 5 lines /mm, the number of the 
band-limited signals is 5. This is because no band-limited 
signal having a frequency peak at 0.03 cycles /mm exists for the 
resolutions of 6.7 lines/imn and 5 lines /mm. The case where 10 
lines/mm is employed as the reference resolution, the 
transformation functions defined by the combinations of the 
parameters X and Y shown in table 5 are employed as the reference 
transformation functions, and the resolution of the original 
image signal to be processed is 200dpi will be described, 
hereinbelow. 
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Table 5 



10 lines /mm 


6.7 lines/imti 


5 lines /mm 


X 


Y 


X 


Y 


X 


Y 


15 


1.00 


20 


1.00 


20 


1.00 


20 


0.99 


30 


0.83 


30 


0.81 


30 


0.80 


40 


0.67 


40 


0.66 


40 


0.66 


50 


0.43 


50 


0.42 


50 


0.41 


60 


0.32 


60 


0.29 


60 


0.29 











Figure 31 shows the frequency response characteristics 
of band-limited signals obtained from an original image signal 
S^^g having a resolution of 200dpi. As can be seen from Figure 
31, the original image signal S^^^at 2 00dpi is about 4.0 cycles /mm 
(more strictly 3.937 cycles /mm) in Nyquist frequency, which 
falls on the peak frequency of the band-limited signal in the 
highest frequency band. The peak frequency of the band-limited 
signal in the second highest frequency band is 0.8 cycles/mm, 
1/5 of the Nyquist frequency. The peak frequency of the 
band-limited signal is reduced 1/2 by 1/2, e.g., 0.4 cycles/mm, 
0 . 2 cycles /mm/ 0 . 1 cycles /mm and 0 . 05 cycles/mm, as the frequency 
band of the band-limited signal becomes lower. This relation 
substantially conforms to the relation for the original image 
signal S^^^ at 10 lines/mm shown in Figure 24, and holds 
irrespective of the resolution of the original image signal S^^g 
so long as the filters used in making the low-resolution image 
signals are not changed. Figure 32 shows the frequency response 
characteristics of the processed image signal Sp^^^ obtained by 
processing the original image signal S^^^at 200dpi by use of the 
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reference transformation functions (broken line) and the 
processed image signal Sp^^^O obtained by processing the reference 
original image signal S^^^O at the reference resolution by use 
of the reference transformation functions (solid line) • As can 
5 be seen from Figure 32, the processed image signal S^^^^ is 
enhanced in frequency bands lower than the processed image signal 

Accordingly, as in the third embodiment, the 
transformation functions for transforming the band-limited 
,,,10 signals obtained from an original image signal S^^gat 2 00dpi are 
If obtained by correcting the reference transformation functions 

St on the basis of the resolution of the original image signal S^^g* 

^ In this embodiment, the peak frequencies of the band-limited 

"7 signals obtained from the reference original image signal S^^-g 

rl5 at 10 lines/mm, the parameters of the reference transformation 

functions (the values of Y in table 5) , and the peak frequencies 
of the band-limited signals obtained from the original image 
signal S^^g at 200dpi, and the number of the transformation 
functions to be obtained (six) are known. Accordingly, 
2 0 responses at peak frequencies of the band-limited signals 
obtained from the original image signal S^^g at 200dpi are first 
obtained according to the following formulae (9) on the basis 
of the frequency response characteristics of the band-limited 
signals obtained from the reference original image signal S^^g 
25 at 10 lines/mm, the parameters Y of the reference transformation 
functions and the peak frequencies of the band-limited signals 
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obtained from the original image signal S^^^ at 200dpi. 

R1=X1[F1]*Y1+X2 [Fl]*Y2+X3[Fl]*y3+X4 [Fl ] *Y4+X5 [Fl ] *Y5+X6 [Fl]*Y6 
R2=X1 [F2 ]*Y1+X2 [F2 ] *y2+x3 [F2 ] *y3+X4 [F2 ]*Y4+X5 [F2 ]*Y5+X6 [F2 ]*Y6 
R3=X1 [F3 ] *Y1+X2 [F3 ] *Y2+X3 [F3 ] *Y3+X4 [F3 ]*Y4+X5 [F3 ]*Y5+X6 [F3 ]*Y6 
5 R4=X 1 [ F4 ] *Y1+X2 [ F4 ] *Y2+X3 [ F4 ] *Y3+X4 [ F4 ] *Y4+X5 [ F4 ] *Y5+X6 [ F4 ] *Y6 

R5=X 1 [ F5 ] *Y1+X2 [ F5 ] *Y2+X3 [ F5 ] *Y3+X4 [ F5 ] *Y4+X5 [ F5 ] *Y5+X6 [ F5 ] *Y6 
R6=X1[F6]*Y1+X2 [F6]*Y2+X3 [F6 ] *Y3+X4 [F6 ]*Y4+X5 [F6 ]*Y5+X6 [F6 ]*Y6 

(9) 

wherein Rl to R6 represent responses at peak frequencies (200dpi) , 
0 XI [Fl ] to X6 [F6 ] represent responses of the band-limited signals 

obtained from the reference original image signal at frequencies 

of Fl to f6, Yl to Y6 represent the parameters of the reference 

transformation functions • 

Then by solving the following simultaneous equations 10, 
5 parameters Al to A6 of the transformation functions for the 

original image signal at 200dpi are obtained. 

Rl=Zl[Fl]*Al+Z2 [Fl]*A2 + Z3 [ Fl ] *A3 + Z4 [Fl ] *A4 + Z 5 [ Fl ] *A5 + Z 6 [Fl ] *A6 
R2=Z1[F2 ]*A1+Z2 [F2 ]*A2 + Z3 [F2 ]*A3 + Z4 [F2 ]*A4 + Z5 [F2 ]*A5 + 26 [F2 ]*A6 
R3=Z1[F3]*A1+Z2 [F3]*A2+Z3 [F3]*A3 + Z4 [F3 ] *A4 + Z 5 [F3 ] *A5 + Z 6 [F3 ] *A6 
0 R4 = Z1[F4]*A1+Z2 [F4]*A2 + Z3 [ F4 ] *A3 + Z4 [ F4 ] *A4 + Z 5 [ F4 ] *A5 + Z 6 [F4 ] *A6 

R5=Z1[F5]*AH-Z2 [F5 ]*A2+Z3 [F5 ]*A3+Z4 [F5 ]*A4+Z5 [F5 ]*A5+Z6 [F5 ]*A6 
R6=Z1 [F6 ]*A1+Z2 [F6 ] *A2+Z3 [F6 ] *A3+Z4 [ F6 ] *A4 + Z5 [ F6 ] *A5+Z 6 [F6 ] *A6 

(10) 

wherein Rl to R6 represent responses at peak frequencies (200dpi) 
5 obtained according to formulae (9), Zl[Fl] to Z6[F6] represent 

responses of the band-limited signals obtained from the original 
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image signal at 200dpi at frequencies of Fl to F6, Al to A6 
represent the parameters of the transformation functions for 
the original image signal at 200dpi. 

When the frequency enhancement processing is carried out 
5 on the original image signal at 200dpi by use of the 
transformation functions thus obtained, a processed image 
signal Sp^^^, which is substantially constant in frequency response 
characteristic can be obtained irrespective of the resolution 
of the original image signal, 
Glo A method of calculating Zl[Fl] to Z6[F6] in formulae (10) 

J: will be described, hereinbelow. As described above, the 

Of band-limited signals obtained from the reference original image 

signal S^^gO at 10 lines/mm have characteristics shown in Figure 
24. The characteristic of the band-limited signal in the 
fl5 highest frequency band shown by the solid line can be obtained 

4- by obtaining the frequency response characteristic of the 

unsharp image signal by multiplying Fourier-transformed value 
of the coefficient of filter of the one-dimensional filter shown 
in Figure 3A by Fourier-transformed value of the coefficient 
2 0 of filter of the one-dimensional filter shown in Figure 5, and 

subtracting the frequency response characteristic of the 
unsharp image signal thus obtained from the frequency response 
characteristic of the reference original image signal S^^gO (=1 
over the entire frequency bands ) . The characteristic of 
25 band-limited signals in lower frequency bands can be obtained 

by obtaining the frequency response characteristic of the 
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unsharp image signal in the similar manner on the basis of the 
coefficient of filter for obtaining the low-resolution image 
signal corresponding to the frequency band and the coefficient 
of filter for obtaining the unsharp image signal, and by 
subtracting the frequency response characteristic of the 
unsharp image signal thus obtained from the frequency response 
characteristic of the unsharp image signal in the frequency band 
higher by one stage. 

When the response of the band-limited signal in the highest 
frequency band of the reference original image signal S^^gO is 
represented by a function XI [Fi] (Fi being a frequency) and the 
values obtained by sampling the response of the band-limited 
signal in the highest frequency band of the reference original 
image signal at 2048 frequencies are represented by a function 
XI [ 1 ] , the relation between Fi and ± is as shown by the following 
formula ( 11 ) . 

i=2047xFi/fnq (11) 
wherein fnq represents a Nyquist frequency and 2047^1^0. This 
is shown in Figure 33. When the function Xl[i] is handled as 
a table, i represents addresses* By thus obtaining responses 
for the band-limited signals in all the frequency bands as 
functions XI [i] to X6[i] and handling the functions XI [i] to 
X6[i] as tables, values of Xl[Fl] to X6[F6] in the aforesaid 
formulae (9) and (10) can be easily obtained. 

The response Z1[F1] to Z6[F6] are obtained on the basis 
of the functions XI [i] to X6[i]. The response Zl[Fi] at a given 
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frequency Fi is obtained from the function XI [i] in the following 
manner. The response Zl[Fi] is the response of the band-limited 
signal in the highest frequency band in the band-limited signals 
obtained from the original image signal at 200dpi. The Nyquist 
frequency of the original image signal at 200dpi is 4 cycles/mm, 
and when the values obtained by sampling, at 2048 frequencies, 
the response of the band-limited signal at the frequency band 
at which the response is Zl[Fi] are represented by a function 
Zl[i] , the relation between Fi and i is as shown by the aforesaid 
formula (11) as in the case of XI [i]. Figure 34 shows the 
relation between the functions XI [i] and Zl[i]. Though, in 
concept, the functions XI [i] and Zl[i] are different functions 
which are different in Nyquist frequency, the response Zl[Fi] 
can be approximated by use of the function XI [i] . For example, 
the response Zl[Fi] at 2 cycles/mm is obtained in the following 
manner. That is, since 1=2047x2/4=1024 from the aforesaid 
formula (11), the response Zl[Fi] at 2 cycles/mm is given as 
the value of the function XI [i] at address of 1024. 

By handling the functions XI [i] to X6[i] as tables and 
obtaining the address i on the basis of formula (11) and the 
Nyquist frequency of the original image signal on which the 
frequency enhancement processing is to be carried out, and 
referring to the tables (the functions XI [i] to X6[i] ) according 
to the address i, the responses Z1[F1] to Z6[F6] can be 
approximated . 

When the frequency enhancement processing is to be carried 
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out on an original image signal S^^^ at 300dpi, the responses at 
frequencies corresponding to the peak frequencies of the 
band-limited signals obtained from the original image signal 
at 300dpi are obtained according to the aforesaid formulae (8) 
on the basis of the frequency response characteristics of the 
band-limited signals obtained from the reference original image 
signal, the peak frequencies of the band-limited signals 
obtained from the original image signal at 300dpi and the 
parameters Y of the reference transformation functions. Then 
the parameters for the transformation functions are obtained 
according to the aforesaid formulae ( 10 ) . In the original image 
signal at 3 00dpi, the Nyquist frequency is 5.9 cycles /mm and 
higher than the Nyquist frequency of the reference original image 
signal S^rgO • Since Xl[Fl] to X6[F6] in the aforesaid formulae 
(9) are functions for the range of 0 to 5 cycles /mm, no value 
is given for 5.9 cycles /mm by formulae (9) and Rl cannot be 
obtained. In this case, the peak frequencies of the band- 
limited signals obtained from the reference original image 
signal are employed as the peak frequencies of the band-limited 
signals obtained from the original image signal at 3 00dpi. 

The case where the frequency enhancement processing is 
to be carried out on a original image signal S^^^ at 100dpi will 
be described, hereinbelow. Since the Nyquist frequency of the 
original image signal at 100dpi is about 2.0 cycles /mm (strictly 
1*97 cycles/mm), the number of the band-limited signals is 5. 
Accordingly, the aforesaid formulae (9) and (10) are rewritten 
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as the following formulae (12) and (13) and five parameters for 
the transformation functions are obtained according to formulae 
(12) and (13). 

R1=X1[F1]*Y1+X2 [Fl]*Y2+X3[Fl]*y3+X4 [Fl ] *Y4+X5 [Fl ] *Y5+X6 [ Fl ] *Y6 
5 R2=X1[F2 ]*Y1+X2 [F2 ]*Y2+X3 [F2 ]*Y3+X4 [F2 ]*Y4+X5 [F2 ]*Y5+X6 [F2 ] *Y6 

R3=X1 [ F3 ] *Y1+X2 [ F3 ] *Y2+X3 [ F3 ] *Y3+X4 [ F3 ] *Y4+X5 [ F3 ] *Y5+X6 [ F3 ] *Y6 
R4=X1 [ F4 ] *Y1+X2 [F4 ] *Y2+X3 [F4 ] *Y3+X4 [F4 ] *Y4+X5 [ F4 ] *Y5+X6 [F4 ]*y6 
R5=X1[F5]*Y1+X2 [F5]*Y2+X3 [F5]*Y3+X4 [F5 ] *Y4+X5 [ F5 ]*Y5+X6[F5 ]*Y6 

(12) 

0 wherein Rl to R5 represent responses at peak frequencies (100dpi) , 
XI [Fl] toX6[F5] represent responses of the band-limited signals 
obtained from the reference original image signal at frequencies 
of Fl to F5, Yl to Y6 represent the parameters of the reference 
transformation functions . 

5 Rl=Zl[Fl]*Al + Z2 [Fl]*A2 + Z3 [F1]*A3+Z4 [Fl ] *A4 + Z5 [ Fl ] *A5 

R2=Z1[F2 ]*A1+Z2 [F2]*A2+Z3 [F2 ]*A3+Z4 [F2 ]*A4+Z5 IF2 ]*A5 
R3=Z1[F3]*A1+Z2 [F3]*A2+Z3 [F3]*A3+Z4 [F3 ]*A4+Z5 [F3 ] *A5 
R4=Z1 [F4 ]*A1+Z2 [F4 ] *A2+Z3 [F4 ] *A3+Z4 [F4 ]*A4+Z5 [F4 ]*A5 
R5=Zl[F5 ]*A1+Z2 [F5 ] *A2+Z3 [F5 ] *A3+Z4 [F5 ]*A4+Z5 [F5 ]*A5 

0 (13) 

wherein Rl to R5 represent responses at peak frequencies (100dpi) 
obtained according to formulae (12)^ Z1[F1] to Z5[F5] represent 
responses of the band-limited signals obtained from the original 
image signal at 100dpi at frequencies of Fl to F5, Al to A5 
5 represent the parameters of the transformation functions for 

the original image signal at 100dpi. 
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Though, in the third and fourth embodiments described 
above, the frequency enhancement processing represented by the 
aforesaid formulae (2) is carried out as a transformation 
processing, a dynamic range compression processing and a 

5 frequency enhancement processing may be executed simultaneously 
as a transformation processing. Figure 35 is a view showing a 
system for executing a frequency enhancement processing and a 
dynamic range compression processing as a transformation 
processing in an image processing system in accordance with a 

0 fifth embodiment of the present invention. As shown in Figure 

35, differences between two of the original image signal S^^^ and 
the unsharp image signals S^^k made by the filtering means 110 
and the interpolation means 111 are taken, whereby a plurality 
of band-limited signals (e-g-, S^^g-S^.l, S^^lS^^l and the like) 

5 respectively representing frequency components in limited 

frequency bands of the original image signal are made. The 
band-limited signals thus obtained are input into first and 
second transformation sections 103a and 103b and processed by 
respective transformation means 122a and 122b in the first and 

0 second transformation sections 103a and 103b. 

Transformation by the transformation means 122a in the 
first transformation section 103a is executed by use of the 
transformation functions calculated on the basis of the 
resolution of the original image signal S^^g as described above* 

5 For example, the transformation functions shown in Figure 36 

or 37 or combinations of the transformation functions shown in 
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Figures 36 and 37 are used as the reference transformation 
functions and transformation functions are calculated on the 
basis of the reference transformation functions according to 
the resolution of the original image signal S^^gto be processed 
5 so that the frequency response characteristic of the processed 

image signal S^^^^ becomes substantially constant irrespective 
of the resolution of the original image signal S^^^. 

In accordance with the transformation functions shown in 
Figure 36, the band-limited signals are transformed so that those 
|o which are large in amplitude are suppressed and the degree of 

^ suppression is increased as the frequency band of the band- 

f limited signal becomes higher. This is for taking into account 

^ the fact that higher frequency components contained in edges 

of an actual radiation image is smaller in amplitude than low 
15 frequency components. In actual radiation images, even a 

= ■ substantially sharp edge is not in the form of a correct step 

'i and the amplitude often becomes smaller as the frequency becomes 

higher. Accordingly, it is preferred that suppression be made 
from a smaller amplitude as the frequency of the band-limited 
20 signals becomes higher. The functions shown in Figure 36 serve 

for this purpose. 

The transformation functions shown in Figure 37 are for 
transforming the band-limited signals to those having values 
not larger than absolute values of the band-limited signals, 
25 which values are determined on the basis of the absolute values 
of the band-limited signals. As the frequency of the frequency 
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band to be processed by the function becomes lower, the absolute 
value of a transformed image signal obtained by transforming 
a band-limited signal whose absolute value is in a predetermined 
range near 0 becomes smaller. That is, the functions all pass 
through the origin and all have inclinations smaller than 1. 
Further, the inclination near 0 is smaller as the frequency of 
the frequency band to be processed by the function becomes lower. 
When an add signal obtained by adding up the transformed image 
signals is added to the original image signal S^^g, these 
functions contribute to smoothen the joint between the original 
image signal S^^^ and the add signal, that is, rise of the signals 

Similarly transformation by the transformation means 22b 
in the second transformation section 103b is executed by use 
of the transformation functions shown in Figure 3 8 or 37 or 
combinations of the transformation functions shown in Figures 
38 and 37. The transformed image signals output from the first 
and second transformation sections 103a and 103b are 
respectively input into operators 123a and 123b. The operator 
123a executes operation for making signals necessary for the 
frequency enhancement processing, and the operator 123b 
executes operation for making signals necessary for the dynamic 
range compression processing. 

The operator 123a executes a frequency enhancement 
processing similar to that in the third embodiment described 
above. That is, the band-limited signals transformed by the 
first transformation section 103a are added up, and an add signal 



thus obtained is multiplied by a coefficient of enhancement p 
which is determined according to the value of the original image 
signal S^^g. 

The operator 123b executes a dynamic range compression 
5 processing in the following manner. The band-limited signals 
transformed by the second transformation section 103b are added 
up^ and an add signal thus obtained is subtracted from the 
original image signal S^^g. Then the differential signal thus 
obtained is transformed on the basis of a transformation function^ 
0 thereby obtaining a coefficient of dynamic range compression. 

The signals respectively obtained by the operators 123a 
and 123b are added to the original image signal S^^g by an adder 
128^ whereby a processed image signal S^^^^ is obtained. 

The processing described above is represented by the 
5 following formula (3). 

+ D{S,,g-F^,,(S,,g, 8,3!, S^,2r S,,N)} (3) 

= { f ul ( Sorg-S.3l ) +f U2 ( ^^^1-8^^2 ) + 

0 + fu.(S..k-l-S,3k)+ - +f.N(SusN-l-S,,N) 

={ f di ( So,,-S„3l ) +f^, ( S„,l-S,,2 ) + - 

+ f,,(S„,k-l-S„3k)+ - +fd»(S„3N-l-S„,N) 
wherein S^^^^ is a processed image signal, S^rg is an original image 
5 signal, S^^k (k=l to N) is an unsharp image signal, f„k(k=l to 

N) is a transformation function used in the first transformation 
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section, fdk(k=l to W) is a transformation function used in the 
second transformation section, p{S^^g) is a coefficient of 
enhancement determined on the basis of the original image signal, 

and D{S,,,-Fd,,{S,,g, S.^l, S,,2, S^,N)} is a coefficient of dynamic 

5 range compression determined on the basis of the low frequency 

component signal, D being a function for transforming DiS^^^- 

In the image processing system shown in Figure 35, the 
add signal used in the frequency enhancement processing and that 

0 used in the dynamic range compression processing are separately 

obtained on the basis of the original image signal S^^g and are 
finally added. However, it is possible to first carry out on 
one of the frequency enhancement processing and the dynamic range 
compression processing on the original image signal S^^g and to 

5 carry out the other processing on the resultant signal. 

Generally the lower density portion of a radiation image includes 
a relatively large amount of noise since the lower density 
portion is exposed to less radiation upon taking the radiation 
image. Since the frequency enhancement processing is a 

0 density-dependent processing where the degree of enhancement 

is increased with increase in density, the low density portion 
is not enhanced and accordingly noise is not enhanced when the 
frequency enhancement processing is directly carried out on the 
original image signal. To the contrast, when the frequency 

5 enhancement processing is carried out after the dynamic range 

compression processing is carried out on the original image 
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signal, the low density portion is increased in its density by 
the dynamic range compression processing, and accordingly, the 
low density portion is enhanced by the following frequency 
enhancement processing, whereby noise is also enhanced. 
5 Accordingly, it is preferred that each signal be made on the 
basis of the original image signal S^^g. Further also from the 
viewpoint of saving the processing time, it is preferred that 
the two processings be executed in parallel* 

By calculating the transformation functions so that the 

CiO frequency response characteristic of the processed image signal 

Sp^oc becomes constant irrespective of resolution of the original 

nj image signal S^^^, the processed image signal Sp^oc can have a 

substantially constant frequency response characteristic 

- irrespective of the resolution of the original image signal S^^g. 

hts Though, in the third to fifth embodiments, the 

4- transformation functions are calculated by correcting the 

C:J reference transformation functions on the basis of the 

resolution of the original image signal S^^^ input, the 
transformation functions obtained may be further corrected on 

20 the basis of information on the response characteristics of the 
original image signal S^^g. That is, the original image signal 
So3,g is obtained from various apparatuses such as of a 
semiconductor sensor other than the radiation image read-out 
apparatus where a radiation image stored on a stimulable phosphor 

25 sheet is read out, and the response of the original image signal 
S^^g differs by apparatus to apparatus. Accordingly, by 
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correcting the transformation functions on the basis of 
information on the response characteristics of the original 
image signal S^^g, a processed image signal Sp^^^ having a constant 
frequency response characteristic can be obtained irrespective 
of the apparatus by which the original image signal S^^g is 
obtained. For example^ when the response at 2 cycles /mm of an 
original image signal S^^g obtained by a certain apparatus is 
inferior to a reference response by 40%, the fact is input into 
the image processing system, and the image processing system 
corrects the transformation functions so that response at 2 
cycles /mm becomes equal to the reference response. 

Further, though in the fourth and fifth embodiments 
described above, non-linear functions are used and a non-linear 
processing is carried out on the band-limited signals, linear 
functions or constants may be employed as the transformation 
functions . 

Further, though, in the third to fifth embodiments 
described above, the unsharp image signals are obtained from 
the original image signal by filtering and 
interpolation/enlargement and the band-limited signals are 
obtained from the original image signal and the unsharp image 
signals, the band-limited signals may be made, for instance, 
by transforming the original image signal to multiple resolution 
image signals by a wavelet transformation or a Laplacian pyramid, 
making the unsharp image signals from the image signals at the 
respective resolutions and making the band-limited signals from 
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the unsharp image signals. 

When the transformation functions for transforming 
band-limited signals obtained from an original image signal 
which is different from the reference original image signal in 
resolution are obtained by correcting the reference 
transformation functions, it is generally impossible to conform 
the frequency response characteristic of the processed image 
signal Sp,,^^ to the reference processed image signal S^^^^O (the 
processed image signal Sp^^^ obtained by transforming the 
band-limited signals from the reference image signal by use of 
the reference transformation functions) over the entire 
frequency bands. In an image obtained by reproducing a 
processed image signal, the low frequency components are 
visually more prominent than the high frequency components. 
Accordingly, it is preferred that the transformation functions 
for band-limited signals for the original image signal to be 
processed be calculated so that the frequency response 
characteristic of the processed image signal Sp^^^ conforms to 
that of the reference processed image signal Sp^^^O in frequency 
bands not lower than 1/5, more preferably not lower than 1/2, 
of the Nyquist frequency of the original image signal to be 
processed. For example, when the reference resolution is 5 
lines/mm and the resolution of the original image signal S^^gto 
be processed is 4 lines /mm, it is preferred that the 
transformation functions for transforming the band-limited 
signals obtained from the original image signal S^^g to be 
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processed be calculated on the basis of reference transformation 
functions for 5 lines /mm in the following manner. 

That is, since the image signal is reduced 1/2 by 1/2 when 
a plurality of band-limited signals are obtained, the band- 
5 limited signal in the highest frequency band is at least 1.0 
to 2.0 cycles /mm in characteristic width. Accordingly, in the 
image processing system in accordance with the third to fifth 
embodiments of the present invention, where the gains of the 
band-limited signals are controlled by the transformation 
QO functions, it is impossible to finely control the response in 

=C a frequency band between 1.0 and 2.0 cycles /mm and accordingly, 

IIJ when the frequency response characteristic sharply changes in 

€i the frequency band between 1.0 and 2.0 cycles /mm, it is 

impossible to obtain the transformation functions for the 4 
hl5 lines /mm which can precisely conform the frequency response 

£ characteristic of the processed image signal Sp^^^ to that of the 

0 reference processed image signal Sp^^^^O in the frequency band 

between 1.0 and 2.0 cycles /mm. Accordingly, it is preferred 
that the transformation functions for band-limited signals for 
20 the original image signal (at 4 lines/mm) to be processed be 

calculated so that the frequency response characteristic of the 
processed image signal Sp^^^ conforms to that of the reference 
processed image signal Sp^^^O in frequency bands not lower than 
1/5, more preferably not lower than 1/2, of the Nyquist frequency 
25 of the original image signal to be processed, as shown in Figures 

39A and 39B. In Figure 3 9A, the broken line shows the frequency 
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response characteristic of the reference processed image signal 
Sp^^^O obtained by transforming the band-limited signals from the 
reference image signal (having a picture element density of 5 
lines /mm) by use of reference transformation functions which 
5 are set so that relatively low frequency components are enhanced 
and the solid line shows the frequency response characteristic 
of the processed image signal Sp^^^ obtained by transforming the 
band-limited signals from the original image signal S^^g to be 
processed by use of the transformation functions which are 
ClO obtained by correcting the reference transformation functions 

^ according to the resolution (4 lines/mm) of the original image 

ni signal S^^g to be processed. In Figure 39B, the broken line shows 

^ the frequency response characteristic of the reference 

= processed image signal Sp^^^O obtained by transforming the 

Hi 5 band-limited signals from the reference image signal (having 

4- resolution of 5 lines/mm) by use of reference transformation 

O functions which are set so that relatively high frequency 

components are enhanced and the solid line shows the frequency 
response characteristic of the processed image signal Sp^^^ 
20 obtained by transforming the band-limited signals from the 

original image signal S^^g to be processed by use of the 
transformation functions which are obtained by correcting the 
reference transformation functions according to the resolution 
(4 lines/mm) of the original image signal S^^g to be processed. 
25 An image processing system in accordance with a sixth 

embodiment of the present invention will be described. 
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hereinbelow. The image processing system is for carrying out 
a frequency enhancement processing by use of unsharp image 
signals on an original image signal obtained by reading out a 
radiation image of a human body recorded on a stimulable phosphor 
sheet so that an image suitable for diagnosis can be obtained. 
An image reproduced on the basis of the processed image signal 
is mainly recorded on photographic film and used in diagnosis. 

In Figure 40, an image processing system 201 in accordance 
with a sixth embodiment of the present invention comprises an 
unsharp image signal making means 202 which makes a plurality 
of unsharp image signals from multiple resolution image signals 
S^^g' which are obtained from an original image signal S^^g by 
transforming the original image signal S^^g to multiple 
resolution space and encoding the transformed image signals, 
and a frequency enhancement processing means 203 which carries 
out a frequency enhancement processing for enhancing a 
particular frequency and obtains an intermediate processed 
image signal Sp,,^^'. The image processing system 201 further 
comprises a parameter setting means 204, a scale input means 
205 and an enlarging/contracting means 206. The parameter 
setting means 204 is a means which sets transformation functions 
which the frequency enhancement processing means 203 uses in 
the frequency enhancement processing on the basis of scale 
information S input from the scale input means 205. The scale 
input means 205 is a means for obtaining scale information S 
on the original image signal S^^g. The scale information S may 



be input as a value by the operator by a keyboard or by the operator 
selecting a value from a plurality of scales displayed on a 
control screen* The enlarging/contracting means 206 enlarges 
or contracts the intermediate processed image signal Sp^^c' 
5 obtained by the frequency enhancement processing means 203 on 
the basis of the scale information thereby making a final 
processed image signal Sp^^^. 

The original image signal S^^^is transformed into multiple 
resolution space and the transformed image signals are encoded 
C|o in the following manner- That is^ as shown in Figure 41A, the 

4- original image signal S^^g is subjected to wavelet transformation 

RJ and is decomposed into four pieces of data LLl, HLO, LHO and 

HHO by resolution. The data LLl represents an image obtained 
by reducing the original image to 1/2 in both longitudinal and 
Hb lateral directions. The data HLO, LHO and HHO represents a 

4" longitudinal edge, a lateral edge and an oblique edge. When the 

O four pieces of data LLl, HLO, LHO and HHO are subjected to inverse 

wavelet transformation, the original image signal S^^g is 
obtained. Then the data LLl is further subjected to wavelet 
20 transformation and decomposed into four pieces of data LL2, HLl, 
LHl and HHl as shown in Figure 4 IB. The data LL2 represents an 
image obtained by reducing data LLl to 1/2 in both longitudinal 
and lateral directions. The data HLl, LHl and HHl represents 
a longitudinal edge component, a lateral edge component and an 
25 oblique edge component of the data LLl. When the four pieces 
of data LL2, HLl, LHl and HHl are subjected to inverse wavelet 
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transformation, a low resolution image signal having a 
resolution equal to 1/2 of that of the original image signal 
is obtained- Further wavelet transformation is repeated a 
desired number of times on data LL obtained each time the wavelet 
5 transformation is carried out, thereby obtaining a plurality 
of pieces of data which are different in resolution . Thereafter , 
data at each resolution is encoded as shown in Figure 41C and 
the multiple resolution image signals S^^g' is obtained. 

By decoding only image signals up to a desired resolution 
Clo in the multiple resolution image signals S^^g' and carrying out 

4" inverse wavelet transformation on the image signals, a low 

[U resolution image signal representing an image at a resolution 

^3 of (k being the desired resolution) of that of the original 

image can be obtained. 
1^15 Making the unsharp image signals will be described, 

4^ hereinbelow. As shown in Figure 42, the unsharp image signal 

O making means 202 comprises first to n-th filtering means 210. 

The first filtering means 210 carries out a filtering processing 
on the original image signal S^^g, reconstructed from the multiple 
20 resolution image signals S^^g' , in x- and y-directions and makes 

a low resolution image signal Bl which is lower than the original 
image signal S^^g in resolution. The second filtering means 210 
carries out a similar filtering processing on the low resolution 
image signal Bl thus obtained and makes a low resolution image 
25 signal B2 which is lower than the low resolution image signal 

Bl in resolution, and the third filtering means 210 carries out 
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a similar filtering processing on the low resolution image signal 
B2 thus obtained and makes a low resolution image signal B3 which 
is lower than the low resolution image signal B3 in resolution. 
In this manner, the n-th filtering means 210 carries out a similar 
filtering processing on the low resolution image signal Bn- 
1 and makes a low resolution image signal Bn which is lower than 
the low resolution image signal Bn-1 in resolution. The unsharp 
image signal making means 202 further comprises first to n- 
th interpolation means 211. The interpolation means 211 carry 
out an interpolation/enlargement processing on the low 
resolution image signals Bl to Bn obtained at the respective 
stages of filtering, and makes a plurality of unsharp image 
signals S^.k (S„,l to S„,N) which are different in sharpness. 

In this particular embodiment, filters which 
substantially conform to one-dimensional Gaussian distribution 
are used in the filtering processing. That is, coefficients of 
the filters are determined according the following formula (4) 
which relates to a Gaussian signal. 

f(t) = e^" (4) 
This is because the Gaussian signal is good in localization in 
both a real space and a freguency space. For example, the 5x1 
one-dimensional filters are as shown in Figure 43 when a=l in 
formula ( 4 ) . 

The filtering is carried out on the original image signal 
S or the low resolution image signals every second picture 
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element as shown in Figure 44. When such filtering is carried 
out in both x- and y- directions, the number of picture elements 
in obtained low resolution image signal is reduced to 1/4 of 
the preceding image signal (e.g., the original image signal S^^g 
in the case of the low resolution image signal Bl, and the low 
resolution image signal Bl in the case of the low resolution 
image signal B2). That is, the number of picture elements in 
each of the low resolution image signals Bk(k stands for 1 to 
n) is 1/2^ of the original image signal S^rg- 

The interpolation/enlargement processing to be carried 
out on the low resolution image signals Bk thus obtained will 
be described, hereinbelow. Though various methods of 
interpolation such as B-spline can be employed, a Gaussian signal 
is also employed in the interpolation in this particular 
embodiment since low-pass filters based on Gaussian signals are 
employed in the filtering processing. Specifically, an 
approximation, a=2^"\ is employed in the following formula (5) . 

When interpolating the low resolution image signal Bl, 
a=l since k=l. In this case, a one-dimensional filter of 5x1 
such as shown in Figure 45 is employed in the interpolation. 
In this interpolation, a picture element of 0 in value is 
interpolated every second picture element in the low resolution 
image signal Bl, whereby the low resolution image signal Bl is 
enlarged to a size equal to the size of the original image, and 
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the interpolated low resolution image signal Bl is subjected 
to a filtering processing using the one-dimensional filter shown 
in Figure 45. 

This interpolation/enlargement processing is carried out 
on all the low resolution image signals Bk (Bl to Bn). When 
interpolating a low resolution image signal Bk, a filter which 
is 3x2^-1 in length is prepared according to formula (5) and 2^-1 
picture elements of 0 in value are interpolated between each 
pair of adjacent picture elements, whereby the low resolution 
image signal Bk is enlarged to a size equal to the size of the 
original image signal. Then the interpolated low resolution 
image signal Bk is subjected to a filtering processing using 

the filter which is 3x2^^-1 in length. 

The frequency enhancement processing to be carried out 
by use of the unsharp image signals S^^k thus obtained will be 
described, hereinbelow. Figure 46 shows a system for carrying 
out a frequency enhancement processing together with the unsharp 
image signal making means 202. As shown in Figure 46, unsharp 
image signals S^3k are made from the original image signal S^^g 
and differences between the original image signal S^^g and the 
unsharp image signals S.^k are taken by respective subtracters 
221, whereby band-limited signals (S^^g-S^^l, S^3l-S^32, and the 
like) which are components in limited frequency bands of the 
original image signal S^^g are made. 

The band-limited signals are transformed to predetermined 
amplitudes with different transformation functions f^ to f^;, in 



respective transformation circuits 222 and the transformed 
band-limited signals are added up and further added to the 
original image signal S„g by an operator 223 according to the 
following formulae (2'), whereby an intejrmediate processed 
5 image signal Sp^„„' in which a desired frequency components is 

enhanced to a required degree is generated. 

Fu,.(S„,,, S.^l, S„,2, S„,N) 

=f , ( S„,,-S,,l ) +f , ( S,3l-S„,2 ) + - 
eiO + f,,(S„3k-l-S,,k)+ - +f„(S,,N-l-S„,N) (2') 

4= wherein Sp„c' is an intermediate processed image signal in which 

nj the high-frequency components are enhanced, S^j-g is an original 

S image signal, S^^K (k=l to N) are unsharp image signals, fk(k=l 

to N) is a transformation function, and p(S„^,) is a coefficient 
ri5 of enhancement determined on the basis of the original image 

=1= signal. 

O The intermediate processed image signal Sp„„' is thus made. 

The problem to be solved by this embodiment and the manner of 
solving the problem in this embodiment will be described through 

20 an example, hereinbelow. As described above, though an 
intermediate processed image signal Sp^,,' representing an image 
equivalent to the original image signal in resolution (either 
may be the same image as the original image or an image different 
from the original image) can be obtained by carrying out 

25 frequency enhancement processing by use of transformation 
functions f^ determined according to the frequency bands of the 
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band-limited signals, an intermediate processed image signal 
Sp^oc' representing an image lower than the original image signal 
in resolution (either may be the same image as the original image 
or an image different from the original image) must be obtained 
5 by carrying out frequency enhancement processing by use of 
transformation functions fj, determined according to the 
resolution of the intermediate processed image signal Sp^^c to 
be obtained. When such transformation functions are stored for 
the respective resolutions, the number of transformation 
ClO functions to be stored becomes too large and management of the 

4^ transformation functions becomes to troublesome- When the 

nj frequency enhancement processing is carried out on an original 

image signal S^^gWith the original image signal S^^^ divided into 
six frequency bands, the frequency characteristics of the 
Hl5 band-limited signals have six peaks as shown in Figure 47. The 

JS frequencies of the respective peaks are 5 cycles/mm, 1 • 0 cycle/mm, 

O 0.5 cycles/mm, 0.25 cycles/mm, 0.12 cycles/mm and 0.06 cycles/mm. 

In an image signal representing an image whose resolution is 
1/2 of that of the original image, the peak frequencies of the 
20 band-limited signals are 2.5 cycles/mm, 0.5 cycles/mm, o.25 

cycles/mm, 0.125 cycles/mm and 0.06 cycles/mm as shown in Figure 
48. Further as shown in Figures 49 to 51, in image signals 
representing image whose resolutions are respectively 1/4, 1/8 
and 1/16 of that of the original image, the peak frequencies 
25 of the band-limited signals are equal to each other except their 

respective highest frequency bands. 
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In this embodiment, when a processed image signal 
representing an image having a desired resolution lower than 
the resolution of the original image is to be obtained, one of 
the low resolution signals representing a low resolution image 

5 closest to the desired resolution is employed as a reference 
low resolution image signal, low resolution band-limited 
signals are obtained on the basis of the reference low resolution 
image signal, and the low resolution band-limited signals are 
transformed into transformed image signals on the basis of 

0 transformation functions corresponding to the frequency bands 
of the band-limited signals, whereby an intermediate processed 
image signal S^^^^' is obtained. Then the intermediate processed 
image signal Sp^oc' is enlarged or contracted according to the 
scale information S and a processed image signal Sp^^^ 

5 representing an image having the desired resolution is obtained. 

Operation of the image processing system 201 of this 
embodiment when a processed image signal S^^^^ representing an 
image whose resolution is 1/2 of the original image signal is 
to be obtained will be described by way of example with reference 

0 to the flow chart shown in Figure 52 . It is assumed here that 

the band-limited signals obtained from the original image signal 
S^^g are in six frequency bands as shown in Figure 47. 

The user first inputs a desired scale through the scale 
input means 205. (step S21) Then the scale input means 205 

5 inputs a scale information S into the unsharp image signal making 

means 202, the parameter setting means 204 and the 
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enlarging/contracting means 206^ and the unsharp image signal 
making means 202 reconstructs, on the basis of the scale 
information S, a low resolution image signal representing an 
image having a resolution closest to the desired scale in the 

5 low resolution image signals Bj, ( step S22 ) , and makes low 

resolution unsharp image signals S^Jz. from the low resolution 
image signals (step S23) . In this example, since the resolution 
of the image represented by the processed image signal Sp^^c to 
be obtained (to be referred to as "the desired resolution", 

0 hereinbelow) is 1/2 of that of the image represented by the 

original image signal S^^g, unsharp image signals S^^l to S^^S are 
made. The correspondence between the unsharp image signals 
obtained from the original image signal S^^g (including the 
original image signal S^^^) and those obtained from a given low 

5 resolution image signal is as shown in Figure 53 . As can be seen 

from Figure 53, in the case where the desired resolution is 1/2 
of that of the original image, unsharp image signals S^gk are 
made with the unsharp image signal S^^l which is at the highest 
hierarchical level in the unsharp image signals obtained from 

0 the reference original image signal S^^gtaken as an original image 

signal S^^^^^^- Accordingly, in this case, the low resolution 
original image signal corresponds to the unsharp image 

signal S^^l obtained from the reference original image signal 
S^^g, the unsharp image signals S^^li/s corresponds to the unsharp 

5 image signals S^^2 obtained from the reference original image 

signal S^^g, the unsharp image signals S^^2^/2 corresponds to the 
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unsharp image signals S^J obtained from the reference original 
image signal S^^g, and so on. Then low resolution band-limited 
signals are made on the basis of the original image signal S^^^^^^ 
and the low resolution unsharp image signals S^gli/s to S^^5^^2' 
(step S24 ) The correspondence between the band-limited signals 
obtained from the reference original image signal S^^g and those 
obtained from a given low resolution image signal is as shown 
in Figure 54 . Since even if the band-limited signals not higher 
than 0.06 cycles/mm are processed, the processed image signal 
Sp^^, is not substantially improved in quality, the low frequency 
band-limited signals not higher than 0.06 cycles/mm are not 
processed. 

The parameter setting means 204 sets the transformation 
functions fj, for transforming the low resolution band-limited 
signals. (step S25) Since the peak frequencies of the low 
resolution band-limited signals are equal to the band-limited 
signals obtained from the reference original image signal S^^g 
except their respective highest frequency bands as shown in 
Figures 47 and 48, transformation functions fj to f^ are employed. 
The correspondence between the transformation functions for the 
band-limited signals obtained from the reference original image 
signal S^^g and those for the band-limited signals obtained from 
a given low resolution image signal is as shown in Figure 55. 
When the desired resolution is 1/2'' of the resolution of the 
reference original image signal, the transformation functions 
on the line lower than the transformation functions for the 
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reference original image signal S^^g by n lines in Figure 55 are 
employed. The following formula (14) represents Fusm in the 
aforesaid formulae (2') in carrying out frequency enhancement 
processing on the reference original image signal S^^g and the 
following formula ( 15 ) represents Fusm in the aforesaid formulae 
(2') in carrying out frequency enhancement processing on the 
original image signal S^^^^j^ at a resolution of 1/2 of that of 
the reference original image signal S^^g. 

F.sn.(Soxgf S,,l, S^^r S,,6) 

=f 1 ( S,,g-S,3l ) +f 2 ( S,3l-S,,2 ) + - 

+ f 3 (S,34-S,35)+f,( 3,35-5,36) (14) 

Fusm( Sorgl/2 / ^us^l/2f S,s2i/2/ ^us^l/l) 
= f 1 ( Sorgl/2~^usi 1/2 ) 2 ( ^usl l/2"S,s2 ) 

+ f,( 3,33-3,34 )+f 5 (S,.4-S,35) (15) 

Then low resolution transformed image signals are made 
from the low resolution band-limited signals (step 326), and 
a low resolution add signal is obtained- Then the low resolution 
add signal is multiplied by the coefficient of enhancement 
P(Sorgi/2) arid the product is added to the low resolution original 
image signal whereby a low resolution intermediate 

processed image signal Sp^^^i/2' is obtained. The coefficient of 
enhancement ^{S,,^^^^) is set according to the low resolution 
original image signal 3^^gi/2* 

The low resolution intermediate processed image signal 
3 wo' thus obtained is input into the enlarging/contracting 
means 206, and the enlarging/contracting means 206 enlarges or 



contracts the low resolution intermediate processed image 
signal Sp^,,!// on the basis of the scale information S input from 
the scale input means 205 so that a low resolution image can 
be reproduced at the desired scale, thereby obtaining a final 
5 low resolution processed image signal Sp^^^i/j. In this particular 

example, since the desired resolution is 1/2, the 
enlarging/contracting means 206 need not enlarge or contract 
the low resolution intermediate processed image signal Sp^^^i/a' . 
As can be understood from the description above, since, 
Cio in this embodiment, the low resolution band-limited signals 
4= obtained from the low resolution original image signal are 

OJ transformed on the basis of a part of the transformation 

^1 functions for transforming the band-limited signals obtained 

from the reference original image signal which correspond to 
Hb the low resolution band-limited signals in frequency bands, the 

£ image reproduced on the basis of the low resolution processed 

O image signal Sp^^.^/j substantially conforms in frequency response 

characteristic to the image reproduced on the basis of the 
processed image signal Sp^^^ obtained from the reference original 
20 image signal, whereby a processed image signal which can 

reproduce an image which is substantially constant in frequency 
response characteristic irrespective of resolution can be 
obtained. 

Further since a part of the transformation functions which 
25 are employed for frequency enhancement processing of the 
reference original image signal are employed, transformation 
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functions need not be prepared for each resolution, which 
simplifies the arrangement of the image processing system and 
saves trouble to manage a large number of transformation 
functions . 

5 When an image signal at a desired resolution is obtained 

by contracting a processed image signal 3^^^^ obtained from the 
reference original image signal S^^g, making unsharp image 
signals, making band-limited signals and transformation of the 
band-limited signals must be carried out using all the low 
lo resolution image signals. To the contrast, in this embodiment, 

^ since only the low resolution image signals up to that 

f representing an image having a resolution closest to the desired 

? resolution have to be reconstructed, operation time required 

for the processing can be shortened. 
l5 Further, though, in the sixth embodiment described above, 

^ the unsharp image signals are obtained from the original image 

^ signal (the reference original image signal or the low resolution 

original image signal) by filtering and 

interpolation/enlargement and the band-limited signals are 
20 obtained from the original image signal and the unsharp image 
signals, the band-limited signals may be made, for instance, 
by transforming the original image signal to multiple resolution 
image signals by a wavelet transformation or a Laplacian pyramid, 
making the unsharp image signals from the image signals at the 
25 respective resolutions and making the band-limited signals from 
the unsharp image signals. 
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Though, in the sixth embodiment described above, frequency 
enhancement processing is carried out on the encoded multiple 
resolution image signals S,^/ , the frequency enhancement 
processing may be carried out on an original image signal S^^g 
5 which is not encoded. 

An image processing system in accordance with a seventh 
embodiment of the present invention will be described, 
hereinbelow. In the sixth embodiment, low resolution image 
signals up to the resolution closest to a desired resolution 
10 are reconstructed from an encoded multiple resolution image 

^ signals S^^g' and unsharp image signals are made on the basis of 

the reconstructed low resolution image signals . To the contrast, 
i in the seventh embodiment, unsharp image signals are made on 

the basis of low resolution image signals from the lowest 
15 resolution to a resolution closest to the resolution 

- corresponding to a desired scale. The seventh embodiment 

i differs from the sixth embodiment only in the processing executed 

by the unsharp image signal making means, and accordingly the 
processing executed by the unsharp image signal making means 
20 in this embodiment will be mainly described, hereinbelow. 

Figure 56 shows a unsharp image signal making means 202' 
in the image processing system of the seventh embodiment. As 
shown in Figure 56, the unsharp image signal making means 202' 
comprises a plurality of encoders 230 and a plurality of 
25 interpolation means 231 . That is, a low resolution image signal 

representing an image at a resolution of 1/2" of that of the 
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original image represented by the reference original image 
signal S^^^ is first made by decoding an image signal representing 
an image which the lowest in resolution in the images represented 
by the multiple resolution image signals S^^g' and carrying out 
an inverse wavelet transformation on the decoded image signals 
Then a low resolution image signal B„_i representing an image at 
a resolution of 1/2""^ of that of the original image is made on 
the basis of the low resolution image signal B^. In this manner, 
low resolution image signals Bj, whose resolution are 1/2^ (k=l 
to n) of that of the original image signal S^^g are made as well 
as the original image signal S^^g- The low resolution signals 
Bj, correspond to data LLk obtained by carrying out a wavelet 
transformation on the original image signal S^^g- Then the 
interpolation means 231 carries out an 

interpolation/enlargement processing on each of the low 
resolution image signal B^ and a plurality of unsharp image 
signals S^^k (k=l to N) which are different in sharpness are made. 

Then band-limited signals are obtained from the unsharp 
image signals S^3k and an intermediate processed image signal 
Sproc' is obtained according to the aforesaid formulae (2'). 

Operation of the image processing system of this 
embodiment when a processed image signal S^^^^ representing an 
image whose resolution is 1/2 of the original image signal is 
to be obtained will be described by way of example. It is assumed 
here that the band-limited signals obtained from the original 
image signal S^^g are in six frequency bands. 



The user first inputs a desired scale through the scale 
input means 2 05, Then the scale input means 205 inputs a scale 
information S into the unsharp image signal making means 202', 
the parameter setting means 204 and the enlarging/contracting 
5 means 206, and the unsharp image signal making means 202' 
reconstructs, on the basis of the scale information S, low 
resolution image signals Bj, from the lowest resolution to a 
resolution closest to the resolution corresponding to the 
desired scale and makes low resolution unsharp image signals 
131 0 S^gk. In this example, since the resolution of the image 

=p represented by the processed image signal Sp^^^ to be obtained 

HJ (to be referred to as "the desired resolution", hereinbelow) 

m is 1/2 of that of the image represented by the original image 

^ signal S^^^, unsharp image signals S^^l to S^^5 are made. The 

hi 5 correspondence between the unsharp image signals obtained from 

4^ the original image signal S^^g (including the original image 

O signal S^^g) and those obtained from the low resolution image 

signals Bj, from the lowest resolution to a resolution closest 
to the resolution corresponding to the desired scale is as shown 
20 in Figure 53. The correspondence between the band-limited 

signals obtained from the reference original image signal S^^g 
and those obtained from the low resolution image signals Bj, from 
the lowest resolution to a resolution closest to the resolution 
corresponding to the desired scale is as shown in Figure 54. 
25 The parameter setting means 204 sets the transformation 

functions fj, for transforming the low resolution band-limited 
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signals in the same manner as in the sixth embodiment* When the 
desired resolution is 1/2" of the resolution of the reference 
original image signal, the transformation functions on the line 
lower than the transformation functions for the reference 
original image signal S^^^ by n lines in Figure 55 are employed. 

Then low resolution transformed image signals are made 
from the low resolution band-limited signals, and a low 
resolution add signal is obtained. Then the low resolution add 
signal is multiplied by the coefficient of enhancement p(S^^gi/2) 
and the product is added to the low resolution original image 
signal S^^gi/2f whereby a low resolution intermediate processed 
image signal Sp^-^^i/a' obtained. 

The low resolution intermediate processed image signal 
Sproci/2' thus obtained is input into the enlarging/contracting 
means 206, and the enlarging/contracting means 206 enlarges or 
contracts the low resolution intermediate processed image 
signal S^^^^i// on the basis of the scale information S input from 
the scale input means 2 05 so that a low resolution image can 
be reproduced at the desired scale, thereby obtaining a final 
low resolution processed image signal m this particular 

example, since the desired resolution is 1/2, the 
enlarging/contracting means 206 need not enlarge up or contract 
the low resolution intermediate processed image signal Sp^^^i/s' . 

As can be understood from the description above, since, 
also in this embodiment, the low resolution band-limited signals 
obtained from the low resolution original image signal are 
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transformed on the basis of a part of the transformation 
functions for transforming the band-limited signals obtained 
from the reference original image signal which correspond to 
the low resolution band-limited signals in frequency bands, the 
5 image reproduced on the basis of the low resolution processed 
image signal Sp^^^,/2 substantially conforms in frequency response 
characteristic to the image reproduced on the basis of the 
processed image signal Sp^^^ obtained from the reference original 
image signal, whereby a processed image signal which can 
C|0 reproduce an image which is substantially constant in frequency 

± response characteristic irrespective of resolution can be 

^Jf obtained. 

^•1 Though, the above description has been made by way of 

i, example on the case where the desired resolution is 1/2 of that 

rl5 of the reference original image signal, in the case where the 

± desired resolution is, for instance, 1/3 of that of the reference 

K original image signal, the unsharp image signal making means 

202 and the frequency enhancement processing means 203 carry 
out the frequency enhancement processing on the low resolution 
20 image signal representing an image whose resolution (1/4 of 

the original image signal) is closest to 1/3, thereby making 
an intermediate processed image signal Sp^^^i// . Then a final 
processed image signal Sp^^^^^^ is obtained by enlarging the 
intermediate processed image signal Sp^^^^/^' 4/3 times by the 
25 enlarging/contracting means 206. When the desired resolution 

is 1/2'' of that of the original image, the enlarging/contracting 
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means 206 may be eliminated, whereby the system can be 
simplified. 

Though, in the sixth and seventh embodiments, the 
transformed image signals are made by transforming the 
band-limited signals on the' basis of the transformation 
functions f^ according to the aforesaid formulae (2'), the 
transformed image signals may be obtained by enhancing the 
band-limited signals by use of coefficients of enhancement a^, 
(k stands for 1 to N) in place of the transformation functions 
fj, as represented by the following formulae (16). 

Fusm(Sorg. S,,l , 8,32, S.^N ) 

=cx, ( S^^g-S^.l ) ( 3,31-8,32 ) + - 

+ aj,(S,,k-l-S,,k)+ - +a^(S,3N-l-S,,N) (16) 
When an original image signal to be processed is higher 
than the reference original image signal S^^g in resolution, the 
number of the band-limited signals and the number of the 
transformation functions for transforming the band-limited 
signals are larger than those for the reference original image 
signal S^^,^. Figure 57 shows the correspondence between the 
band-limited signals obtained from the reference original image 
signal S^^g and those obtained from original image signals which 
are twice and four times the reference original image signal 
S^^g in resolution. When the object original image signal to be 
processed has a resolution twice that of the reference original 
image signal S^^g (this object original image signal is denoted 
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by Sorg2 in Figure 57, and the unsharp image signals obtained from 
the object original image signal S^^^^ are denoted by S^.^k ( k stands 
for an integer larger than 0) ) , there is no band-limited signal, 
in the band-limited signals obtained from the reference original 
5 image signal S^^,, corresponding to band-limited signal S^^^^- 

S^^jl in the highest frequency band in the band-limited signals 
obtained from the object original image signal S<,rg2f and 
accordingly there is no transformation function for 
transforming the band-limited signal S^^g^-S^.^l in the reference 
10 transformation functions for transforming the band-limited 

i signals obtained from the reference original image signal S^^,. 

J Similarly when the object original image signal to be processed 

i has a resolution four times that of the reference original image 

signal S„^g (this object original image signal is denoted by S^^^^ 
U5 in Figure 57, and the unsharp image signals obtained from the 

object original image signal S^,^ are denoted by S^.^k (k stands 
for an integer larger than 0) ) , there is no band-limited signal, 
in the band-limited signals obtained from the reference original 
image signal S^j-g/ corresponding to band-limited signals S^^^^-S^^il 
20 and S„,4l- S,,42 in the highest and second highest frequency bands 
in the band-limited signals obtained from the object original 
image signal S^^g^, and accordingly there is no transformation 
function for transforming the band-limited signals S^^^i-S^^^l and 
Sus4l~ Sus42 in the reference transformation functions for 
25 transforming the band-limited signals obtained from the 

reference original image signal S^^^. In such cases, the 
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transformation function f 1 for the band-limited signal 
in the highest frequency band in the band-limited signals 
obtained from the reference original image signal S^^g is used 
as the transformation functions for transforming the band- 
limited signal(s) in the highest and second highest frequency 
bands obtained from the object original image signal as shown 
in Figure 58, 

When there are prepared a plurality of reference original 
image signals S,,,, it is preferred that one of the reference 
original image signal 3,,^ closest to the object original image 
signal in resolution be taken as the reference original image 
signal S^^g in the manner described above in conjunction with the 
first and second embodiments. It is preferred that the 
transformation functions which are determined for the band- 
limited signals obtained from the object original image signal 
be related to the object original image signal and stored 
together with the resolution of the reference original image 
signal. With this arrangement, the image represented by the 
object original image signal can be constantly reproduced with 
the same characteristic. 

An image processing system in accordance with a seventh 
embodiment of the present invention will be described, 
hereinbelow. The image processing system is for carrying out 
a frequency enhancement processing by use of unsharp image 
signals on an original image signal obtained by reading out a 
radiation image of a human body recorded on a stimulable phosphor 



sheet so that an image suitable for diagnosis can be obtained. 
An image reproduced on the basis of the processed image signal 
is mainly recorded on photographic film and used in diagnosis. 

In Figure 59, an image processing system 301 in accordance 
with the seventh embodiment of the present invention comprises 
a picture element density transformation means 304, an unsharp 
image signal making means 302, a frequency enhancement 
processing means 303, an enlarging/contracting means 305 and 
a parameter correcting means 306. When frequency enhancement 
processing is to be carried out on an original image signal S^^g 
itself, the picture element density transformation means 304 
passes the original image signal S^^g as it is to the unsharp image 
signal making means 302 which makes a plurality of unsharp image 
signals from the original image signal S^^g. The unsharp image 
signal making means 302 inputs the unsharp image signals into 
the frequency enhancement processing means 303. The frequency 
enhancement processing means 303 carries out the frequency 
enhancement processing for enhancing a particular frequency to 
obtain an intermediate processed image signal S^^^^' and inputs 
the intermediate processed image signal Sp^^^' into the 
enlarging/contracting means 305. In this case, the 
enlarging/contracting means 305 outputs the intermediate 
processed image signal Sp^,,' as it is as a final processed image 
signal Sp^^^. To the contrast, when the frequency enhancement 
processing is to be carried out on an original image signal S^^g 
which is obtained by picture element density transformation of 



a base original image signal S^^/ which is read out by an image 
read-out apparatus or the like and has a reference predetermined 
picture element density (10 lines /mm in this particular 
embodiment), the picture element density transformation means 
5 304 transforms the picture element density of the base original 
image signal S^^/ to a desired picture element density, e.g., 
6.7 lines /mm or 5 lines/mm, thereby obtaining an original image 
signal S^^g at the desired picture element density. Then the 
picture element density transformation means 3 04 inputs the 
ao picture element density-transformed original image signal S^^g 
£ into the unsharp image signal making means 302. The unsharp 

nJ image signal making means 3 02 makes a plurality of unsharp image 

ul signals from the original image signal S^^g and inputs the unsharp 

^ image signals into the frequency enhancement processing means 

1-15 303 . The frequency enhancement processing means 303 carries out 

i~" the frequency enhancement processing for enhancing a particular 

P frequency to obtain an intermediate processed image signal Sp^^^' 

and inputs the intermediate processed image signal Sp^^^' into 
the enlarging/contracting means 3 05. In this case, the 
20 enlarging/contracting means 305 enlarges the intermediate 
processed image signal S^^,,' by interpolation at a desired 
magnification and outputs the enlarged intermediate processed 
image signal Sp^^^' as a final processed image signal Sp^,^^. In 
this specification, the term ^^magnification" should be broadly 
25 interpreted to include a magnification smaller than 1. The 

desired picture element density may be input through a keyboard 
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or the like as a value or may be selected from a plurality of 
picture element densities displayed on a control screen. 
Similarly the desired magnification may be input through a 
keyboard or the like as a value or may be selected from a plurality 
of magnifications displayed on a control screen. 

The picture element density transformation processing 
will be described, hereinbelow. In the aforesaid radiation 
image read-out system using a stimulable phosphor sheet, the 
read density or the picture element density differs depending 
on the size of the stimulable phosphor sheet and may be freely 
changed as required by the user. In this embodiment, a base 
original image signal S,,/ read at 10 lines/mm is transformed 
to an original image signal 8,,^ at 6.7 lines/mm or 5 lines/mm 
by carrying out a filtering processing on the base original image 
signal S^^g' the picture element density transformation means 
304. In this specification, the original image signal S^^g 
obtained by carrying out the picture element density 
transformation on the base original image signal S^^g' will be 
referred to as "the picture element density- trans formed 
original image signal S,^/' , hereinbelow, in order to clearly 
distinguish the former from the latter. 

The unsharp image signal making means 302 makes unsharp 
image signals S^.k in the same manner as described above and the 
frequency enhancement processing means 303 carries out the 
frequency enhancement processing by use of the unsharp image 
signals S^.k in the same manner as described, whereby the 
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intermediate processed image signal Sp^^^' is obtained. 

The enlarging/contracting means 305 enlarges the 
intermediate processed image signal Sp^^^' by interpolation using 
an interpolation filter at a desired magnification, thereby 
5 obtaining a final processed image signal Sp„c' 

Desired frequency components of the original image signal 
S„rg can be enhanced by thus carrying out on the original image 
signal S^^g a frequency enhancement processing by use of 
transformation functions f„ which are set for the respective 
10 frequency bands of the band-limited signals. However, when the 

P frequency enhancement processing is carried out on the picture 

IJ element density-transformed original image signal S„gby use of 

0 the same transformation functions f^ as those used for carrying 

out the frequency enhancement processing on the base original 
15 image signal S„g' > the image reproduced on the basis of the 

P processed image signal Sp^„„ obtained from the picture element 

J density- trans formed original image signal S„^g becomes 

different in frequency response characteristic from that 
reproduced on the basis of the processed image signal Sp„^ 
obtained from the base original image signal S,,^,' and the former 
image becomes less sharp than the latter image though the images 
are the same in size. This is due to the picture element density 
transformation which is carried out when the picture element 
density- trans formed original image signal S^^, is obtained and 
the interpolation processing carried out to enlarge the 
intermediate processed image signal Sp„„'. 
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Accordingly, in this embodiment, the transformation 
functions for transforming the band-limited signals obtained 
from the picture element density-transformed original image 
signal S„^g are obtained by correcting the parameters of the 
transformation functions f^ for transforming the band-limited 
signals obtained from the base original image signal S„^g' 
according to the filtering characteristic Fl of the picture 
element density transformation filter which is employed to carry 
out the picture element density transformation processing by 
the picture element density transformation means 304 and the 
filtering characteristic F2 of an interpolation filter which 
is employed to carry out the interpolation by the 
enlarging/contracting means 305 so that the frequency response 
characteristic of the image reproduced on the basis of the 
processed image signal Sp^„^ obtained from the picture element 
density- trans formed original image signal S^^g becomes equal to 
that reproduced on the basis of the processed image signal Sp^„^ 
obtained from the base original image signal S^^^' . 

Correction of the parameters of the transformation 
functions by the parameter correcting means 3 06 will be described, 
hereinbelow. In the description below. 

Figure 60 shows the filtering characteristic Fl of the 
picture element density transformation filter which is employed 
by the picture element density transformation means 304 to 
transform the base original image signal S^^,' having a picture 
element density of 10 lines /mm to an picture element 
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density- trans formed original image signal S^^g having a picture 
element density of 5 lines/mm, and Figure 61 shows the filtering 
characteristic F2 of the interpolation filter which is employed 
by the enlarging/contracting means 305 to double the picture 
5 element density-transformed original image signal S^^g. The 

picture element density-transformed original image signal S^^g 
obtained by transforming, by use of the picture element density 
transformation filter whose filtering characteristic is Fl, the 
base original image signal S^^^g' having a picture element density 
fflo of 10 lines /mm to an picture element density-transformed 
3^ original image signal S^^g having a picture element density of 

Jjf 5 lines /mm will be referred to as "the 5-line transformed 

original image signal S^^^5" hereinbelow. Figure 62 shows the 
frequency response characteristic of the processed image signal 
ft 5 Sp^^^ obtained by carrying out the frequency enhancement 

processing by use of the transformation functions f^ set for the 
frequency enhancement processing of the base original image 
signal S^^g' (these transformation functions will be referred to 
as "the base transformation functions fy,^" hereinbelow) on the 
20 5-line transformed original image signal S^^gS and by doubling 

the obtained intermediate processed image signal S^^^^' by use 
of the interpolation filter whose filtering characteristic is 
F2 . The processed image signal Sp^^^^thus obtained from the 5-line 
transformed original image signal S^^gS will be referred to as 
25 "the 5-line processed image signal S^^^^S" hereinbelow. Further 
the processed image signal Sp^^^ obtained by frequency-enhancing 
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the base original image signal S^^^g' by use of the base 
transformation functions fj^^will be referred to as the "the base 
processed image signal S^^^^B" , hereinbelow. In Figure 62, the 
frequency response characteristic of the 5-line processed image 
signal Sp^^^S is shown with the frequency response characteristic 
of the base processed image signal S^^^^B taken as 1 • As can be 
seen from Figure 62, the frequency response characteristic of 
the 5-line processed image signal S^^^^5 can be easily calculated 
on the basis of the filtering characteristic Fl of the picture 
element density-transformation filter and the filtering 
characteristic F2 of the interpolation filter. Further, the 
difference in frequency response characteristic between the 
base processed image signal Sp^^^B and the 5-line processed image 
signal S^^^^S can be easily obtained from Figure 62 • 

For example, in order to equalize aspect of the image 
reproduced on the basis of the 5-line processed image signal 
Sproc^ to that of the image reproduced on the basis of the base 
processed image signal Sp^^^^B at 2 cycles /mm, it is necessary to 
double the frequency response at 2 cycles /mm of the 5-line 
processed image signal Sp^^^S. Accordingly, in this embodiment, 
the transformation function for the band-limited signal in the 
frequency band corresponding to 2 cycles /mm is set by correcting 
the parameters for the base transformation function fj,j^ for the 
band-limited signal in the frequency band corresponding to 2 
cycles /mm so that the degree of enhancement is doubled. As shown 
in Figure 63, the band-limited signal in the frequency band 



corresponding to 2 cycles /mm is in the highest frequency band, 
and accordingly the parameter of the transformation function 
fi for the band-limited signal in the highest frequency band is 
corrected so that the degree of enhancement is doubled. More 
specifically, when the transformation function is a function 
having an inclination, the inclination is doubled, and when the 
transformation function is a constant, the value of the constant 
is doubled. The frequency component at which aspect of the image 
reproduced on the basis of the 5-line processed image signal 
SprocS is to be equalized to that of the image reproduced on the 
basis of the base processed image signal Sp^^^B may be any 
frequency and may be two or more frequencies. Further the 
frequency component may be input by the user through an input 
means not shown* 

Operation of the image processing of this embodiment will 
be described with reference to the flow chart shown in Figure 
64, hereinbelow. It is assumed that a 5-line processed image 
signal Sp^^^^S is obtained by transforming a base original image 
signal S^^^' having a picture element density of 10 lines/mm to 
a 5-line transformed original image signal S^^^S having a picture 
element density of 5 lines /mm, carrying out the frequency 
enhancement processing on the 5-line transformed original image 
signal S^^^S , thereby obtaining an intermediate processed image 
signal S^^^^', and doubling the intermediate processed image 
signal Sp,,^^' . The base original image signal S^^g' is first input 
into the image processing system 301. (step S31) Then the 



picture element density transformation means 304 transforms the 
base original image signal S^^^' to the 5-line transformed 
original image signal S^^g5. (step S3 2) The 5-line transformed 
original image signal S^^g5 is input into the unsharp image signal 
5 making means 302 and the unsharp image signal making means 302 
makes unsharp image signals S^^k from the 5-line transformed 
original image signal S^^^gS. (step S33) The filtering 
characteristics Fl and F2 are input into the parameter correcting 
means 306. Further information on the frequency component at 
^4o which aspect of the image reproduced on the basis of the 5- 
J" line processed image signal S^^^^5 is to be equalized to that of 

vt the image reproduced on the basis of the base processed image 

signal S^^^^B is input into the parameter correcting means 306 
(the frequency component will be referred to as "the target 
ri5 frequency component", hereinbelow) . (step 534) The parameter 

2! correcting means 306 corrects the parameters of the base 

transformation functions f^^, on the basis of the filtering 
characteristics Fl and F2 and the information on the target 
frequency component. (step S3 5 ) Steps S4 and S5 may be executed 
20 before steps S2 and S3 or in parallel thereto. 

Then the frequency enhancement processing means 3 03 makes 
band-limited signals on the basis of the unsharp image signals 
Sygk and carries out the frequency enhancement processing on the 
unsharp image signals S^^k by use of the transformation functions 
25 defined by the corrected parameters, thereby obtaining the 
intermediate processed image signal S^^^^ ' . ( step 83 6 ) Then the 
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enlarging/contracting means 305 enlarges the intermediate 
processed image signal Sp^.^^', thereby obtaining the final 5- 
line processed image signal Sp^^^S. (step S37) Finally, an image 
is reproduced by a printer or the like on the basis of the 5-line 

5 processed image signal Sp^o^5. (step S38) 

Thus, in this embodiment, the parameters of the 
transformation functions for carrying out the frequency 
enhancement processing on the picture element density- 
transformed original image signal S^rg are determined by 

0 correcting the parameters of the base transformation functions 
ffcb (the transformation functions for carrying out the frequency 
enhancement processing on the base original image signal S^^g') 
on the basis of the filtering characteristics Fl and F2 so that 
the frequency response characteristic of the processed image 

5 signal Sp^,^^ becomes substantially equal to that of the base 

processed image signal Sp^^^B, and accordingly the image 
reproduced on the basis of the processed image signal S^^^^ 
obtained from the picture element density- trans formed original 
image signal S^^g can be equivalent in sharpness to the image 

0 reproduced on the basis of the base processed image signal Sp^^^B 

obtained from the base original image signal S^^^' . Further since 
the frequency response characteristic is corrected 
simultaneously with the frequency enhancement processing, the 
operation time can be saved as compared with when correction 

5 of the frequency response characteristic is carried out on the 

processed image signal Sp^^^ which has been subjected to the 
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frequency enhancement processing. 

Though, in the seventh embodiment, the parameters of the 
base transformation functions f^t are corrected according to the 
filtering characteristic Fl of the picture element density- 

5 transformation filer and the filtering characteristic F2 of the 
interpolation filter, it is possible to prepare in advance a 
table in which the kind of the picture element density 
transformation to be carried out by the picture element density 
transformation means 3 04 and the kind of the interpolation filter 

0 to be used by the enlarging/contracting means 305 are related 
to the degree of correction of the parameters and to correct 
the parameters with reference to the table on the basis of the 
picture element density, the magnification and the target 
frequency component input from the input means . 

5 Further though, in the seventh embodiment, the frequency 

response characteristic of the 5-line processed image signal 
Sproc^ is obtained on the basis of the filtering characteristics 
Fl and F2 shown in Figures 58 and 59 and the parameters of the 
transformation function for transforming the band-limited 

0 signal corresponding to the target frequency component are 

corrected so that frequency response of the processed image 
signal Sp^,^^ obtained from the 5-line transformed original image 
signal S^^^5 becomes equal to that of the base processed image 
signal Sp^^^B, when the filtering characteristics Fl and F2 are 

5 represented by formulae, the frequency response characteristic 
of the 5-line processed image signal Sp^^^^S may be obtained on 



136 



the basis of the formulae. 

Further^ though, in the seventh embodiment, the base 
original image signal S^^^g' is transformed to the picture element 
density- trans formed original image signal S^^^gby filtering the 

5 base original image signal S^^g' by use of a picture element 
density-transformation filter, the base original image signal 
S^^g' may be transformed to the picture element density- 
transformed original image signal S^^g by other various methods. 
For example, linear interpolation or spline interpolation or 

0 thinning of picture elements after filtering may be employed. 
In such a case, the frequency response characteristic of the 
picture element density-transformed original image signal S^^g 
is calculated on the basis of the coefficient of contraction 
by the interpolation, or by Fourier transformation of the picture 

5 element density-transformed original image signal S^^g, and the 

parameters of the base transformation functions f are corrected 
on the basis of the frequency response characteristic. 

Also when a processed image signal Sp^.^^ is obtained by 
enlarging an intermediate processed image signal S^^^^' obtained 

0 from the base original image signal S^^^g' as it is by the 

enlarging/contracting means 305, or when a processed image 
signal Sp^^^ is obtained from a picture element density- 
transformed original image signal S^^^ without enlarging the 
intermediate processed image signal Sp^^^', the processed image 

5 signal Sp^^^can differ in frequency response characteristic from 

the base processed image signal S^^^^B. In such a case, the base 
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transformation functions t^^vaay be corrected on the basis of one 

of the filtering characteristics Fl and F2* 

Further, though, in the seventh embodiment described above, 

the unsharp image signals are obtained from the original image 
5 signal by filtering and interpolation/enlargement and the 

band-limited signals are obtained from the original image signal 

and the unsharp image signals, the band-limited signals may be 

made, for instance, by transforming the original image signal 

to multiple resolution image signals by a wavelet transformation 
0 or a Laplacian pyramid, and the image signals at the respective 

resolutions may be taken as the band-limited signals. 

An image processing system 331 in accordance with an eighth 

embodiment of the present invention will be described, 

hereinbelow. 

5 In Figure 65, the image processing system 331 in accordance 

with the eighth embodiment of the present invention comprises 
an unsharp image signal making means 332 which makes a plurality 
of unsharp image signals from multiple resolution image signals 
S^j.g'' which are obtained from a base original image signal S^^g' 

0 by transforming the base original image signal S^rg' to multiple 
resolution space and encoding the transformed image signals, 
and a frequency enhancement processing means 333 which carries 
out a frequency enhancement processing for enhancing a 
particular frequency and obtains an intermediate processed 

5 image signal Sp^.^^'. The image processing system 331 further 

comprises an enlarging/contracting means 335 which enlarges the 
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intermediate processed image signal Sp^^^', and a parameter 
correcting means 336. The processing performed by the unsharp 
image signal making means 332, the frequency enhancement 
processing means 333 and the enlarging/contracting means 335 

5 is the same as that in the seventh embodiment, and accordingly 
will not be described here. In this embodiment, the parameter 
correcting means 33 6 corrects the parameters of the base 
transformation functions fj.^ according to the characteristics 
of the wavelet transformation functions which are employed when 

0 transforming the original image signal S^^^^ to multiple 
resolution space. 

The base original image signal S^^^' is transformed into 
multiple resolution space and the transformed image signals are 
encoded in the following manner. That is, as shown in Figure 

5 66A, the base original image signal S^^g' is subjected to wavelet 

transformation and is decomposed into four pieces of data LLl, 
HLO, LHO and HHO by resolution. The data LLl represents an image 
obtained by reducing the original image to 1/2 in both 
longitudinal and lateral directions. The data HLO, LHO and HHO 

0 represents a longitudinal edge, a lateral edge and an oblique 

edge. When the four pieces of data LLl, HLO, LHO and HHO are 
subjected to inverse wavelet transformation, the original image 
signal S^^g is obtained. Then the data LLl is further subjected 
to wavelet transformation and decomposed into four pieces of 

5 data LL2 , HLl, LHl and HHl as shown in Figure 66B. The data LL2 

represents an image obtained by reducing data LLl to 1/2 in both 
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longitudinal and lateral directions. The data HLl, LHl and HHl 
represents a longitudinal edge component, a lateral edge 
component and an oblique edge component of the data LLl . When 
the four pieces of data LL2, HLl, LHl and HHl are subjected to 

5 inverse wavelet transformation, a low resolution image signal 
having a resolution equal to 1/2 of that of the original image 
signal is obtained* Further wavelet transformation is repeated 
a desired number of times on data LL obtained each time the 
wavelet transformation is carried out, thereby obtaining a 

0 plurality of pieces of data which are different in resolution. 
Thereafter, data at each resolution is encoded as shown in Figure 
66C and the multiple resolution image signals S^^^' ' is obtained. 

By decoding only image signals up to a desired resolution 
in the multiple resolution image signals S^^^' ' and carrying out 

5 inverse wavelet transformation on the image signals, a low 

resolution image signal representing an image at a resolution 
of 1/2^ (k being the desired resolution) of that of the original 
image can be obtained. 

As described above, by restoring low resolution image 

0 signals up to a desired resolution in the multiple resolution 

image signals S^rg' ' 9 ^ low resolution image signal representing 
an image at a resolution of 1/2^ (k being the desired resolution) 
of that of the original image can be obtained. However, the image 
reproduced on the basis of the low resolution image signal 

5 becomes unsharp depending on wavelet transformation functions 
used in the wavelet transformation. Figure 67 shows frequency 
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response characteristics of the low resolution image signals 
which are restored up to a resolution equal to 1/2 of that of 
the base original image signal S^^g' and up to a resolution equal 
to 1/4 of that of the base original image signal S^^^' when the 

5 base original image signal S^^g' has been decomposed by use of 
wavelet transformation functions having first coefficients of 
filter, and Figure 68 shows the same when the base original image 
signal S^^g' has been decomposed by use of wavelet transformation 
functions having second coefficients of filter. As can be seen 

0 from Figures 65 and 66, as the response deteriorates, high 

frequency components deteriorate and the image becomes less 
sharp. In Figures 65 and 66, the frequency response of the base 
original image signal S^^^' is 1 over the entire frequency bands. 

The multiple resolution image signals S^rg' ' obtained by 

5 decomposing a base original image signal S^^g' wavelet 

transformation are attached with information on the size of the 
wavelet transformation functions and/or the coefficients of 
filter thereof since the wavelet transformation functions to 
be employed in the wavelet transformation should be in one to 

0 one correspondence to the wavelet transformation functions to 

be employed in the inverse wavelet trans formation . Accordingly , 
in the eighth embodiment, when a low resolution image signal 
at a desired resolution restored from multiple resolution image 
signals S^^-g' ' obtained by decomposing a base original image 

5 signal S^^g' is taken as an object original image signal S^^g, the 

difference in response at a target frequency component between 
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the processed image signal Sp^^^ which is obtained by carrying 
out the frequency enhancement processing on the object original 
image signal S^^^ by use of the base transformation functions f^t 
and the base processed image signal S^^^^B obtained by carrying 
out the frequency enhancement processing on the base original 
image signal S^^g' by use of the base transformation functions 
t^^ is obtained on the basis of the information on the wavelet 
transformation functions attached to the multiple resolution 
image signals S^^g' ' and the parameters of the transformation 
functions for transforming the object original image signal S^^^ 
are corrected by the parameter correcting means 336 so that the 
difference is compensated for* When information on the wavelet 
transformation functions is not attached to the multiple 
resolution image signals S^^^'', the information is manually 
input into the image processing system 331, 

For example, when a low resolution image signal restored 
up to a resolution equal to 1/4 of that of a base original image 
signal S^^^ from multiple resolution image signals S^^g' ' obtained 
by decomposing the base original image signal S^^g' is taken as 
an object original image signal S^^g, aspect of the image 
reproduced on the basis of the 1/4 resolution image signal at 
a frequency of 1 cycle/mm is equalized to that of the image 
reproduced on the basis of the base original image signal S^^g' 
in the following manner. The frequency at which aspect of the 
image reproduced on the basis of the 1/4 resolution image signal 
is to be equalized to that of the image reproduced on the basis 



of the base original image signal S^^/ is input in advance through 
an input means by the user. Since the response of the image 
reproduced on the basis of the 1/4 resolution image signal at 
1 cycle/mm is 0.9 whereas the response of the image reproduced 
5 on the basis of the base original image signal S„g' at 1 cycle/mm 
is 1, it is necessary to multiply the frequency response at 1 
cycle/mm of the 1/4 resolution image signal by 1/0.9 in order 
to equalize aspect of the image reproduced on the basis of the 
1/4 resolution image signal to that of the image reproduced on 
C|o the basis of the base original image signal S^^/ at 1 cycle/mm. 

=4= Accordingly, in this embodiment, the transformation function 

W for the band-limited signal in the frequency band corresponding 

«3 to 1 cycle/mm is set by correcting the parameters for the base 

transformation function f^b for the band-limited signal in the 
Hs frequency band corresponding to 1 cycle/mm so that the degree 

4= of enhancement becomes 1/0.9 times . More specifically, when the 

O transformation function is a function having an inclination, 

the inclination is multiplied by 1/0.9, and when the 
transformation function is a constant, the value of the constant 
20 is multiplied by 1 /O . 9 . The frequency component at which aspect 

of the image reproduced on the basis of the 1/4 resolution image 
signal is to be equalized to that of the image reproduced on 
the basis of the base original image signal S^^^' may be any 
frequency and may be two or more frequencies. Further the 
25 frequency component may be input by the user through an input 

means not shown. 
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Thus in this embodiment, when frequency enhancement 
processing is to be carried out on an image signal which is lower 
in resolution than a base original image signal S^^^' and is 
obtained from multiple resolution image signals S^rg'' obtained 

5 by wavelet transformation of the base original image signal S^^^ ' , 
the transformation functions for transforming the band-limited 
signals obtained from the low resolution image signal are 
corrected on the basis of the wavelet transformation functions 
employed in the wavelet transformation of the base original image 

0 signal S^^g so that the image reproduced on the basis of the 
processed image signal obtained from the low resolution image 
signal becomes substantially equal to the image reproduced on 
the basis of the processed image signal obtained from the base 
original image signal in frequency response characteristic at 

5 a desired frequency. Accordingly, a processed image signal 
which is substantially equal in frequency response 
characteristic to a processed image signal obtained from the 
base original image signal can be constantly obtained 
irrespective of the resolution of the low resolution image signal . 

0 Further since the frequency response characteristic is 
corrected simultaneously with the frequency enhancement 
processing, the operation time can be saved as compared with 
when correction of the frequency response characteristic is 
carried out on the processed image signal Sp^^^ which has been 

5 subjected to the frequency enhancement processing. 

In the case where the processed image signal Sp^^^ is to be 
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enlarged by the enlarging/contracting means 335 in the eighth 
embodiment, the parameters of the base transformation functions 
f^b are corrected on the basis of the filtering characteristic 
of the filter employed by the enlarging/contracting means 335 
in addition to the characteristics of the wavelet transformation 
functions . 

In the seventh and eighth embodiments described above, 
it is preferred that the object original image signal S^^g and 
the parameters for the object original image signal S^^^ obtained 
by correcting the base transformation functions according to 
the characteristic of enlargement or the characteristic of 
contraction be stored correlated to each other. With this 
arrangement, when an image is output (output as a hard copy on 
photographic film or displayed on a CRT or the like) again on 
the basis of the object original image signal S^^g, an image at 
a quality equivalent to that of a previously output image can 
be obtained. 

Further, it is preferred that the base original image 
signal S^^g' from which the object original image signal S^^g is 
obtained, the characteristic of enlargement or the 
characteristic of contraction with which the base original image 
signal S^^^' is enlarged or contracted to the object original 
image signal S^^g and the parameters of the reference 
transformation functions be stored correlated to each other. 
With this arrangement, enlarged or contracted images can be 
repeatedly reproduced on the basis of the base original image 



signal S^^g' at the same image quality. 

The characteristic of enlargement or the characteristic 
of contraction may be represented by a space frequency 
characteristic and may be stored as table data such as shown 
5 in Figure 69. Otherwise the space frequency characteristic 
representing the characteristic of enlargement or the 
characteristic of contraction may be approximated by a function 
and parameters of the function may be stored as the 
characteristic of enlargement or the characteristic of 

Wo contraction. As such a function, a Gaussian function such as 

4; shown below may be employed. 

IzJ Res=exp(-a^/a^) 

wherein a=0 when Fq-m<0 and otherwise a=Fq-m/ Res represents 
1^^ the response and Fq represents the space frequency. In this case, 

fls a and m may be employed as the parameters representing the space 

4; frequency characteristic which in turn represents the 

characteristic of enlargement or the characteristic of 
contraction as shown in Figures 7 OA and 7 OB. 
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What is claimed is; 

!• An image processing method for obtaining a processed 
image signal from an original image signal representing an 
5 original image having a certain picture element density, in which 
a plurality of intermediate image signals which are different 
in frequency band are made on the basis of the original image 
signal, a plurality of transformed image signals are obtained 
by carrying out a transformation processing on the respective 
Ho intermediate image signals on the basis of respective 

transformation functions, and a processed image signal is 
IJf obtained from the transformed image signals, wherein the 

;f improvement comprises the step of 

defining said transformation functions by determining 
fl5 transformation function defining parameters for the 
3:: transformation functions on the basis of the picture element 

density of the original image. 

2. An image processing method as defined in Claim 1 in 
which said plurality of intermediate image signals are 
2 0 band-limited signals which are made by carrying out on the 
original image signal a filtering processing by use of filters 
whose coefficients of filter are determined on the basis of the 
picture element density of the original image signal, thereby 
making a plurality of unsharp image signals which are different 
25 in frequency response characteristic, and making a plurality 

of band-limited signals representing the signals in the 
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respective frequency bands of the original image signal on the 
basis of the unsharp image signals and the original image signal. 

3. An image processing method as defined in Claim 1 in 
which said predetermined transformation functions are non- 

5 linear functions ♦ 

4. An image processing method as defined in Claim 1 in 
which said transformation processing is a frequency enhancement 
processing. 

5. An image processing method as defined in Claim 1 in 
0 which said transformation processing is a dynamic range 

compression processing. 

6. An image processing method as defined in Claim 1 
further comprising the step of storing the transformation 
function defining parameters determined for the original image 

5 signal with the parameters related to the original image signal. 

7. An image processing method as defined in Claim 1 in 
which said step of defining said transformation functions by 
determining transformation function defining parameters for the 
transformation functions on the basis of the picture element 

0 density of the original image comprises the steps of preparing 
transformation function defining parameters for at least two 
reference picture element densities, comparing the picture 
element density of the original image with the reference picture 
element densities, and determining the transformation function 

5 defining parameters for one of the reference picture element 
densities closest to the picture element density of the original 
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image as the transformation function defining parameters for 
the original image signal. 

8. An image processing system for obtaining a processed 
image signal from an original image signal representing an 
original image having a certain picture element density 
comprising an intermediate image signal making means which makes 
a plurality of intermediate image signals which are different 
in frequency band on the basis of the original image signal, 
and a transformation means which obtains a plurality of 
transformed image signals by carrying out a transformation 
processing on the respective intermediate image signals on the 
basis of respective transformation functions, and obtains a 
processed image signal from the transformed image signals, 
wherein the improvement comprises that 

a transformation function defining means which determines 
transformation function defining parameters for the 
transformation functions on the basis of the picture element 
density of the original image and defines the transformation 
functions . 

9. An image processing system as defined in Claim 8 in 
which the intermediate image signal making means comprises an 
unsharp image signal making means which makes a plurality of 
unsharp image signals which are different in frequency response 
characteristic by carrying out on the original image signal a 
filtering processing by use of filters whose coefficients of 
filter are determined on the basis of the picture element density 
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of the original image, and a band-limited signal making means 
which makes, as the intermediates image signals, a plurality 
of band-limited signals representing the signals in the 
respective frequency bands of the original image signal on the 
basis of the unsharp image signals and the original image signal • 

10. An image processing system as defined in Claim 8 in 
which said predetermined transformation functions are non- 
linear functions. 

11. An image processing system as defined in Claim 8 in 
which said transformation processing is a frequency enhancement 
processing. 

12. An image processing system as defined in Claim 8 in 
which said transformation processing is a dynamic range 
compression processing. 

13. An image processing system as defined in Claim 8 
further comprising a means for storing the transformation 
function defining parameters determined for the original image 
signal with the parameters related to the original image signal. 

14. An image processing system as defined in Claim 8 in 
which said transformation function determining means stores 
therein transformation function defining parameters for at 
least two reference picture element densities, compares the 
picture element density of the original image represented by 
the original image signal with the reference picture element 
densities, and determines the transformation function defining 
parameters for one of the reference picture element densities 
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closest to the picture element density of the original image 
as the transformation function defining parameters for the 
original image signal. 

15. A computer-readable recording medium loaded with 
program for causing a computer to perform an image processing 
method for obtaining a processed image signal from an original 
image signal representing an original image having a certain 
picture element density in which a plurality of intermediate 
image signals which are different in frequency band are made 
on the basis of the original image, a plurality of transformed 
image signals are obtained by carrying out a transformation 
processing on the respective intermediate image signals on the 
basis of respective transformation functions, and a processed 
image signal is obtained from the transformed image signals, 
wherein the improvement comprises that 

the program includes the step of defining said 
transformation functions by determining transformation 
function defining parameters for the transformation functions 
on the basis of the picture element density of the original image. 

16. A computer-readable recording medium as defined in 
Claim 15 in which the step of making said plurality of 
intermediate image signals comprises the steps of carrying out 
on the original image signal a filtering processing by use of 
filters whose coefficients of filter are determined on the basis 
of the picture element density of the original image, thereby 
making a plurality of unsharp image signals which are different 



in frequency response characteristic, and making a plurality 
of band-limited signals representing the signals in the 
respective frequency bands of the original image signal, as said 
intermediate image signals, on the basis of the unsharp image 
5 signals and the original image signal. 

17. A computer-readable recording medium as defined in 
Claim 15 in which said predetermined transformation functions 
are non-linear functions. 

18. A computer-readable recording medium as defined in 
0 Claim 15 in which said transformation processing is a frequency 

enhancement processing. 

19. A computer-readable recording medium as defined in 
Claim 15 in which said transformation processing is a dynamic 
range compression processing. 

5 2 0, A computer-readable recording medium as defined in 

Claim 15 in which the program further includes the step of storing 
the transformation function defining parameters determined for 
the original image signal with the parameters related to the 
original image signal. 

0 21. A computer-readable recording medium as defined in 

Claim 15 in which said step of defining said transformation 
functions by determining transformation function defining 
parameters for the transformation functions on the basis of the 
picture element density of the original image comprises the steps 

5 of preparing transformation function defining parameters for 
at least two reference picture element densities, comparing the 
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picture element density of the original image represented by 
the original image signal with the reference picture element 
densities , and determining the transformation function defining 
parameters for one of the reference picture element densities 
closest to the picture element density of the original image 
signal as the transformation function defining parameters for 
the original image signal. 

22. An image processing method for obtaining a processed 
image signal from an original image signal representing an 
original image having a certain resolution, in which a plurality 
of intermediate image signals which are different in frequency 
band are made on the basis of the original image signal, a 
plurality of transformed image signals are obtained by carrying 
out a transformation processing on the respective intermediate 
image signals on the basis of respective transformation 
functions, and a processed image signal is obtained from the 
transformed image signals, wherein the improvement comprises 
that 

the transformation functions for processing the 
intermediate image signals obtained from an object original 
image signal to be processed are determined by correcting, 
according to the resolution of the object original image signal, 
reference transformation functions for transforming the 
intermediate image signals obtained from a reference original 
image signal representing an image having a reference resolution, 
the reference transformation functions having been determined 
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and stored in advance. 

23. An image processing method as defined in Claim 22 in 
which said plurality of intermediate image signals are 
band-limited signals which are made by carrying out on the 

5 original image signal a filtering processing by use of filters 
having predetermined coefficients of filter, thereby making a 
plurality of unsharp image signals which are different in 
frequency response characteristic, and making a plurality of 
band-limited signals representing the signals in the respective 

0 frequency bands of the original image signal on the basis of 

the unsharp image signals and the original image signal. 

24. An image processing method as defined in Claim 22 in 
which said transformation functions are non-linear functions. 

25. An image processing method as defined in Claim 22 in 
5 which said transformation processing is a frequency enhancement 

processing. 

26. An image processing method as defined in Claim 22 in 
which said transformation processing is a dynamic range 
compression processing. 

0 27. An image processing method as defined in Claim 22 in 

which the determined transformation functions are stored with 
the functions related to the original image signal. 

28. An image processing method as defined in Claim 22 in 
which the transformation functions for band-limited signals for 

5 the original image signal to be processed are calculated so that 
the frequency response characteristic of the processed image 
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signal Sp^^c conforms to that of the reference processed image 
signal Sp^^^O at least in frequency bands not lower than 1/5 of 
the Nyquist frequency of the original image signal to be 
processed. 

29. An image processing method as defined in Claim 28 in 
which the transformation functions for band-limited signals for 
the original image signal to be processed are calculated so that 
the frequency response characteristic of the processed image 
signal S^^^^ conforms to that of the reference processed image 
signal S^^^^O at least in frequency bands not lower than 1/2 of 
the Nyquist frequency of the original image signal to be 
processed. 

30. An image processing system comprising an 
intermediate image signal making means which makes a plurality 
of intermediate image signals which are different in frequency 
band on the basis of an original image signal representing an 
original image having a predetermined resolution, and a 
transformation means which obtains a plurality of transformed 
image signals by carrying out a transformation processing on 
the respective intermediate image signals on the basis of 
respective transformation functions, and obtains a processed 
image signal from the transformed image signals, wherein the 
improvement comprises that 

a transformation function calculating means which 
calculates the transformation functions for processing the 
intermediate image signals obtained from an object original 
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image signal to be processed by correcting, according to the 
resolution of the object original image signal, reference 
transformation functions for transforming the intermediate 
image signals obtained from a reference original image signal 
5 having a reference resolution, the reference transformation 
functions having been determined and stored in advance. 

31. An image processing system as defined in Claim 30 in 
which the intermediate image signal making means comprises an 
unsharp image signal making means which makes a plurality of 

©) unsharp image signals which are different in frequency response 

4= characteristic by carrying out on the original image signal a 

[y filtering processing by use of filters whose coefficients of 

filter are determined on the basis of the resolution of the 
^ original image signal, and a band-limited signal making means 

its which makes, as the intermediates image signals, a plurality 

of band-limited signals representing the signals in the 
O respective frequency bands of the original image signal on the 

basis of the unsharp image signals and the original image signal. 

32. An image processing system as defined in Claim 30 in 
20 which said predetermined transformation functions are non- 
linear functions. 

33. An image processing system as defined in Claim 30 in 
which said transformation processing is a frequency enhancement 
processing. 

25 34. An image processing system as defined in Claim 30 in 

which said transformation processing is a dynamic range 
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compression processing* 

35* An image processing system as defined in Claim 30 
further comprising a means for storing the determined 
transformation functions with the functions related to the 
original image signal. 

36. An image processing system as defined in Claim 30 in 
which the transformation function calculating means calculates 
the transformation functions for processing the intermediate 
image signals obtained from an object original image signal to 
be processed by correcting^ according to the resolution of the 
object original image signal, reference transformation 
functions so that the frequency response characteristic of the 
processed image signal conforms to that of the reference 
processed image signal at least in frequency bands not lower 
than 1/5 of the Nyquist frequency of the original image signal 
to be processed. 

37. An image processing system as defined in Claim 36 in 
which the transformation function calculating means calculates 
the transformation functions for processing the intermediate 
image signals obtained from an object original image signal to 
be processed by correcting, according to the resolution of the 
object original image signal, reference transformation 
functions so that the frequency response characteristic of the 
processed image signal conforms to that of the reference 
processed image signal at least in frequency bands not lower 
than 1/2 of the Nyquist frequency of the original image signal 



to be processed. 

38. A computer-readable recording medium loaded with 
program for causing a computer to perform an image processing 
method in which a plurality of intermediate image signals which 
are different in frequency band are made on the basis of an 
original image signal representing an original image having a 
predetermined resolution, a plurality of transformed image 
signals are obtained by carrying out a transformation processing 
on the respective intermediate image signals on the basis of 
respective transformation functions, and a processed image 
signal is obtained from the transformed image signals, wherein 
the improvement comprises that 

the program includes the step of calculating the 
transformation functions for processing the intermediate image 
signals obtained from an object original image signal to be 
processed by correcting, according to the resolution of the 
object original image signal, reference transformation 
functions for transforming the intermediate image signals 
obtained from a reference original image signal having a 
reference resolution, the reference transformation functions 
having been determined and stored in advance. 

39. A computer-readable recording medium as defined in 
Claim 3 8 in which the step of making said plurality of 
intermediate image signals comprises the steps of carrying out 
on the original image signal a filtering processing by use of 
filters whose coefficients of filter are determined on the basis 
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of the resolution of the original image signal, thereby making 
a plurality of unsharp image signals which are different in 
frequency response characteristic, and making a plurality of 
band-limited signals representing the signals in the respective 
5 frequency bands of the original image signal, as said 
intermediate image signals, on the basis of the unsharp image 
signals and the original image signal. 

40. A computer-readable recording medium as defined in 
Claim 3 8 in which said predetermined transformation functions 
10 are non-linear functions. 

" 41. A computer-readable recording medium as defined in 

J Claim 38 in which said transformation processing is a frequency 

enhancement processing. 

42. A computer-readable recording medium as defined in 
15 Claim 3 8 in which said transformation processing is a dynamic 
- range compression processing. 

3 43. A computer-readable recording medium as defined in 

Claim 3 8 in which said program further comprises the step of 
storing the determined transformation functions with the 

20 functions related to the original image signal. 

44. A computer-readable recording medium as defined in 
Claim 3 8 in which the reference transformation functions are 
corrected according to the resolution of the object original 
image signal so that the frequency response characteristic of 

25 the processed image signal conforms to that of the reference 

processed image signal at least in frequency bands not lower 



159 



than 1/5 of the Nyquist frequency of the original image signal 
to be processed. 

45. A computer-readable recording medium as defined in 
Claim 44 in which the reference transformation functions are 
5 corrected according to the resolution of the object original 
image signal so that the frequency response characteristic of 
the processed image signal Sp^-^^ conforms to that of the reference 
processed image signal Sp^^^O at least in frequency bands not lower 
than 1/2 of the Nyquist frequency of the original image signal 
iJQ to be processed. 

'%t 46. An image processing method for obtaining a processed 

image signal from an original image signal representing an 
original image having a certain resolution, in which a plurality 
of intermediate image signals which are different in frequency 

1^ band are made on the basis of the original image signal, a 
'J. plurality of transformed image signals are obtained by carrying 
out a transformation processing on the respective intermediate 
image signals on the basis of transformation functions which 
are set according to the frequency bands of the respective 

20 intermediate image signals, and a processed image signal is 

obtained from the transformed image signals, wherein the 
improvement comprises the steps of 

preparing reference transformation functions for 
transforming the intermediate image signals obtained from a 

25 reference original image signal representing a reference 

original image having a reference resolution which are set 
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according to the frequency bands of the respective intermediate 
image signals, and, when an object original image signal to be 
processed represents an original image having a resolution lower 
than the reference resolution, setting the transformation 
5 functions for transforming the intermediate image signals 
obtained from the object original image signal to be equal to 
the respective reference transformation functions for the 
intermediate image signals in the frequency bands not higher 
than the frequency band corresponding to the resolution of the 
original image represented by the object original image signal. 
4= 47. An image processing method as defined in Claim 46 in 

which said plurality of intermediate image signals are 
^ band-limited signals which are made by making, on the basis of 

: the original image signal, a plurality of unsharp image signals 

which are different in frequency response characteristic, and 
making a plurality of band-limited signals representing the 
signals in the respective frequency bands of the original image 
signal on the basis of the unsharp image signals and the original 
image signal. 

48. An image processing method as defined in Claim 46 in 
which information on the resolution of the object original image 
signal is obtained and the step of setting the transformation 
functions for transforming the intermediate image signals 
obtained from the object original image signal is executed on 

25 the basis of the information. 

49. An image processing system for obtaining a processed 



20 
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image signal from an original image signal representing an 
original image having a certain resolution comprising an 
intermediate signal making means which makes a plurality of 
intermediate image signals which are different in frequency band 
on the basis of the original image signal, and a transformation 
processing means which makes a plurality of transformed image 
signals by carrying out a transformation processing on the 
respective intermediate image signals on the basis of 
transformation functions which are set according to the 
frequency bands of the respective intermediate image signals 
and obtains a processed image signal from the transformed image 
signals, wherein the improvement comprises that 

the transformation processing means prepares reference 
transformation functions for transforming the intermediate 
image signals obtained from a reference original image signal 
representing a reference original image having a reference 
resolution which are set according to the frequency bands of 
the respective intermediate image signals, and, when the object 
original image signal to be processed represents an original 
image having a resolution lower than the reference resolution, 
sets the transformation functions for transforming the 
intermediate image signals obtained from the object original 
image signal to be equal to the respective reference 
transformation functions for the intermediate image signals in 
the frequency bands not higher than the frequency band 
corresponding to the resolution of the original image 



represented by the object original image signal. 

50. An image processing system as defined in Claim 49 in 
which the intermediate image signal making means comprises an 
unsharp image signal making means which makes a plurality of 
unsharp image signals which are different in frequency response 
characteristic on the basis of the original image, and a 
band-limited signal making means which makes, as the 
intermediates image signals, a plurality of band-limited 
signals representing the signals in the respective frequency 
bands of the original image signal on the basis of the unsharp 
image signals and the original image signal. 

51. An image processing system as defined in Claim 49 in 
which a resolution information obtaining means for obtaining 
information on the resolution of the original image represented 
by the object original image signal is provided and the 
transformation processing means sets the transformation 
functions for transforming the intermediate image signals 
obtained from the object original image signal on the basis of 
the information. 

52. A computer-readable recording medium loaded with 
program for causing a computer to perform an image processing 
method for obtaining a processed image signal from an original 
image signal representing an original image having a certain 
resolution, in which a plurality of intermediate image signals 
which are different in frequency band are made on the basis of 
the original image signal, a plurality of transformed image 
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signals are obtained by carrying out a transformation processing 
on the respective intermediate image signals on the basis of 
transformation functions which are set according to the 
frequency bands of the respective intermediate image signals, 
5 and a processed image signal is obtained from the transformed 

image signals, wherein the improvement comprises that 

the program includes the steps of preparing reference 
transformation functions for transforming the intermediate 
image signals obtained from a reference original image signal 
C3.0 representing a reference original image having a reference 

4^ resolution which are set according to the frequency bands of 

rlj the respective intermediate image signals, and, when the object 

€i original image signal to be processed represents an original 

3 image having a resolution lower than the reference resolution, 

K5 setting the transformation functions for transforming the 

JS intermediate image signals obtained from the object original 

O image signal to be equal to the respective reference 

transformation functions for the intermediate image signals in 
the frequency bands not higher than the frequency band 
20 corresponding to the resolution of the original image 

represented by the object original image signal. 

53. A computer-readable recording medium as defined in 
Claim 52 in which said plurality of intermediate image signals 
are band-limited signals which are made by making, on the basis 
25 of the original image signal, a plurality of unsharp image 

signals which are different in frequency response 
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characteristic, and making a plurality of band-limited signals 
representing the signals in the respective frequency bands of 
the original image signal on the basis of the unsharp image 
signals and the original image signal. 
5 54. A computer-readable recording medium as defined in 

Claim 52 in which information on the resolution of the object 
original image signal is obtained and the step of setting the 
transformation functions for transforming the intermediate 
image signals obtained from the object original image signal 

10 is executed on the basis of the information. 

55. An image processing method for making a 
contraction-processed image signal representing an image having 
a desired resolution lower than that of an image represented 
by a reference processed image signal which is obtained by making 

15 first to n-th reference low resolution image signals from a 

reference original image signal representing a reference 
original image, the reference low resolution image signals being 
different from each other in frequency band and respectively 
representing images whose resolutions are 1/2*" (k=l to n) of the 

20 resolution of the reference original image signal; making 

reference low-resolution band-limited signals on the basis of 
the reference low resolution image signals; obtaining a 
plurality of reference transformed image signals by carrying 
out a predetermined transformation processing on the reference 

25 low resolution band-limited signals on the basis of a plurality 

of reference transformation functions which are set according 
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to the respective frequency bands of the reference low resolution 
band-limited signals; and carrying out a predetermined 
processing on the reference transformed image signals, the image 
processing method characterized by the steps of 
5 taking as an original image signal one of said first to 

n-th reference low resolution image signals representing an 
image whose resolution is the closest to said desired resolution 
in the images represented by the reference low resolution image 
signals; making first to m-th low resolution image signals from 
^~flO the original image signal, the low resolution image signals being 

2r different from each other in frequency band and respectively 

^ representing images whose resolutions are 1/2^' (k'=l to m) of 

S the resolution of the original image signal; making a low- 

1^ resolution band-limited signals on the basis of the low 

FT 15 resolution image signals; obtaining a plurality of transformed 

21 image signals by carrying out said predetermined transformation 

processing on the low resolution band-limited signals on the 
basis of a part of said reference transformation functions which 
are set for the respective frequency bands of the reference low 
20 resolution band-limited signals not higher than the frequency 

band corresponding to the resolution closest to said desired 
resolution; carrying out said predetermined processing on the 
transformed image signals, thereby obtaining an intermediate 
processed image signal; enlarging or contracting the 
25 intermediate processed image so that the resolution of the image 

represented by the intermediate processed image signal becomes 
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equal to the desired resolution; and taking the enlarged or 
contracted intermediate processed image signal as the 
contraction-processed image signal. 

56. An image processing method as defined in Claim 55 in 
5 which said low resolution band-limited signals are made by making, 

on the basis of the original image signal, a plurality of unsharp 
image signals which are different in frequency response 
characteristic, and making a plurality of band-limited signals 
representing the signals in the respective frequency bands of 
0 the original image signal on the basis of the unsharp image 

signals and the original image signal. 

57. An image processing system for making a 
contraction-processed image signal representing an image having 
a desired resolution lower than that of an image represented 

5 by a reference processed image signal which is obtained by making 

first to n-th reference low resolution image signals from a 
reference original image signal representing a reference 
original image, the reference low resolution image signals being 
different from each other in frequency band and respectively 

0 representing images whose resolutions are 1/2^ (k=l to n) of the 

resolution of the original image represented by the reference 
original image signal; making reference low-resolution 
band-limited signals on the basis of the reference low resolution 
image signals; obtaining a plurality of reference transformed 

5 image signals by carrying out a predetermined transformation 

processing on the reference low resolution band-limited signals 
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on the basis of a plurality of reference transformation functions 
which are set according to the respective frequency bands of 
the reference low resolution band-limited signals; and carrying 
out a predetermined processing on the reference transformed 
5 image signals, the image processing system characterized by 

having 

a selection means which selects, as an original image 
signal, one of said first to n-th reference low resolution image 
signals representing an image whose resolution is the closest 

10 to said desired resolution in the images represented by the 

reference low resolution image signals; a low resolution image 
signal making means which makes first to m-th low resolution 
image signals from the original image signal, the low resolution 
image signals being different from each other in frequency band 

15 and respectively representing images whose resolutions are 1/2^' 

(k'=l to m) of the resolution of the original image signal; a 
low resolution band-limited signal making means which makes a 
low resolution band-limited signals on the basis of the low 
resolution image signals; a transformation processing means 

20 which obtains a plurality of transformed image signals by 
carrying out said predetermined transformation processing on 
the low resolution band-limited signals on the basis of a part 
of said reference transformation functions which are set for 
the respective frequency bands of the reference low resolution 

25 image signals not higher than the frequency band corresponding 

to the resolution closest to said desired resolution and carries 
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out said predetermined processing on the transformed image 
signals, thereby obtaining an intermediate processed image 
signal; and an enlarging/contracting means which enlarges or 
contracts the intermediate processed image so that the 
5 resolution of the image represented by the intermediate 
processed image signal becomes equal to the desired resolution. 

58. An image processing system as defined in Claim 57 in 
which said low resolution band-limited signal making means makes 
on the basis of the original image signal, a plurality of unsharp 
£|L0 image signals which are different in frequency response 

J; characteristic, and makes a plurality of band-limited signals 

representing the signals in the respective frequency bands of 
the original image signal on the basis of the unsharp image 
signals and the original image signal. 
f"15 59. A computer-readable recording medium loaded with 

|: program for causing a computer to perform an image processing 

method for making a contraction-processed image signal 
representing an image having a desired resolution lower than 
that of an image represented by a reference processed image 
2 0 signal which is obtained by making first to n-th reference low 

resolution image signals from a reference original image signal 
representing a reference original image, the reference low 
resolution image signals being different from each other in 
frequency band and respectively representing images whose 
2 5 resolutions are 1/2^ (k=l to n) of the resolution of the reference 

original image signal; making reference low-resolution 
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band-limited signals on the basis of the reference low resolution 
image signals; obtaining a plurality of reference transformed 
image signals by carrying out a predetermined transformation 
processing on the reference low resolution band-limited signals 
5 on the basis of a plurality of reference transformation functions 

which are set according to the respective frequency bands of 
the reference low resolution band-limited signals; and carrying 
out a predetermined processing on the reference transformed 
image signals, the improvement comprises that 
cjLO the program includes the steps of taking as an original 

image signal one of said first to n-th reference low resolution 
riJ image signals representing an image whose resolution is the 

J3 closest to said desired resolution in the images represented 

~ by the reference low resolution image signals; making first to 

Ml5 m-th low resolution image signals from the original image signal, 

4Z the low resolution image signals being different from each other 

□ in frequency band and respectively representing images whose 

resolutions are 1/2^' (k'=l tom) of the resolution of the original 
image signal; making a low-resolution band-limited signals on 
2 0 the basis of the low resolution image signals; obtaining a 

plurality of transformed image signals by carrying out said 
predetermined transformation processing on the low resolution 
band-limited signals on the basis of a part of said reference 
transformation functions which are set for the respective 
25 frequency bands of the reference low resolution band-limited 

signals not higher than the frequency band corresponding to the 
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resolution closest to said desired resolution; carrying out said 
predetermined processing on the transformed image signals, 
thereby obtaining an intermediate processed image signal; 
enlarging or contracting the intermediate processed image so 
that the resolution of the image represented by the intermediate 
processed image signal becomes equal to the desired resolution; 
and taking the enlarged or contracted intermediate processed 
image signal as the contraction-processed image signal. 

60. A computer-readable recording medium as defined in 
Claim 59 in which said low resolution band-limited signals are 
made by making, on the basis of the original image signal, a 
plurality of unsharp image signals which are different in 
frequency response characteristic, and making a plurality of 
band-limited signals representing the signals in the respective 
frequency bands of the original image signal on the basis of 
the unsharp image signals and the original image signal. 

61. An image processing method for obtaining a processed 
image signal from an original image signal representing an 
original image having a certain resolution, in which a plurality 
of intermediate image signals which are different in frequency 
band are made on the basis of the original image signal, a 
plurality of transformed image signals are obtained by carrying 
out a transformation processing on the respective intermediate 
image signals on the basis of transformation functions which 
are set according to the frequency bands of the respective 
intermediate image signals, and a processed image signal is 



obtained from the transformed image signals, wherein the 
improvement comprises the steps of 

preparing reference transformation functions for 
transforming the intermediate image signals obtained from a 
reference original image signal representing a reference 
original image having a reference resolution which are set 
according to the frequency bands of the respective intermediate 
image signals, and, when an object original image signal to be 
processed represents an original image having a resolution 
higher than the reference resolution, setting the 
transformation functions for transforming the intermediate 
image signals obtained from the object original image signal 
in frequency bands not higher than the frequency band 
corresponding to the reference resolution to be equal to the 
respective reference transformation functions and the 
transformation functions for transforming the intermediate 
image signals obtained from the object original image signal 
in frequency bands higher than the frequency band corresponding 
to the reference resolution to be equal to the reference 
transformation function for transforming the intermediate image 
signal obtained from the reference original image signal in the 
highest frequency band. 

62. An image processing method as defined in Claim 61 in 
which said plurality of intermediate image signals are 
band-limited signals which are made by making, on the basis of 
the original image signal, a plurality of unsharp image signals 



which are different in frequency response characteristic, and 
making a plurality of band-limited signals representing the 
signals in the respective frequency bands of the original image 
signal on the basis of the unsharp image signals and the original 
image signal. 

63. An image processing method as defined in Claim 61 in 
which information on the resolution of the object original image 
signal is obtained and the step of setting the transformation 
functions for transforming the intermediate image signals 
obtained from the object original image signal is executed on 
the basis of the information. 

64. An image processing method as defined in Claim 61 
further comprising the step of storing the transformation 
function defining parameters determined for the original image 
signal with the parameters related to the original image signal. 

65 . An image processing system for obtaining a processed 
image signal from an original image signal representing an 
original image having a certain resolution comprising an 
intermediate signal making means which makes a plurality of 
intermediate image signals which are different in frequency band 
on the basis of the original image signal, and a transformation 
processing means which makes a plurality of transformed image 
signals by carrying out a transformation processing on the 
respective intermediate image signals on the basis of 
transformation functions which are set according to the 
frequency bands of the respective intermediate image signals 
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and obtains a processed image signal from the transformed image 
signals, wherein the improvement comprises that 

the transformation processing means prepares reference 
transformation functions for transforming the intermediate 
image signals obtained from a reference original image signal 
representing a reference original image having a reference 
resolution which are set according to the frequency bands of 
the respective intermediate image signals, and, when an object 
original image signal to be processed represents an original 
image having a resolution higher than the reference resolution, 
sets the transformation functions for transforming the 
intermediate image signals obtained from the object original 
image signal in frequency bands not higher than the frequency 
band corresponding to the reference resolution to be equal to 
the respective reference transformation functions and the 
transformation functions for transforming the intermediate 
image signals obtained from the object original image signal 
in frequency bands higher than the frequency band corresponding 
to the reference resolution to be equal to the reference 
transformation function for transforming the intermediate image 
signal obtained from the reference original image signal in the 
highest frequency band. 

66. An image processing system as defined in Claim 65 in 
which the intermediate image signal making means comprises an 
unsharp image signal making means which makes a plurality of 
unsharp image signals which are different in frequency response 



characteristic on the basis of the original image, and a 
band-limited signal making means which makes, as the 
intermediates image signals, a plurality of band-limited 
signals representing the signals in the respective frequency 
bands of the original image signal on the basis of the unsharp 
image signals and the original image signal. 

67. An image processing system as defined in Claim 65 in 
which a resolution information obtaining means for obtaining 
information on the resolution of the original image represented 
by the object original image signal is provided and the 
transformation processing means sets the transformation 
functions for transforming the intermediate image signals 
obtained from the object original image signal on the basis of 
the information. 

68. An image processing system as defined in Claim 65 
further comprising a means for storing the transformation 
function defining parameters determined for the original image 
signal with the parameters related to the original image signal. 

69. A computer-readable recording medium loaded with 
program for causing a computer to perform an image processing 
method for obtaining a processed image signal from an original 
image signal representing an original image having a certain 
resolution, in which a plurality of intermediate image signals 
which are different in frequency band are made on the basis of 
the original image signal, a plurality of transformed image 
signals are obtained by carrying out a transformation processing 
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on the respective intermediate image signals on the basis of 
transformation functions which are set according to the 
frequency bands of the respective intermediate image signals, 
and a processed image signal is obtained from the transformed 
image signals, wherein the improvement comprises that 

the program includes the steps of preparing reference 
transformation functions for transforming the intermediate 
image signals obtained from a reference original image signal 
representing a reference original image having a reference 
resolution which are set according to the frequency bands of 
the respective intermediate image signals, and, when an object 
original image signal to be processed represents an original 
image having a resolution higher than the reference resolution, 
setting the transformation functions for transforming the 
intermediate image signals obtained from the object original 
image signal in frequency bands not higher than the frequency 
band corresponding to the reference resolution to be equal to 
the respective reference transformation functions and the 
transformation functions for transforming the intermediate 
image signals obtained from the object original image signal 
in frequency bands higher than the frequency band corresponding 
to the reference resolution to be equal to the reference 
transformation function for transforming the intermediate image 
signal obtained from the reference original image signal in the 
highest frequency band. 

70. A computer-readable recording medium as defined in 



Claim 69 in which said plurality of intermediate image signals 
are band-limited signals which are made by making, on the basis 
of the original image signal, a plurality of unsharp image 
signals which are different in frequency response 
characteristic, and making a plurality of band-limited signals 
representing the signals in the respective frequency bands of 
the original image signal on the basis of the unsharp image 
signals and the original image signal. 

71. A computer-readable recording medium as defined in 
Claim 69 in which information on the resolution of the object 
original image signal is obtained and the step of setting the 
transformation functions for transforming the intermediate 
image signals obtained from the object original image signal 
is executed on the basis of the information. 

72. A computer-readable recording medium as defined in 
Claim 69 in which the program further includes the step of storing 
the transformation function defining parameters determined for 
the original image signal with the parameters related to the 
original image signal. 

73 . An image processing method for obtaining a processed 
image signal from an original image signal representing an 
original image having a certain picture element density, in which 
a plurality of intermediate image signals which are different 
in frequency band are made on the basis of the original image 
signal, a plurality of transformed image signals are obtained 
by carrying out a transformation processing on the respective 



intermediate image signals on the basis of respective 
transformation functions, and a processed image signal in which 
a predetermined frequency component is enhanced is obtained from 
the transformed image signals, wherein the improvement 
comprises that 

when an object original image signal to be processed is 
a contracted image signal obtained by carrying out a picture 
element density transformation processing on a reference 
original image signal representing an image having a reference 
picture element density, the transformation functions for 
processing the intermediate image signals obtained from the 
object original image signal are determined by correcting, 
according to the characteristic of the picture element density 
transformation processing, reference transformation functions 
for transforming the intermediate image signals obtained from 
the reference original image so that the frequency response 
characteristic of a desired frequency component of the processed 
image signal obtained from the object original image signal 
becomes equivalent to that of the processed image signal obtained 
from the reference original image signal, the reference 
transformation functions having been determined and stored in 
advance. 

74. An image processing method as defined in Claim 73 in 
which when the processed image signal obtained from the 
contracted image signal is enlarged at a magnification, the 
transformation functions for processing the intermediate image 
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signals obtained from the object original image signal are 
determined by correcting the reference transformation functions 
according to both the characteristic of the picture element 
density transformation processing and the characteristic of the 
enlargement. 

75. An image processing method as defined in Claim 73 in 
which said plurality of intermediate image signals are 
band-limited signals which are made by making from the original 
image signal a plurality of unsharp image signals which are 
different in frequency response characteristic, and making a 
plurality of band-limited signals representing the signals in 
the respective frequency bands of the original image signal on 
the basis of the unsharp image signals and the original image 
signal . 

76. An image processing method as defined in Claim 73 in 
which the object original image signal and the parameters for 
the object original image signal obtained by correcting the 
reference transformation functions according to the 
characteristic of contraction are stored correlated to each 
other. 

77. An image processing method as defined in Claim 73 in 
which the reference original image signal from which the object 
original image signal is obtained, the characteristic of 
contraction with which the reference original image signal is 
contracted to the object original image signal and the parameters 
of the reference transformation functions be stored correlated 



to each other. 

78. An image processing method for obtaining a processed 
image signal from an original image signal representing an 
original image, in which a plurality of intermediate image 
signals which are different in frequency band are made on the 
basis of the original image signal, a plurality of transformed 
image signals are obtained by carrying out a transformation 
processing on the respective intermediate image signals on the 
basis of respective transformation functions, and a processed 
image signal in which a predetermined frequency component is 
enhanced is obtained from the transformed image signals, wherein 
the improvement comprises that 

when an enlarged processed image signal is to be obtained 
by enlarging the processed image signal at a desired 
magnification by carrying out an enlargement transformation on 
the processed image signal, the transformation functions for 
transforming the respective intermediate image signals are 
corrected according to the characteristic of the enlargement 
transformation so that the frequency response characteristic 
of a desired frequency component of the enlarged processed image 
signal becomes a predetermined frequency response 
characteristic . 

79. An image processing method as defined in Claim 7 8 in 
which said plurality of intermediate image signals are 
band-limited signals which are made by making from the original 
image signal a plurality of unsharp image signals which are 
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different in frequency response characteristic^ and making a 
plurality of band-limited signals representing the signals in 
the respective frequency bands of the original image signal on 
the basis of the unsharp image signals and the original image 
signal* 

80. An image processing method as defined in Claim 78 in 
which the object original image signal and the parameters for 
the object original image signal obtained by correcting the 
reference transformation functions according to the 
characteristic of enlargement are stored correlated to each 
other. 

81. An image processing method as defined in Claim 78 in 
which the reference original image signal from which the object 
original image signal is obtained, the characteristic of 
enlargement with which the reference original image signal is 
enlarged to the object original image signal and the parameters 
of the reference transformation functions be stored correlated 
to each other. 

82. An image processing system for obtaining a processed 
image signal from an original image signal representing an 
original image having a certain picture element density, in which 
a plurality of intermediate image signals which are different 
in frequency band are made on the basis of the original image 
signal, a plurality of transformed image signals are obtained 
by carrying out a transformation processing on the respective 
intermediate image signals on the basis of respective 
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transformation functions, and a processed image signal in which 
a predetermined frequency component is enhanced is obtained from 
the transformed image signals, wherein the improvement 
comprising 

a transformation function correcting means which, when 
an object original image signal to be processed is a contracted 
image signal obtained by carrying out a picture element density 
transformation processing on a reference original image signal 
representing an image having a reference picture element density, 
determines the transformation functions for processing the 
intermediate image signals obtained from the object original 
image signal by correcting, according to the characteristic of 
the picture element density transformation processing, 
reference transformation functions for transforming the 
intermediate image signals obtained from the reference original 
image so that the frequency response characteristic of a desired 
frequency component of the processed image signal obtained from 
the object original image signal becomes equivalent to that of 
the processed image signal obtained from the reference original 
image signal, the reference transformation functions having 
been determined and stored in advance. 

83. An image processing system as defined in Claim 82 in 
which, when the processed image signal obtained from the 
contracted image signal is enlarged at a magnification, the 
transformation function correcting means determines the 
transformation functions for processing the intermediate image 
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signals obtained from the object original image signal by 
correcting the reference transformation functions according to 
both the characteristic of the picture element density 
transformation processing and the characteristic of the 
enlargement. 

84. An image processing system as defined in Claim 82 in 
which said plurality of intermediate image signals are 
band-limited signals which are made by making from the original 
image signal a plurality of unsharp image signals which are 
different in frequency response characteristic, and making a 
plurality of band-limited signals representing the signals in 
the respective frequency bands of the original image signal on 
the basis of the unsharp image signals and the original image 
signal . 

85. An image processing method as defined in Claim 82 
further comprising a means for storing the object original image 
signal and the parameters for the object original image signal 
obtained by correcting the reference transformation functions 
according to the characteristic of contraction correlated to 
each other. 

86. An image processing method as defined in Claim 82 
further comprising a means for storing the reference original 
image signal from which the object original image signal is 
obtained, the characteristic of contraction with which the 
reference original image signal is contracted to the object 
original image signal and the parameters of the reference 
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transformation functions correlated to each other. 

87. An image processing system for obtaining a processed 
image signal from an original image signal representing an 
original image, in which a plurality of intermediate image 
signals which are different in frequency band are made on the 
basis of the original image signal, a plurality of transformed 
image signals are obtained by carrying out a transformation 
processing on the respective intermediate image signals on the 
basis of respective transformation functions, and a processed 
image signal in which a predetermined frequency component is 
enhanced is obtained from the transformed image signals , wherein 
the improvement comprising 

a transformation function correcting means which, when 
an enlarged processed image signal is to be obtained by enlarging 
the processed image signal at a desired magnification by carrying 
out an enlargement transformation on the processed image signal, 
corrects the transformation functions for transforming the 
respective intermediate image signals according to the 
characteristic of the enlargement transformation so that the 
frequency response characteristic of a desired frequency 
component of the enlarged processed image signal becomes a 
predetermined frequency response characteristic. 

88. An image processing system as defined in Claim 87 in 
which said plurality of intermediate image signals are 
band-limited signals which are made by making from the original 
image signal a plurality of unsharp image signals which are 
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different in frequency response characteristic, and making a 
plurality of band-limited signals representing the signals in 
the respective frequency bands of the original image signal on 
the basis of the unsharp image signals and the original image 
signal . 

89. An image processing method as defined in Claim 87 
further comprising a means for storing the object original image 
signal and the parameters for the object original image signal 
obtained by correcting the reference transformation functions 
according to the characteristic of enlargement correlated to 
each other* 

90. An image processing method as defined in Claim 87 
further comprising a means for storing the reference original 
image signal from which the object original image signal is 
obtained, the characteristic of enlargement with which the 
reference original image signal is enlarged to the object 
original image signal and the parameters of the reference 
transformation functions correlated to each other. 

91. A computer-readable recording medium loaded with 
program for causing a computer to perform an image processing 
method for obtaining a processed image signal from an original 
image signal representing an original image having a certain 
picture element density, in which a plurality of intermediate 
image signals which are different in frequency band are made 
on the basis of the original image signal, a plurality of 
transformed image signals are obtained by carrying out a 
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transformation processing on the respective intermediate image 
signals on the basis of respective transformation functions, 
and a processed image signal in which a predetermined frequency 
component is enhanced is obtained from the transformed image 
signals, wherein the improvement comprises that 

the program includes the step of, when an object original 
image signal to be processed is a contracted image signal 
obtained by carrying out a picture element density 
transformation processing on a reference original image signal 
representing an image having a reference picture element density, 
determining the transformation functions for processing the 
intermediate image signals obtained from the object original 
image signal by correcting, according to the characteristic of 
the picture element density transformation processing, 
reference transformation functions for transforming the 
intermediate image signals obtained from the reference original 
image so that the frequency response characteristic of a desired 
frequency component of the processed image signal obtained from 
the object original image signal becomes equivalent to that of 
the processed image signal obtained from the reference original 
image signal, the reference transformation functions having 
been determined and stored in advance. 

92. A computer-readable recording medium as defined in 
claim 91 in which said program further includes the step of, 
when the processed image signal obtained from the contracted 
image signal is enlarged at a magnification, determining the 
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transformation functions for processing the intermediate image 
signals obtained from the object original image signal by 
correcting the reference transformation functions according to 
both the characteristic of the picture element density 
transformation processing and the characteristic of the 
enlargement. 

93 . A computer-readable recording medium as defined in 
Claim 91 in which said plurality of intermediate image signals 
are band-limited signals which are made by making from the 
original image signal a plurality of unsharp image signals which 
are different in frequency response characteristic, and making 
a plurality of band-limited signals representing the signals 
in the respective frequency bands of the original image signal 
on the basis of the unsharp image signals and the original image 
signal. 

94. A computer-readable recording medium as defined in 
Claim 91 in which said program further includes the step of 
storing the object original image signal and the parameters for 
the object original image signal obtained by correcting the 
reference transformation functions according to the 
characteristic of contraction correlated to each other. 

95. A computer-readable recording medium as defined in 
Claim 91 in which said program further includes the step of 
storing the reference original image signal from which the object 
original image signal is obtained, the characteristic of 
contraction with which the reference original image signal is 
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contracted to the object original image signal and the parameters 
of the reference transformation functions correlated to each 
other. 

96. A computer-readable recording medium loaded with 
program for causing a computer to perform an image processing 
method for obtaining a processed image signal from an original 
image signal representing an original image, in which a plurality 
of intermediate image signals which are different in frequency 
band are made on the basis of the original image signal, a 
plurality of transformed image signals are obtained by carrying 
out a transformation processing on the respective intermediate 
image signals on the basis of respective transformation 
functions, and a processed image signal in which a predetermined 
frequency component is enhanced is obtained from the transformed 
image signals, wherein the improvement comprises that 

the program includes the step of, when an enlarged 
processed image signal is to be obtained by enlarging the 
processed image signal at a desired magnification by carrying 
out an enlargement transformation on the processed image signal, 
correcting the transformation functions for transforming the 
respective intermediate image signals according to the 
characteristic of the enlargement transformation so that the 
frequency response characteristic of a desired frequency 
component of the enlarged processed image signal becomes a 
predetermined frequency response characteristic. 

97. A computer-readable recording medium as defined in 
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Claim 96 in which said plurality of intermediate image signals 
are band-limited signals which are made by making from the 
original image signal a plurality of unsharp image signals which 
are different in frequency response characteristic, and making 
a plurality of band-limited signals representing the signals 
in the respective frequency bands of the original image signal 
on the basis of the unsharp image signals and the original image 
signal. 

98. A computer-readable recording medium as defined in 
Claim 9 6 in which said program further includes the step of 
storing the object original image signal and the parameters for 
the object original image signal obtained by correcting the 
reference transformation functions according to the 
characteristic of enlargement correlated to each other. 

99. A computer-readable recording medium as defined in 
Claim 96 in which said program further includes the step of 
storing the reference original image signal from which the object 
original image signal is obtained, the characteristic of 
enlargement with which the reference original image signal is 
enlarged to the object original image signal and the parameters 
of the reference transformation functions correlated to each 
other. 
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ABSTRACT OF THE DISCLOSURE 
In an image processing method for obtaining a processed 
image signal from an original image signal representing an 
original image having a certain picture element density, a 
plurality of intermediate image signals which are different in 
frequency band are made on the basis of the original image signal . 
A plurality of transformed image signals are obtained by carrying 
out a transformation processing on the respective intermediate 
image signals on the basis of respective transformation 
functions, and a processed image signal is obtained from the 
transformed image signals- The transformation functions are 
defined by determining transformation function defining 
parameters for the transformation functions on the basis of the 
picture element density of the original image. 
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Declaration and Power of Attorney for Patent Application 
Japanese Language Declaration 



C T IS ft ^ ^ - ^ ) ^ L < t ± ^ ^0 7i ^ o ^ di ^ 0^ # T! 



As a beiow named inventor, I hereby declare that: 

Masahiko Yamada 

My residence, post office address and citizenship are as 
stated next to my name, c/o Fuji Photo Film 
Co., Ltd., 798 Miyanodai, Kaisei-machi , 
Ashigarakami-gun , Kanagawa-ken , Japan 
I believe I am the original, first and sole inventor (if only one 
name is listed below) or an original, first and joint inventor (if 
plural names are listed below) of the subject matter which is 
claimed and for which a patent is sought on the invention 
entitled 



"METHOD OF AND SYSTEM FOR IMAGE 



PROCESSING AND RECORDING MEDIUM 



FOR CARRYING OUT THE METHOD" 



the specification of which is attached hereto unless the 
following box is checked: 



□ was filed on 



as United States Application Number or 
PCT International Application Number 



and was amended on 



^ (if applicable). 



I hereby state that I have reviewed and understand the 
contents of the above-identified specification, including the 
claims, as amended by any amendment referred to above. 



I acknowledge the dtrty to disclose information which is 
material to patentability as defined in Title 37, Code of 
Federal Regulations, Section 1 .56. 
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Japanese Language Declaration 



Prior Foreign Applications 

^mxco^ntt\m 



(patent) 7293/1999 
(Number) 

(#^) 

(pat ent) 15904/1999 
(Number) 

(pat ent) 17250/1999 
(Number) 

(paterft) 26323/1999 



Japan 



(Country) 
Japan 



(Country) 
Japan 



(Country) 



Japan 

mt. M2sm^mf^mii9^ie)m{zm-^\,^xTm<7:>^mi^nmm 



I hereby claim foreign pnority under Title 35, United States 
Code, Section 119(a)-(d) or 365(b) of any foreign app!ication(s) 
for patent or inventor's certificate, or 365(a) of any PCT 
Internationai application which designated at least one country 
other than the United States, listed below and have also 
identified below, by checking the box, any foreign application for 
patent or inventor's certificate, or PCT International application 
having a filing date before that of the application on which 
priority is claimed. 



14/01/1999 



Priority Not Claimed 
□ 



(Day/MonthA'ear Filed) 

25/01/1999 □ 

(Day/Month/Year Filed) 

26/01/1999 □ 

(Day/Month/Year Filed) 

03/02/1999 
I hereby claim the benefit under Title 35, United States Code, 
Section 119(e) of any United States provisional application(s) 
listed below. 



(Application No.) 



(Filing Date) 



56 Ji -c ^ s ^ n rc# ff 5S[ IS ^ M -r -sil ^ if ^ o V ^ 
xm^i^m^^^h^zt^mmLx^^-^-^-r, 



(Application No.) 



(Filing Date) 



1 hereby claim the benefit of Title 35. United States Code 
Section 120 of any United States application(s), or 365(c) of 
any PCT International application designating the United States, 
listed below and, insofar as the subject matter of each of the 
claims of this application is not disclosed in the prior United 
States or PCT Internationai application In the manner provided 
by the first paragraph of Title 35, United States Code Section 
112, I acknowledge the duty to disclose any material 
information which is material to patentability as defined in Title 
37, Code of Federal Regulations, Section 1 .56 which became 
available between the filing date of the prior application and the 
national or PCT International filing date of this application: 



(Application No.) (R(ing Date) 



(Application No.) (Filing Date) 

%\t. fi i ^<o^m^zm^^^x:^'g:m^xfi^^rro^m7!}^ 
mm ^^±xMnv^^tm\ix^^^zt. ^hi,c^M^zfj:^tifz 



(Status: Patented, Pending, Abandoned) 



(Status: Patented, Pending, Abandoned) 

I hereby declare that all statements made herein of my own 
knowledge are true and that ait statements made on Information 
and belief are believed to be true; and further that these 
statements were made with the knowledge that willful false 
statements and the like so made are punishable by fine or 
imprisonment, or both, under Section 1CX}1 of Title 18 of the 
United States Code and that such willful false statements may 
jeopardize the validity of the application or any patent issued 
thereon. 
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Japanese Language Declaration 



POWER OF ATTORNEY; As a named inventor, I hereby 
appoint the following attorney(s) and/or agent(s) to 
prosecute this application and transact all business in the 
Patent and Trademark Office connected therewith {list 
name and registration number) 



John H. Mion, Reg. No. 18,879; Thomas J. Macpeak, Reg. No. 19,292; Robert J. Seas, Jr., Reg No 21 092- 
Darryl Mexic, Reg. No. 23,063; Robert V. Sloan, Reg, No. 22,775; Peter D. Olexy, Reg. No, 24,513; J Frank 
Osha, Reg. No. 24,625; Waddell A. Biggart, Reg, No. 24,861; Louis Gubinsky, Reg. No. 24.835; Neil B 
Siegel, Reg. No. 25,200; David J. Gushing, Reg. No. 28,703; John R, Inge, Reg. No. 26,916; Joseph J. Ruch, 
Jr., Reg. No. 26,577; Sheldon I. Landsman, Reg. No, 25,430; Richard C, Turner, Reg. No. 29,710; Howard l' 
Bernstein, Reg. No. 25,665; Alan J. Kasper, Reg. No, 25,426; Kenneth J, Burchfiel, Reg. No. 31,333; Gordon 
Kit. Reg. No. 30,764; Susan J, Mack, Reg. No. 30,951; Frank L. Bernstein, Reg, No. 31,484; Mark Boland, 
Reg. No, 32,197; WiUiam H. Mandir, Reg. No. 32,156; Scott M. Daniels, Reg. No. 32,562; Brian W. Hannon 
Reg. No. 32,778; Abraham J. Rosner, Reg, No. 33,276; Bruce E. Kramer, Reg. No. 33,725; Paul F NeUs Reg 
No. 33, 102 and Brett S. Sylvester, Reg. No, 32,765 



'MM^j^H9c: Send Correspondence to: 

SUGHRUE, MION, ZINN, MACPEAK & SEAS, PLLC 
2100 Pennsylvania Avenue, N.W., Washington, D.C. 20037-3202 



mmmt^mm^: i^^Ramm^^r) Direct Telephone Calls to: (name and telephone number) 

I (202)293-7060 







Full name of sole or first inventor 
Masahiko Yamada 




0 ft 


inventor's signature Date _ 

/ / ^ Decembei 






Residence ' ^ . 

Kaisei-machi , Japan 


[mm 




Citizenship 

Japan 




Post office address Fuji Photo Film Co., 
Litd. , 798 Miyanodai, Kaisei-machi , 




?^shigarakami-gun, Kanagawa-ken , 
Japan 




Full nanne of second joint inventor, if any 




B ft 


Second inventor's signature Date 




Residence 


mm 


Citizenship 




Post office address 







im ^\;x^(D^m^m^ \cr>^-Xi,m m{zmmL. m^^i-^z {Supply similar information and signature for third and 
^ subsequent joint inventors.) 
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